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Several  problems  related  to  model  reference  adaptive  systems  (MRAS) 
are  investigated.  The  MRAS  finds  application  both  in  control  system  design 
in  an  uncertain  environment  and  in  parameter  identification.  The  development 
here  is  based  on  a parallel  MRAS  structure  whose  adaptation  design  guarantees 
that  the  MRAS  is  hyperstable. 
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For  MRAS  applied  to  parameter  identification,  we  consider  two  problems.  Necess- 
ary and  sufficient  conditions  guaranteeing  asymptotic  stability  of  the  param- 
eter-error are  found.  For  a specific  discrete-time  hyperstable  MRAS,  we  con- 
sider the  difficulty  in  realizing  the  hyperstable  design  due  to  the  nature  of 
the  identification  problem.  An  alternate  design  is  proposed  and  a strategy 
for  investigating  its  stability  properties  is  described.  For  a continuous- 
time MRAS  control  system  with  discontinuous  adaptation,  variable  structure 
system  theory  is  used  to  establish  its  stability. 

A new  MRAS  structure  containing  multiple  reference  models  is  intro- 
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CHAPTER  1 


INTRODUCTION 

1.1  The  Model-Reference  Adaptive  System 

The  field  of  model-reference  adaptive  systems  (MRAS)  has  developed 
during  the  1960's  as  an  approach  to  the  design  of  control  systems.  It  pro- 
vides an  alternative  to  the  optimal  control  approach,  a viewpoint  which  has 
received  the  predominant  attention  from  control  system  researchers  during 
the  past  decade.  The  term  model  reference  adaptive  control  system  (MRACS) 
refers  to  a variety  of  system  structures  which  have  the  following  common 
features.  The  desired  control  system  performance  is  specified  by  the  state 
or  output  response  of  a MRACS  subsystem  which  we  denote  the  M-subsystem. 

The  control  system  itself  is  part  of  the  MRACS  structure,  and  we  denote  it 
as  the  S-subsystem.  It  is  designed  to  accept  either  on-line  changes  to  its 
controller  parameters  or  an  augmented  input  signal.  The  S-subsystem  parameter 
adaptation  or  augmented  input  signal  is  generated  by  a third  subsystem  which 
we  denote  the  A-subsystem.  This  subsystem  generates  the  parameter  adjustment 
or  augmented  input  signal  based  on  a comparison  between  the  output  or  state 
responses  of  the  M-  and  S-subsystems.  The  overall  objective  of  the  MRACS 
is  for  the  S-subsystem  output  or  state  trajectory  to  converge  to  the  corres- 
ponding M-subsystem  trajectory. 

In  this  thesis,  attention  will  be  focused  on  the  so-called  "parallel" 
MRAS  structure,  depicted  in  Figure  1.1-1.  The  name  has  resulted  from  the 
parallel  connection  between  the  M-  and  S-subsystems. 


In  addition  to  the  control  system  context,  the  field  of  MRAS  may 
be  directed  toward  the  system  identification  problem.  In  this  case,  the 
M-subsystem  corresponds  to  the  (possibly  controlled)  plant  while  the  S-subsystem 
corresponds  to  an  identification  model  with  adjustable  parameters. 

Motivation  for  application  of  the  MRAS  concept  to  control  system 
design  is  based  on  the  realization  that  the  plant  model  used  in  the  design 
of  a control  system  (either  by  optimal  techniques  or  otherwise)  is  only  an 
approximate  description  of  the  actual  plant.  For  instance,  a linear  model 
may  be  used  to  approximate  what  is  actually  a nonlinear  plant;  model  dynamic 
order  may  be  lower  than  that  of  the  actual  plant;  or  the  values  assigned  to 
model  parameters  may  be  only  approximations  to  those  in  the  plant,  and  the 
actual  parameter  values  may  vary  with  time  in  a nondeterministic  way.  These 
plant  uncertainties  may  result  in  unsatisfactory  performance  of  the  control 
system  when  implemented  with  the  actual  plant.  Thus  a need  exists  to  be  able 
to  adapt  the  controller  on-line  and  automatically  when  the  control  system 
is  sensed  to  be  deviating  from  its  intended  performance.  In  the  MRAS,  the 
reference  model  is  introduced  to  specify  the  desired  performance  and  an  adap- 
tation subsystem  is  provided  to  direct  the  necessary  compensating  adjustments. 

The  optimal  control  and  MRAS  approaches  to  control  system  design 
may  be  applied  as  complementary  design  techniques.  The  MRAS  may  be  viewed 
as  compensating  for  the  uncertainties  ignored  in  the  optimization  process. 

This  viewpoint  can  be  simply  illustrated  by  considering  the  simple  linear- 
quadratic  regulator  optimal  control  problem  [3].  Consider  the  linear  time- 
invariant  plant  model  and  quadratic  performance  index  defined  by: 
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x ■ Ax  + B u 
o o 


(1.1-1) 


j = i P (x^Qx  + u^Ru)dt 
2 Jo 


(1.1-2) 


where  Q s 0,  R > 0,  and  (Ao>  B ) provides  nominal  plant  parameter  values. 
The  optimal  control  is  given  by: 


u*(t)  = -R_1BT  P x(t) 
o o 


(1.1-3) 


where  Pq  is  the  unique  positive-definite  symmetric  matrix  satisfying: 


T -IT 

-A  P - PA  - Q + PB  R B P = o 

O O 0 0 


(1.1-4) 


Thus  the  theoretically  optimal  closed-loop  system  becomes 


x = (A  - B R_1BT  P ) x 

O O 0 0 


(1.1-5) 


In  order  to  construct  a MRACS  which  complements  this  optimal  design,  the 
M-subsystem  would  be  chosen  as: 


x = (A  - B R_1BT  P ) x 
M o o o o M 


(1.1-6) 


0 

I0 

0 


while  the  actual  plant  feedback  control  would  be  implemented  with  a time- 
varying  adaptable  matrix  P(t),  where  P(o)  = P . The  adjustment  P(t)  from 
the  A-subsystem  would  compensate  for  the  likely  possibility  that  the  actual 
plant  matrices  (A,B)  were  not  exactly  equal  to  the  plant  model  matrices 
(Aq,Bo).  Thus  the  S-subsystem  (actual  controlled  plant)  appears  as 


xs  = (A  - BR_1BTP(t)xs 


(1.1-7) 


with  P(t)  determined  by  the  adaptation  design. 


1 . 2 Historical  Background 

The  design  philosophy  for  MRAS  has  evolved  through  two  distinct 
phases  since  about  1960.  The  earliest  published  work,  which  introduced  the 
first  phase,  is  attributed  to  Whitaker  [62],  describing  an  adaptive  flight 
control  system.  The  objective  in  this  phase  of  MRAS  development  is  to  design 
an  adaptation  algorithm  which  minimizes  a performance  index  measuring  either 
the  state-  or  output-error  between  the  M-  and  S-subsystems.  A gradient- 
related  optimization  method  [44]  is  usually  applied  to  generate  parameter 
adjustment  signals.  Other  representative  work  during  this  phase  is  found 
in  [25,49,52].  One  drawback  for  the  performance  index  approach  is  the  compro- 
mise necessary  between  state-error  convergence  rate  and  MRAS  stability.  The 
gradient  algorithm  gain  (or  step  size  for  discrete-time  problems)  is  required 
to  be  smaller  than  a system-dependent  bound  to  guarantee  MRAS  stability. 
However,  this  bound  also  places  an  upper  limit  on  the  state-error  convergence 
rate. 

An  effort  to  circumvent  the  potential  instability  inherent  in  this 
design  approach,  as  well  as  the  necessity  of  stability  analysis,  resulted 
in  the  second  phase,  where  the  design  philosophy  was  to  guarantee  stability 
of  the  MRAS.  Among  the  first  published  works  during  this  phase  were  [54,48]. 
The  stable  design  is  obtained  by  applying  Lyapunov's  second  method  [37] to 
the  MRAS  structure,  resulting  in  MRAS  state-error  which  is  asymptotically 
stable.  Other  representative  work  includes  [63,53].  A departure  from  the 
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design  approach  using  Lyapunov's  method,  but  still  resulting  in  a guaranteed 
stable  MRAS,  was  introduced  in  1969  by  Landau  [31].  In  this  approach,  the 
design  is  obtained  by  applying  the  Popov  hyperstability  criterion  [51]  to  a 
reconfigured  form  of  the  MRAS.  Representative  work  in  this  area,  done 
exclusively  by  Landau  and  his  coworkers,  includes  [33,35,36], 

Several  other  recent  developments  in  MRAS  theory  are  worth  noting 
here.  The  stable  MRAS  design  discussed  above  assumes  availability  of  both 
the  M-  and  S-subsystem  state  vectors.  Various  approaches  have  been  proposed 
which  deal  with  the  case  when  only  the  subsystem  outputs  are  available.  The 
state  variable  filter  approach  [40],  and  the  augmented  signal  synthesis 
approach  [46]  are  two  such  methods.  In  the  context  of  MRAS  identification, 
the  adaptive  observer  [6,30,47]  has  been  developed  as  an  extension  of  the 
Luenberger  observer  [41].  This  form  of  MRAS  is  designed  to  provide  plant 
identification  in  conjunction  with  state  estimation.  However,  question  has 
been  raised  as  to  whether  the  distinction  between  a standard  MRAS  identifi- 
cation system  and  the  adaptive  observer  really  goes  beyond  the  issue  of 
system  representation  [ 1] . 

Several  excellent  surveys,  providing  extensive  bibliographies  in 
this  field  are  [12,39,34] . 

1 . 3 Summary  of  Results 

In  this  thesis,  we  investigate  a series  of  topics  related  to  the 
parallel  MRAS  structure.  Except  for  Section  2.6,  the  analysis  throughout 
is  based  on  the  discrete-time  single-input  single-output  hyperstable  MRAS 
design  introduced  in  [35].  Chapter  2 considers  problems  for  the  standard 
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form  of  MRAS  introduced  in  Section  1.1.  In  Chapter  3,  a new  MRAS  structure 
is  introduced  which  has  multiple  reference  models. 

Section  2.1  summarizes  the  MRAS  design  presented  in  [35],  which 
will  be  referred  to  here  as  the  Landau  MRAS.  In  Section  2.2  we  propose  a 
non-minimal  state  variable  representation,  named  the  Input-Output  Delay 
representation,  for  the  M-  and  S-subsystems  of  the  Landau  MRAS.  This 
representation  provides  a systematic  basis  for  implementing  a MRAS  simulation. 
Although  implicit  in  Landau's  work,  it  is  not  developed  explicitly.  The 
M-  and  S-subsystems  of  the  Landau  MRAS  also  provide  for  an  assumption  of 
knowledge  about  the  number  of  plant  poles  or  zeroes  located  at  the  origin 
of  the  z-plane.  By  extending  the  dimension  of  the  parameter  vector,  we  are 
able  to  generalize  the  MRAS  to  the  case  where  only  the  dynamic  order  of  the 
linear  plant  is  assumed  known,  with  no  specific  knowledge  about  plant  pole 
or  zero  locations. 

In  Section  2.3,  we  obtain  necessary  and  sufficient  conditions  for 
controllability  of  the  1-0  Delay  state  representation.  We  then  show  how 
these  conditions  may  be  used  to  produce  a reduced-order  M-subsystem.  The 
Landau  MRAS  assumes  both  M-  and  S-subsystems  have  the  same  order.  This 
result  provides  the  capability  for  applying  the  Landau  design  to  MRAS  control 
problems  where  the  reference  model  is  chosen  to  be  of  lower  order  than  the 
plant.  In  such  a case,  we  have  demonstrated  by  simulation  that  the  S-subsystem 
adapts  itself  to  a system  whose  order  is  the  same  as  that  of  the  M-subBystem 
by  producing  pole-zero  cancellations. 
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Section  2.4  considers  the  problem  of  providing  asymptotic  stability 
for  the  parameter-error  in  MRAS  identification,  a property  not  inherent  in 
the  hyperstable  design.  First  we  develop  some  geometric  and  algebraic  repre- 
sentations for  the  equilibrium  subspace  of  the  parameter-error.  Based  on 
these  representations  we  find  necessary  and  sufficient  conditions  for  which 
this  subspace  consists  of  only  the  origin  of  the  parameter -error  space.  Since 
it  turns  out  that  these  conditions  depend  on  the  unknown  plant  parameters, 
however,  we  conclude  that  previously  published  sufficient  condition  results 
for  this  problem  are  the  practical  limit  in  actual  MRAS  design. 

In  Section  2.5,  we  examine  a practical  problem  relating  to  hyper- 
stability of  the  Landau  MRAS  when  applied  to  parameter  identification.  In 
this  case,  the  M- subsystem  (plant)  parameters  are  assumed  unknown,  and  so  a 
direct  solution  to  the  hyperstable  design  problem  does  not  appear  possible. 
Landau  has  suggested  an  approach  to  this  difficulty  which  does  not  resolve 
the  potential  for  MRAS  instability.  We  suggest  another  approach  which 
introduces  a time-varying  property  to  the  linear  subsystem  of  the  MRAS  in 
its  Popov  configuration.  We  suggest  some  approaches  which  might  be  taken 
to  investigate  the  stability  of  this  modified  version  of  the  Landau  MRAS. 

In  concluding,  we  show  that  by  substituting  a series-parallel  MRAS  structure 
for  the  parallel  structure,  the  hyperstable  design  no  longer  depends  on 
parameters  in  the  M-subsystem,  thus  eliminating  this  problem  altogether. 
However,  one  must  pay  a price,  as  Landau  has  shown,  when  operating  in  a 
noisy  environment. 
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In  Section  2.6,  we  consider  a class  of  MRAS  which  uses  discontinuous 
adaptation  algorithms.  Work  in  this  area  has  been  done  primarily  by 
researchers  in  the  Soviet  Union.  Using  results  from  the  theory  of  variable 
structure  systems  (VSS),  we  prove  stability  of  a particular  continuous  time 
MRAS.  This  viewpoint  is  different  from  that  taken  in  previously  published 
work.  In  previous  cases,  Lyapunov's  method  has  been  used  to  prove  stability. 

Here  we  show  an  equivalence  between  the  MRAS  problem  and  a VSS  measurable 
disturbance  rejection  problem,  and  then  apply  existing  stability  results  for 

I 

the  VSS  problem. 

Chapter  3 introduces  a new  MRAS  structure,  the  multi-model  MRAS. 

'j 

This  particular  form  of  MRAS  has  not  been  treated  in  the  literature  up  to 
this  time.  In  Section  3.1  the  structure  is  introduced  and  motivation  for 
interest  in  such  a structure  is  presented.  We  also  jointly  develop  a 
regulator-type  two-model  MRAS  (2M-MRAS)  and  a dynamic  game,  the  goal  being 
to  provide  a basis  for  examining  a possible  correspondence  between  the  two 
structures.  As  part  of  this  joint  development,  we  introduce  alternate  view- 
points which  the  players  might  adopt  in  obtaining  their  reference  models 
(M^- subsystems).  In  Section  3.2  we  examine  the  stability  of  the  2M-MRAS 
where  each  adaptation  subsystem  is  of  the  hyperstable  Landau-type  considered 
in  Chapter  2.  By  transforming  the  2M-MRAS  to  the  Popov  configuration  we 
identify  how  its  structure  differs  from  the  standard  MRAS  in  that  same 
configuration. 

In  Section  3.3  we  jointly  develop  a dynamic  game  and  a 2M-MRAS  of 
the  quadratic  regulator  type.  We  establish  a procedure  for  comparing  the 
equilibrium  solution  of  the  dynamic  game  and  the  parameter  equilibrium  of 
the  2M-MRAS  by  comparing  the  actual  closed-loop  Mash  equilibrium  solution 
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of  the  dynamic  game  with  the  equilibrium  structure  obtained  in  the  2M-MRAS. 

For  the  particular  case  under  consideration,  we  conclude  that  the  two 
equilibrium  solutions  represent  different  strategies  and  objectives.  The 
2M-MRAS  equilibrium  structure  is  then  obtained  explicitly  in  terms  of  the 
parameters  of  the  M^-subsystems . 

Section  3.4  considers  the  quasi -symmetric  character  of  the  2M-MRAS, 
and  we  obtain  results  for  the  regulator-type  2M-MRAS  which  define  how  devia- 
tions from  perfect  symmetry  effect  the  S-subsystem  parameter  equilibrium. 

Chapter  4 contains  conclusions  and  some  suggestions  for  further 
research  on  the  topics  addressed  in  this  thesis. 
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CHAPTER  2 


THE  SINGLE  MODEL  MRAS 

In  this  chapter  we  consider  the  class  of  parallel  MRAS 
introduced  in  Chapter  1.  This  class  is  defined,  by  a discrete- 

time  SISO-LTV  n— -order  subsystem  S,  whose  (2N+1)  parameters  vary  according 
th 

time  SISO-LTV  n — -order  subsystem  S,  whose  (2N+1)  parameters  vary  according 
to  the  output  of  an  adaptation  subsystem  A.  The  adaptation  subsystem  measures 
and  processes  M-  and  S-subsystem  variables  to  determine  S-subsystem  parameter 
changes  at  equally-spaced  time  intervals.  The  goal  of  subsystem  A is  to 
drive  to  zero  state-error,  and  possibly  parameter -error  between  M and  S. 

The  terminology  "single  model"  is  introduced  in  this  chapter  to 
contrast  this  class  of  systems  with  a class  to  be  studied  in  Chapter  3. 

Here,  a single  reference  model  is  assumed  in  the  MRAS,  which  specifies  the 
single  performance  objective  or  goal  for  the  overall  MRAS.  In  Chapter  3, 
we  consider  the  case  where  two  reference  models  are  assumed,  each  of  which 
specifies  a different,  possibly  conflicting,  performance  objective  or  goal 
for  the  overall  MRAS. 

The  adaptation  subsystem  A considered  in  Sections  2. 1-2.5  is 
continuous  with  respect  to  all  input  variables,  in  contrast  to  the  problem 
considered  in  Section  2.6.  There  we  will  study  the  class  of  systems  where 
subsystem  A is  defined  according  to  a difference  equation  with  discontinuous 
right-hand  side. 
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express  the  Landau  MRAS  algorithm  in  a compact  form,  and  provide  a systematic 
basis  for  implementing  the  algorithm  on  a digital  computer.  Next,  the  Landau 
MRAS  algorithm  is  generalized  to  a less  restricted  form,  called  the  Extended 
Landau  MRAS  Algorithm  (ELMA) . This  generalization  eliminates  the  assumption 
by  Landau  of  the  number  of  poles  or  zeroes  occurring  at  the  origin  of  the 
z-plane  in  the  M-subsystem.  Necessary  and  sufficient  conditions  are  then 
found  for  controllability  and  observability  of  the  1-0  delay  state  repre- 
sentation, in  terms  of  the  system  transfer  function.  Lastly,  we  obtain  an 
algebraic  relationship  between  the  1-0  delay  state  and  an  arbitrary  minimal 
state  representing  the  same  transfer  function.  The  motivation  here  is  to 
determine  the  feasibility  of  expressing  the  Landau  MRAS  algorithm  in  terms 
of  a minimal  dimension  (N)  state  vector,  rather  than  the  2N~dimensional 
1-0  Delay  State  vector. 

In  Section  2.4  we  consider  the  relationship  between  the  input 
sequence  u(k)  and  the  equilibrium  for  the  parameter -error  of  the  MRAS. 
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1 

This  topic  is  of  Importance  when  the  MRAS  Is  applied  to  the  Identification 
problem.  Results  relating  to  necessary  and  sufficient  conditions  for  the 
parameter-error  to  converge  to  the  origin  are  obtained. 

In  Section  2.5  we  consider  some  practical  difficulties  relating  to 
the  stability  of  the  MRAS  when  applied  to  the  identification  problem.  The 
positive-real  design  requirement  for  the  MRAS  cannot  be  directly  realized  due 
to  the  nature  of  the  identification  problem.  We  suggest  a design  modification 
to  the  Landau  MRAS  and  describe  methods  to  analyze  stability  properties  of 
this  modified  version. 

In  Section  2.6  the  MRAS  adaptation  algorithm  is  assumed  to  be  a 
discontinuous  function  of  the  measured  system  variables.  Both  parameter 
adaptation  and  signal  synthesis  formulations  are  presented.  Using  results 
on  the  stability  of  variable  structure  systems  designed  for  disturbance 
rejection,  we  establish  the  stability  of  the  MRAS  with  discontinuous 
adaptation. 

Appendix  A contains  definitions  and  notations  which  are  used  in 
the  development  of  this  chapter. 

2.1  The  Landau  MRAS  Design 

In  this  section  we  present  a mathematical  statement  of  the  Landau 
MRAS  design  |i35j,  in  order  to  develop  a basis  in  the  reader  for  understanding 
the  fundamental  problem.  From  this  basis  we  will  be  ready  to  examine  more 
specific  issues  related  to  the  MRAS  problem. 

The  basis  for  his  design  method  is  the  result  by  Popov  on 
necessary  and  sufficient  conditions  for  hyperstability  [50].  By  transforming 
the  parallel  MRAS  structure  to  the  linear -nonlinear  feedback  form  considered 
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by  Popov,  he  was  able  to  find  a particular  adaptation  subsystem  A so  that 
the  overall  MRAS  satisfies  the  Popov  asymptotic  hyperstability  conditions; 
i.e.,  the  MRAS  state-error  asymptotically  approaches  zero.  Figure  2.1-1 
depicts  the  particular  parallel  MRAS  considered  by  Landau,  and  Figure  2.1-2 
depicts  that  system  transformed  to  the  linear-nonlinear  structural  form  used 
for  hyperstability  analysis.  Definition  of  the  equations  in  each  block 
appears  later  in  this  section. 

Define  the  following  (n+rofl)  vectors: 


Figure  2.1-1  Landau  MRAS--  Discrete-Time  Linear  SISO  Subsystems  M and  S 


Linear  Subsystem  (LS) 


FP-S256 


Figure  2.1-2  Transformation  of  Landau  Parallel  MRAS  Structure  for 
Popov  Hyperstability  Analysis 


Then  the  three  MRAS  subsystems  M,  S,  and  A are  defined  to  be  (See  Figure 

2.1-1): 


yM(k) 

" Pm  xM(k) 

(2.1-7) 

ys(k) 

= pT(k+l)  xg(k) 

(2.1-8) 

ys(k) 

= pT(k)  xg(k) 

(2.1-9) 

e(k)  = 

yM(fc)  - ys(k) 

(2.1-10) 

c°(k) 

- yM(k)  * yjjoo 

(2.1-11) 

v (k)  = 

«(k)  + cTe(k) 

(2.1-12) 

v°(k) 

= e°(k)  + cTe(k) 

(2.1-13) 

P(k+1) 

= P(k)  + [F(k)  xs(k)]v(k) 

(2.1-14) 

* P(k)  + [F(k)  x (k)]  - 

1 + Xg(k)  F (k)  xs(k) 

[F(k)  xs(k)][F(k)  xg(k)]T 
1 + x*(k)  F (k)  x (k) 


F(k+1)  - F(k)  - 


(2.1-15) 
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1 

F(0)  > 0;  symmetric 

It  may  be  shown  that  (2.1-12)  and  (2.1-13)  are  related  by: 

v(k)  = 

1 + Xg(k)  F(k)  Xg(k) 

The  variables  which  appear  both  without  and  with  superscript  o 
(e-g*>  yg(k),  y°(k))are  referred  to  by  Landau  as  a posteriori  and  a priori 
variables,  respectively,  since  the  superscripted  variables  are  generated 
before  updating  jT(k)  while  the  non-superscripted  variables  are  generated 
after  p(k)  is  updated  to  p(k+l). 

Note  also  that  no  recursive  expressions  are  written  for  the  'Jc(k) 

vectors,  although  these  expressions  should  be  provided  for  completeness. 

This  aspect  of  the  algorithm  will  be  dealt  with  in  the  next  section. 

The  Popov  structure  in  Figure  2.1-2  is  obtained  from  (2.1-7  - 

2.1-15),  where  6p(k)  = p - p’(k)  is  the  parameter-error  The  state-reali- 

M 

zation  of  the  linear  subsystem  (LS)  is  defined  by: 


= [0  0 1]; 


(2.1-17) 


(2.1-18) 


(2.1-19) 
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i.e.,  this  is  a phase-canonic  (completely  controllable)  realization. 

In  developing  the  necessary  and  sufficient  conditions  for  hyper- 
stability, Popov  assumes  the  nonlinear  subsystem  satisfies  an  inequality 
constraint: 

K 2 

I z(k)  v(k)  £ -y  (2.1-20) 

k=o  ° 

2 

for  all  K,  vo  a positive  real  constant. 

However,  (NS)  is  actually  a linear  time-varying  block  with  respect  to  the 
variables  v(k)  and  z(k),  with  state  vector  an<*  time-varying  coeffi- 

cients computed  from  xg(k)  and  F(k).  Thus,  a sufficient  condition  for  (NS) 
to  satisfy  (2.1-20)  is  that  (NS)  satisfies  the  conditions  of  the  positive 
real  lemma  extended  to  time-varying  discrete  systems  [36]. 

The  transfer  function  of  (LS)  is: 


LS(z)  =1-1 


i?i(aM+Ci)zn"1 


n r 

z 'til  \ 2 


(2.1-21) 


It  is  necessary  and  sufficient  that  c€Rn  be  chosen  so  that  LS(z)  is  a 
discrete  positive  real  transfer  function  [16].  This  discrete  positive  real 
concept  is  defined  in  Section  2.5  This  is  not  a trivial  problem,  and 
unless  proper  attention  is  paid  to  this  design  phase,  the  algorithm  will 
necessarily  lose  its  hyperstability  property.  Setting  c * -aj^  clearly 
satisfies  the  positive  real  condition.  However,  for  the  identification 
problem,  aM  corresponds  to  part  of  the  unknown  plant  parameter  vector,  sb 
setting  c * -aM  is  not  possible.  Landau  suggests  obtaining  an  initial  estimate 
a for  the  S-subsystem  parameter  a(k)  in  (2.1-4),  and  then  choose  c “ -a. 
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Although  this  approach  does  not  guarantee  satisfaction  of  the  positive  real 
condition,  Landau  claims  it  has  not  failed  him  in  achieving  state-error 
convergence  to  zero.  In  Section  2.3  we  add  another  condition,  not  addressed 
by  Landau,  which  c must  satisfy,  based  on  assumptions  made  by  Popov.  In 
Section  2.5  we  consider  in  more  detail  the  effects  of  choosing  the  vector  c. 

2.2  Input-Output  Delay  - A State  Representation  for  the  Landau  MRAS 

The  difference  equation: 

m n m 

y(k)  = p x(k)  = i21  a£  y(k-i)  + tZQ  bi  u(k-i)  (2.2-1) 

representing  either  subsystem  H or  S of  the  Landau  MRAS  in  Section  2.1  may 

s 

be  expressed  both  as  a transfer  function  and  in  state  variable  form.  In 
this  section  we  find  the  transfer  function  for  arbitrary  (n,m),  and  express 
che  difference  equation  in  a state  variable  form.  This  form  provides  a 
basis  from  which  we  will  later  analyze  the  MRAS  subsystems,  it  completes 
the  dynamic  expression  of  the  MRAS  algorithm,  and  provides  a systematic  basis 
for  implementing  the  algorithm  on  digital  computer. 

The  structure  of  (2.2-1)  suggests  the  name  input-output  delay 
difference  equation,  and  henceforth  we  will  refer  to  it  by  that  name. 

Taking  the  Z-transform  of  (2.2-1)  we  obtain: 

n -i  m -i 

(1  - Z a z L)  y(z)  = (Z  b.z  l)  U(z)  (2.2-2) 

i*l  i«*o 
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Then  the  transfer  function  Is: 


£ b z 

i 1 


n 

-i 

1 - £ a . 

z 

i-1  1 

z(n*ra>  I b< 

m-i 

z 

1-0  1 

n J? 

n-i 

z - £ a 
i-1  1 

z 

n 

r k n_i 
£ b z 

i-o  1 

n £ n-i 

z - £ a z 

i-1 


if  n > m 


if  n ■ m 


£ b z 
l-o 


(m-n)  r n J?  _ n-i, 
'■z  ’i-1  1 2 J 


if  n < m 


We  note  that,  when  n t m,  either: 

1)  at  least  (n-m)  zeroes  exist  at  the  origin  out  of  a maximum  of 
(n-1)  possible  zeroes,  or 

2)  at  least  (m-n)  poles  exist  at  the  origin  out  of  a total  of 
m poles 


*1 


I 

I 

I 

I 

r 


0 


i: 


c 


e 

0 


D 


0 
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Denote  N ■ max(n,m)  as  the  order  of  the  transfer  function.  We  may  nov  write 
a single  expression  for  H(z): 


m 

I b. 


N-i 


H(z) 


i*o 


zN-e" 


i-1 


ai  Z~ 


N-i 


(2.2-4) 


Either  n or  m may  be  less  than  N.  However,  there  is  no  loss  in  generality 
.T  p t . T 


by  extending  a ■»  [a  a,]  or  b ■ [b  •••  b.]  to  dimension  N,  and  imposing 

n i m l 

either  [aN  • ••  a^^]  * 0 or  [b^,  ••*,  b j]  - 0,  whichever  is  appropriate. 
Then  the  system: 


N 


y(k)  * 2 (a  y(k-i)  + b u(k-i))  + b u(k) 
i-1  1 1 ° 


(2.2-5) 


N 

2 


H(z) 


i*o 


biZ 


N-i 


N 


N TT 
zW  - 2 


£<bi  + aibo>z 


N-i 


i-1 


V" 


N-i 


N ^ iTT 
z - 2 a z 

i-1  1 


+ b 


(2.2-6) 


is  equivalent  to  (2.2-1,  2.2-4),  with  appropriate  coefficients  set  to  zero. 
We  will  refer  to  (2.2-5,  2.2-6)  as  the  extended  version  of  the  system 
(2.2-1,  2.2-4). 

Consider  now  the  2N-dimensional  state  variable  definition: 


I ,T 


xT(k)  - [x°  (k)  I x1  (k)]  = Cy(k-N)  •••  y(k-l)  | u(k-N)  •••  u(k-l)] 


(2.2-7) 


■ 

i 


|3 

I 


t* 


j 


23 


Then  a state  variable  realisation  of  the  extended  system  (2.2-5,  2.2-6)  Is: 


where 


x(k+l)  - F x(k)  + G u(k) 
y(k)  - pT  x(k)  + bQ  u(k) 


r ><i>  | 


gt  - [gT  !g*]  - O-Ob  jo... 01] 
1 | 2 O I 


- [« 


jbT]  -c.N. 


N 


•bi] 


(2.2-8) 

(2.2-S) 

(2.2-10) 

(2.2-11) 

(2.2-12) 


The  matrices  F,  F(l),  and  I*  are  fully  defined  in  Appendix  A. 

The  validity  of  (2.2-8.-  2.2-12)  as  a 2 N- dimensional  realization 
of  (2.2-5,  2.2-6)  may  be  established  by  inspection  of  the  state  definition, 
and  recognizing  that  x°(k+l)  » y(k).  Clearly  y(k)  is  identically  defined  in 
(2.2-5)  and  (2.2-9).  We  will  call  this  state  representation  the  1-0  Delay 
State  Representation. 

The  motivation  for  extending  both  a and  b to  R1*  from  Rn  and  Rm, 
respectively,  lies  in  the  intention  to  remove  from  the  Landau  algorithm  a 
potentially  restrictive  assumption  which  he  implicitly  makes.  The  system 
(3.3-1,  2.2-3)  assumes,  not  only  the  denominator  and  numerator  dimensions 
(i.e.,  basic  system  structure),  but  further  that  both  subsystems  M and  S 
have  either  (n-m)  zeroes  at  the  origin  (when  n > m) , or  (m-n)  poles  at  the 
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origin  (when  m > n).  By  extending  the  pereaeter  vector*  e end  b to 

R , we  choose  to  assume  basic  structure  only,  allowing  the  identification 

process  to  determine  whether  or  not  there  are  poles  or  zeroes  at  the  origin. 

This  is  accomplished  by  observing  whether  either  [aM(k)...  a (k)]  or 

N 0+1 

cVk>  ■■■  Vi  (k)]  converges  to  zero  as  a result  of  the  adaptation  process. 
We  will  refer  to  the  algorithm  using  the  extended  parameter  set  as  the 
Extended  Landau  MRAS  Algorithm  (ELMA). 

The  1-0  Delay  State  Representation  and  the  associated  realization 
for  (2.2-6)  provide  a systematic  method  for  recursively  generating  x(k). 

We  make  use  of  this  representation  in  stating  the  Landau  MRAS  algorithm  in 
its  extended  form  (ELMA): 


M: 

xM(k+i)  - Fm  xM(k)  + G„  u(k) 

(2.2-13) 

yM(k>  * xM(k>  + bM  °(k) 

0 

(2.2-14) 

S: 

xg(k+l)  - F(k+1)  xg(k)  + G(k+1)  u(k) 

(2.2-15) 

yS(k)  “ PT<k+1>  xS<k>  + b0(k+1)  u(k) 

(2.2-16) 

y°(k)  - pT(k)  xs(k)  + bo(k)  u(k) 

(2.2-17) 

This  subsystem  remains  unchanged,  except  that  we  define 

state-error: 

e(k+l)  - ij  e(k)  + g2  e (k) , 

(2.2-18) 

I 

I 
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xs(k)  la  formed  by  appanding  u(k),  the  currant  input,  and  7(k) 
corraaponda  now  to  tha  axtandad  parameter  sat.  An  altarnate  dynamic  aquation 
for  o(k),  obtainad  diractly  from  (2.2-13,  2.2-15)  ia: 

•(k+l)  - Fm  e(k)  + G2  Xg(k)  6P(k+l)  (2.2-19) 

Va  conclude  this  section  by  presenting  simulation  results  of  two 
examples  of  tha  Extended  Landau  MRAS  Algorithm  (ELMA).  Tha  first  example 
is  of  a 1—  order  problem  (3  parameters),  the  second  of  a 2— ~ order  problem 
(5  parameters).  In  each  example  we  consider  both  pulse  input  and  sinusoidal 
input.  Time  discretization  is  taken  as  0.1  second.  Initial  state  conditions 
for  the  M-  and  S-subsystems  are  assumed  zero. 

Example  2.2-1  - 1—  Order  Problem 
M-subsystem: 

a)  pole  at  z * 0.9 

b)  b - 0.1 

o 

c)  b^  chosen  to  yield  unity  gain  on  |z|  - 1 

Thus  pT  - [0.9  I 0,  0.1] 

m l 

Initial  S-subsystem: 

a)  pole  at  z - 0 

b)  S-transfer  function  numerator  • M-transfer  function 
numerator 

Thus  pT(0)  - [0  j 0.09,  0.1] 
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A-aubaystem: 

■)  F(0)  - diag  [1000,  1000,  1000] 

b)  c (t)  - -a(t) 


Input  u(k): 

case  (a)  - pulse  with  (amplitude,  period)  • (1,5) 
case  (b)  - sine  with  (amplitude,  period)  * (1,4) 


Final  S-subsystem: 

P(10) 

case  (a) 

[0.9  \ 6 x 10‘5,  0.1] 
1 

case  (b) 

[0.904J  -0.012,  0.11] 

Simulation  results  are  shown  in  Figure  2.2-1.  In  both  input  cases, 
state-error  e(k)  converges  to  zero,  in  accordance  with  the  hyperstable  pro- 
perty of  the  MRAS.  Furthermore,  parameter-error,  6p(k)  - p^  - 7(k),  converges 
to  zero. 


n 

c 


Example  2.2-2  2 — Order  Problem 

M-subsystem: 

a)  poles  at  z ■ (0.9,  0.4) 

b)  zero  at  z - 0.8 

c)  b * 0 

o 

d)  b^  chosen  to  yield  unity  gain  on  | z)  * 1 
Thus  p£  - [-0.36,  1.3  ! -.24,  .3,  0] 

M t 


F 


a 


h* 


(3)  S-Pole 


(4)  S-Numerator  Coefficient 

FP-5222 


Figure  2.2-l(a)  ETLMA  Example  2.2-1— Case  a 


(3)  S-Pole 


(4)  S-Numerator  Coefficients 

FP-5223 


Figure  2 .2- 1(b)  E1MA  Example  2.2-1— Case  b 
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Initial  S-subsystem: 

a)  poles  at  z - (0,0) 

b)  S-transfer  function  numerator  - M-transfer  function 
numerator 

Thus  pT(0)  - [0,0  | -0.24,  0.3,  0] 


A-subsystem: 

a)  F(0)  - diag  [1000,  1000,  1000,  1000,  1000] 

b)  c (t)  * - a(t) 

Input  u(k): 

case  (a)  - pulse  with  (amplitude,  period)  * (1,5) 

case  (b)  - sine  with  (amplitude,  period)  * (1,4) 

case  (c)  - sine  with  (amplitude,  period)  * (1,20)  + (10,3) 


Final  S-subsystem: 
case 


PT(12) 


(poles,  zeroes) 


«) 

b) 

c) 


[-0.356,  1.3  j -0.239,  0.3,  10"4] 
[0.105,  0. 754  j -0.263,  0.392,  0.008] 
[-0.361,  1.3  j -0.242,  0.301,  1.6xl0-4] 


(.898,  .391  J .792) 
(-0.12,  0.874  j 0.66) 
(0.901,  0.4  | 0.801) 


Simulation  results  are  shown  in  Figure  2.2-2.  In  all  three  input 
cases,  the  state-error  converges  to  zero.  Notice,  however,  that  parameter 
convergence  fails  for  case  (b) , where  the  sine  input  contains  only  a single 
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Figure  2. 2-2 (a)  EUMA  Example  2.2-2— Case  a 
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frequency.  The  Landau  MRAS  design  guarantees  only  that  the  state-error  is 
asymptotically  hyperstable,  and  claims  nothing  more  than  stability  (bounded- 
ness) for  the  parameter  error  vector  6p(k)  = p„  - p(k).  In  Section  2.4  we 
investigate  this  aspect  of  the  MRAS  problem  in  detail. 


2.3  Structural  Properties  of  the  1-0  Delay  State  Representation  - 
Implications  for  the  MRAS 

In  this  section  we  derive  a series  of  properties  relating  to  the 
structure  of  the  1-0  Delay  state  representation.  First  we  examine  the  con- 
ditions for  controllability  and  observability  of  the  1-0  Delay  state  repre- 
sentation. This  leads  to  consideration  of  the  effects  of  pole-zero  cancella- 
tion (and  loss  of  1-0  Delay  state  controllability)  in  the  M-subsystem  on  the 
performance  of  the  Extended  Landau  MRAS  Algorithm  (ELMA). 

Next  we  find  an  algebraic  relation  between  the  2N-dimensional  1-0 
Delay  state  and  an  arbitrary  N-dimensional  state,  each  of  which  realizes 
the  same  transfer  function  (2.2-6).  This  leads  to  the  problem  of  expressing 
the  1-0  Delay  state  as  a function  of  a minimal  state.  The  motivation  here 
is  to  determine  whether  expressing  the  ELMA  in  terms  of  minimal  state  repre- 
sentations of  subsystems  M and  S will  introduce  computational  economies. 

The  conditions  for  controllability  and  observability  of  the  1-0 
Delay  state  realization  (2.2-8  - 2.2-12)  may  be  determined  by  examining  the 
controllability  and  observability  matrices: 


R = [G  FG 
c 


i2"-1  G] 


(2.3-1) 
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I 

r 

•* 

D 


0 


0 


T ~ 
P-  F 


T ~2N-1 
P F 


(2.3-2) 


Using  (2.2-10  - 2.2-12),  and  Rq  may  be  expressed: 


h(0)  h(l)  ...  h(N-l)^  h(N)  F h(N)  •••  FN-1  h(N)j 

tf  j “ ~0  I 


N 


(2.3-3) 


where  h(k)  = F GL  + F(k)  G2  , k - o,  •••,  N , 


(2.3-4) 


F , F(k),  PN  are  defined  in  Appendix  A,  G^  and  G2  in  (2.2-11), 


and 


/_  l*®.  J 


2N  I N A 
F^”  , FN  F(N) 


(2.3-5) 


The  next  4 lemmas  will  be  required  in  what  follows.  The  first  is 
a useful  and  well-known  fact  from  matrix  algebra  [10]. 


Lemma  2.3.1 

Given  a square  matrix  partitioned  into  4 submatrices 


I 


A 1 A 
_1_  |_  2_ 


A I A 

L_  I 1 


1 * 

t,* 


rss 


» $ 

LT-> 
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assuming  either  or  A^  is  nonsingular,  then 


det(A)  * det(A1)'det(A^  - A^ 

= det(A^)  •det(A1  - A^"^) 


Lemma  2.3.2 


Two  arbitrary  SISO  LTI  systems,  Sj  = (C^A^B^  and  S2  = (C2,A2>B2) 
are  realizations  for  the  same  transfer  function  if  and  only  if 


k k 

C^B^  = C2A2B2  for  all  k » 0,1,2,*  •• 


Proof:  The  relaxed  output  response  to  an  input 


u(k)  = 1 if  k=0 

l 0 otherwise 


. {y(0),  y(i),  y(2) , •••  }, 


called  the  input  weighting  sequence,  for  each  system  is: 

{0,  C^B^,  C^A^B^ , C^AjB^,*‘*J  and  {0,  C^Bj,  ^2^2^2*  ^2^0^2,***^, 
respectively. 

The  relaxed  output  response  from  two  equivalent  system  realizations 
must  be  equal.  Therefore  their  respective  input  weighting  sequences  must  be 
equal.  By  definition  of  the  input  weighting  sequence,  two  equal  input  weight* 
ing  sequences  imply  two  equivalent  systems;  i.e.,  they  realize  the  same 
transfer  function.  g 

In  the  next  lemma,  we  consider  the  SISO  N-dimensional  system: 


w(k+l)  = F w(k)  + h(N)  u(k) 


(2.3-6) 
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y(k)  * c w(k)  + bQ  u(k) 


(2.3-7) 


where  h(N)  is  defined  by  (2.3-4),  F from  Appendix  A. 2,  and  c arbitrary. 


Lemma  2.3.3 

The  system  (2.3-6,  2.3-7)  is  a completely  observable  realization 
of  the  transfer  function  (2.2-6)  if  and  only  if  c^  * [1  0 • • • 0]. 

Proof;  The  1-0  Delay  realization  and  the  above  system  each  produce 
the  following  N-stage  input  weighting  sequences,  respectively: 


T — 
p FG 


T ~N-1~ 
p F G 


cT  h(N) 
cT  F h(N) 


cT  FN-1  h(N) 


= [FN!  F(N)]  G = FNgx  + F(N)  G2  = h(N) 


T N- 1 
c F 


c F b(N)  = K h (N) 


From  Lemma  2.3.2,  the  input  weighting  sequences  are  equal  since 

T 

both  systems  realize  the  same  transfer  function.  Thus  c must  satisfy  the 


equality: 


R - 1° 
o N 


Using  the  definition  for  F in  Appendix  A. 2,  we  see  that  the  equality  holds 
only  if  cT  - [1  0 •••  0]  . If  cT  = [1  0 • • • 0] , Rq  = 1°  , so  that  the  input 
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I 

l 

0 

c 

c 


weighting  sequences  for  the  two  system  realizations  are  equal.  Again  by 
Lemma  2.3.2,  these  two  realizations  must  realize  the  same  transfer  function, 
namely  (2.2-6).  Further,  since  Rq  * 1°  is  nonsingular,  and  is  the  observa- 
bility matrix  for  (2.3-6,  2.3-7),  that  system  is  completely  observable. 


Lemma  2.3.4 


The  system  (2.3-6,  2.3-7)  with  c = [1  0 0] , which  realizes 


.th 


the  N — order  transfer  function  (2.2-6),  is  completely  controllable  if  and 
only  if  N is  the  minimal  system  order;  i.e.,  there  are  no  pole-zero  cancella- 
tions in  (2.2-6). 

Proof:  Since  the  system  is  already  observable,  from  Lemma  2.3.3, 
controllability  implies  the  realization  is  minimal  [22],  which  implies  no 
further  possible  reduction  in  the  order  of  (2.2-6).  A minimal  realization  is 
always  completely  controllable  [22].  ^ 

Having  stated  these  lemmas  we  are  now  prepared  to  state  and  prove 
the  main  results  of  this  section. 

Proposition  2.3.1 

The  2N-dimensional  1-0  Delay  state  realization  (2 .2-8)-(2 .2-12) 
for  the  transfer  function  (2'.2-6)  Is  completely  controllable  if  and  only  if 
N is  the  minimal  system  order;  i.ev,  there  are  no  pole-zero  cancel lat lone  in 
*(2.2-6).  Thus,  controllability  depends  on  the  parameter  vector  ? only  to  the 

r- ■ y ■*  _ 2N+1 

extent  that  the  location  of  p e R doed  not  incroduce  common  roots  between 
the  numer&tor  and  ddhotain&tor  of  (2.2-6). 


Proof:  The  1-0  Delay  state  realization  is  completely  controllable  if  and 

only  if  (2.3-3)  is  a nonsingular  matrix  [3].  Since  the  rank  of  a matrix  is 
unaffected  by  reordering  its  columns,  we  rewrite  (2.3-3)  as: 


N-l 

FW  h(N) 


F h(N)  h(N)  | h(N- 


(N-l)  •••  h(l)  h(0) 


Using  Lemma  2.3.1,  det(Rc)  = det([F  h(N)  •••  F h(N)  h(N)]).  But  non- 
singularity of  is  now  reduced  to  nonsingularity  of  the  controllability 

matrix  for  the  system  of  Lemma  2.3.4.  Therefore  R is  nonsingular  if  and 

c 

only  if  N is  the  minimal  system  order. 


Proposition  2.3.2 

The  1-0  Delay  state  realization  for  the  transfer  function  (2.2-6) 
is  not  completely  observable,  independent  of  the  parameter  vector  f>. 

Proof:  The  1-0  Delay  state  realization  is  completely  observable  if  and  only 

if  (2.3-5)  is  a nonsingular  matrix  [3].  But  Rq  is  clearly  singular  since 
the  lower  half  of  the  matrix  is  a linear  combination  of  the  upper  half: 

[F2N  j fn  F(N)]  = FN[FN  j F(N)]  b 

Given  the  necessary  and  sufficient  conditions  for  controllability 


of  the  1-0  Delay  state  realization  in  Proposition  2.3.1,  consider  the 
following  situation  which  may  occur  in  the  M-subsystem  of  the  MRAS.  Let 
be  such  that  the  corresponding  transfer  function 
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(2.3-7) 


i-1  Mi 

has  at  least  one  pole  and  zero  in  common,  giving  subsystem  M (in  effect)  a 
structure  of  lesser  order  (by  at  least  1)  than  that  of  the  adaptable  subsystem 
S.  We  introduce  the  notation  G(M)<6(S)  to  indicate  this  condition.  For  MRAS 
identification,  this  situation  might  correspond  to  an  identification  model  with 
a higher  order  than  that  of  the  plant  to  be  identified;  for  MRAS  control,  this 
might  correspond  to  a reference  model  with  a lower  order  than  that  of  the 
plant/controller  subsystem  S.  The  MRAS  control  case  is  of  particular  interest. 
If  we  consider  a MIAS  structure  where  the  reference  model  is  characterized  by 


an  1-0  delay  difference  equation  of  the  form  (2.2-1)  with  2NM+1  parameters, 
while  the  plant  is  characterized  by  2N.+1  parameters,  N>NU,  then  it  is  not 
clear  that  the  Landau  MIAS  would  still  retain  its  hyperstable  property.  Thus 
the  pole-zero  cancellation  strategy  in  defining  a lower  order  reference  model 
permits  us  to  retain  the  hyperstable  property  of  the  MRAS. 

We  investigate  the  effect  of  M-subsystem  pole-zero  cancellation  on  the 
movement  of  the  "excess"  adaptable  poles  and  zeroes  in  subsystem  S,  through 
simulation  examples.  Refer  first  to  Example  2.2-2,  case  c.  That  2—  order 
problem,  without  pole-zero  cancellation,  was  successful  in  converging  both  state- 
error  and  parameter-error  to  zero.  Now  consider  the  same  problem,  except  where 
the  zero  is  shifted  to  cancel  first  the  pole  at  z*0.9,  then  the  pole  at  z«0.4. 


Example  2.3-1 

The  problem  is  set  up  exactly  as  in  Example  2.2-2,  case  (c),  except 
that  the  sine  input  is  (amplitude,  period)  ■ (1,20)  + (4,3) , and  the  M-subsystem 
zero  is  chosen  to  cancel  one  of  the  poles.  The  two  cases  considered  here  are: 


- 
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case(a):  M-subsystem  poles  it  z ■ (0.9,  0.4) 

zero  at  z - 0.9 

Thus  p£  * [-0.36,  1.3  S -0.54,  0.6,  0] 

M I 

case(b):  M-subsystem  poles  it  z 1 (0.9,  0.4) 

zero  at  z * (0.4) 

Thus  p£  = [-0.36,  1.3  | -0.04,  0.1,  0] 

The  simulation  results  are  as  follows,  also  depicted  in  Figure  2.3-1. 


Final  S-subsystem: 


In  both  cases  one  S-subsystem  pole  converges  to  the  uncancelled  M-subsystem 

pole;  the  other  S-subsystem  pole  and  the  zero  tend  to  a common  value,  although 

not  the  pole-zero  value  chosen  in  the  M-subsystem.  Thus  both  subsystems' 

s t 

outputs  appear  as  if  from  1—  order  systems.  The  final  common  pole-zero 
value  appears  to  be  influenced  by  the  choice  of  input.  Further,  we  can  show 
through  simulation  that  a single  input  frequency  is  sufficient  to  cause  the 
same  type  of  parameter  convergence  property  as  has  been  demonstrated  here. 
Notice  that,  in  Example  2.2-2,  case  (b)  with  5 parameters,  a single  input 
frequency  was  insufficient  for  parameter  convergence,  although  it  was 
sufficient  for  the  1—  order  problem  in  Example  2.2-1,  case  (b),  with  3 
parameters.  This  fact  reinforces  the  contention  that  an  M-subsystem  pole-zero 
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Figure  2.3-l(a)  EIMA  Example  2.3-l--Case  a 


cancellation  is  recognized  by  the  adaptation  process,  which  then  solves  a 

g £ 

1 — order  problem,  while  causing  an  arbitrary  S-subsystem  pole-zero 


cancellation. 


0 


In  concluding  the  discussion  on  this  point  we  consider  one  more 

aspect  of  pole-zero  cancellation;  in  this  case  how  it  might  affect  the 

hyperstability  properties  of  the  MRAS,  as  established  by  Landau.  Popov 

[50]  assumes  the  linear  subsystem  of  Figure  2.1-3  to  be  non-degenerate; 

that  is,  in  the  MRAS  context,  the  state  error  e(k)  is  assumed  completely 

T 

controllable  through  B„  and  completely  observable  through  Cu,  where 
T 

(Cm*  Fm,  Bm)  corresponds  to  the  phase-canonic  realization  for  (LS),  as 
defined  in  (2.1-17)  - (2.1-19).  We  note  again  that  (LS)  is  a hyperstable 
subsystem  if  and  only  if  c « R1*,  in  (2.1-17),  is  chosen  to  guarantee  (LS) 
has  a positive  real  transfer  function  (2.1-20).  However,  one  further 
assumption  which  must  be  made  in  choosing  c,  in  order  for  (LS)  to  satisfy 
Popov's  non-degeneracy  assumption,  is  that  the  parameter  vector  for  (LS): 


-T  r I ! T!  il 
PLS  ' [aM  !C  1] 


(2.3-8) 
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provide  no  pole-zero  cancellations  in  (2.1-21).  Otherwise  (LS)  will  be 
degenerate.  Landau  does  not  consider  this  aspect  of  selecting  the  vector 
c in  his  MRAS  design. 

We  turn  now  to  the  consideration  of  obtaining  a linear  algebraic 
relation  between  the  2N-dimensional  1-0  Delay  state  and  an  arbitrary 
N-dimensional  state,  each  of  which  realizes  the  same  N---  order  transfer 
function. 


Proposition  2.3.3 

Given  an  arbitrary  N-dimensional  realization  for  the  transfer 
function  (2.2-6): 

w(k+l)  = A w(k)  + B u(k)  (2.3-9) 

y(k)  = cT  w(k)  + bQ  u(k)  (2.3-10) 

and  the  1-0  Delay  state  realization  (2.2.8  - 2.2-12),  the  state  variables 
w(k)  and  x(k)  are  related  by 


R0  w(k) 


T A 

c A 


T . N-l 
c A 


w(k)  = [FN j F(N)]  x(k) 


(2.3-11) 


_N 


where  the  matrices  F , f’(N)  are  defined  in  Appendix  A. 2. 


Proof:  Since  both  realizations  represent  the  same  transfer  function,  the  zero- 

state  output  responses  to  an  arbitrary  input  must  be  the  same  for  both.  Construct 


» 
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the  N-stage  output  sequence  for  each  realization,  Y (N)  * {y(o),  y(l),*' 
y(N-l). 


Y(N) 


c 

CA 

CA 


N-l 


w(o)  + 


o 

CB 

CAB 


OAN-2„ 
CA  B 


o 

CB  b 


o o 
b 


i . , X 

a | b 

aTF  J bTl^  + aT  F (1) 

aV'1  ! b^1  + aT  F(N-l) 

I _ 


^(o) 

x2  (o) 


u(o) 

u(l) 

u (N-l) 


O t 0 0 o 


T ~ 
P G 


X 

p F G 


T ~N-2  ~ 
p F G 


T ~ 
P G 


o 

b 


u(o) 

u(l) 


u(N-l) 


According  to  Lemma  2.3.2,  the  terms  involving  input  must  be  equal.  Using 
the  definitions  of  FN  and  F(N)  in  Appendix  A. 2,  we  obtain  (2.3-11).  I 

Two  specific  cases  of  N-dimensional  SISO  system  realization  for 
(2.2-6)  are  the  completely  controllable  (phase)  canonic  and  completely 
observable  canonic  forms  [3],  In  these  two  cases,  (2.3-9,  2.3-10)  become: 


w (k+1)  + F w (k)  + G0  u(k) 
c c l 

y(k)  = (b1  + bQaT)  wc(k)  + bQ  u(k) 

wo(k+l)  = (P°  FT  P°)  wQ(k)  + P°(b  + boa)  u(k) 

y(k)  = G*  P°  ww(k)  + bQu(k) 
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(2.3-12) 

(2.3-13) 

(2.3-14) 

(2.3-15) 


respectively,  where  a,  b,  bQ,  F,  G2  are  defined  in  (2.2-10  - 2.2-12),  and 


P^  is  defined  in  Appendix  A.l. 


Corollary  2.3.1 

The  completely  controllable  and  completely  observable  states,  w (k) 

c 

and  WQ(k),  respectively  are  related  to  the  1-0  Delay  state  by: 


+ bQFN  ^ wc(k)  - [FN  lF(N)]  x(k) 


,N 


PN  Wo(k)  = £F  I F(N)3  x(k) 


(2.3-16) 


(2.3-17) 


Proof:  This  follows  from  Proposition  2.3-3  by  direct  substitution,  and 

recognizing  from  Appendix  A the  definition  of  FN  and  that  P°  P°  * 1°. 

N N N 


1 


We  next  observe  the  following: 


D 


Proposition  2.3.4 


Given  the  1-0  delay  state  realization  x(k)  and  an  arbitrary  N- 
dimensional  state  realization  w(k)  of  (2.2-6),  there  does  not  exist  R ■ [ 


such  that 


:l(k>l  K 

w(k). 

[2  (k)J  R2 


Proof:  It  is  sufficient  to  consider  only  x^Ck)  = R^wOO.  x^(k)  = 

[u(k-N)  ...  u(k-l)],  and  clearly,  w(k)  above  is  not  sufficient  to  determine 

the  past  N inputs.  The  inputs  necessary  to  reach  w(k)  depend  also  on  the 

starting  point  w(k-N).  ■ 

It  is  possible  to  express  x^(k)  as  a linear  function  of  w(k-N) 

T 

and  (k) • Given  the  N-dimensional  realization  (C  ,A,B),  we  can  show  that : 


Xj  (k) 


CA  w<k4N>  + ctab 


x2(k)  (2.3-18) 


cV1 


T N-2 
CA  B 


However,  no  advantage  in  reducing  the  amount  of  ELMA  computation  is  attained 
by  using  w(k)  instead  of  x^k). 
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2.4  Parameter-Error  Equilibria  in  MRAS  Identification 

Landau's  MRAS  algorithm,  as  well  as  others  designed  from  Liapunov 
stability  theory,  guarantee  only  that  the  parameter -error,  6p(k),  between  the 
M-  and  S-subsystems  is  bounded,  not  asymptotically  stable,  with  respect  to 
the  origin  of  r^N+I.  in  this  section  we  examine  in  more  detail  the  character 
of  the  trajectory  6p(k),  particularly  its  admissible  equilibrium  region. 

What  constraints  may  we  develop  which  characterize  this  region?  Examination 
of  these  questions  provides  further  insight  into  the  stability  properties  of 
the  Landau  MRAS  or  the  ELMA,  as  well  as  into  how  one  might  augment  the  ELMA 
design  (through  appropriate  choice  of  input  signal)  to  extend  the  property 
of  asymptotic  stability  to  the  parameter-error.  A fundamental  issue  in  the 
study  of  parameter  identification  is  that  of  identifiability . If  the  MRAS 
identification  model  contains  more  parameters  than  are  necessary  to  characterize 
the  system  being  identified,  then  it  is  not  possible  for  the  MRAS  parameter- 
error  to  be  guaranteed  asymptotically  stable.  This  property  was  illustrated 
in  the  last  section  in  Example  2.3-1  where  the  identification  model  contained 
five  parameters  while  the  plant  was  capable  of  being  completely  characterized 
by  only  3 parameters.  In  the  remainder  of  this  section  we  will  assume 
identifiability  conditions  are  satisfied. 

We  first  survey  published  results  dealing  with  asymptotic  stability 
of  the  MRAS  parameter-error.  The  problem  has  received  detailed  attention  from 
several  authors,  and  several  different  sufficient  condition  results  on  the 
input  signal  have  been  obtained  which  guarantee  the  asymptotic  stability. 

Then  we  examine  the  implications  on  6p(k)  of  state-error  convergence,  which 
leads  to  a development  of  several  geometric  and  algebraic  representations 


for  the  6p-equi librium  subspace.  Then  we  consider  the  problem  of  guaranteeing 
6p(k)  asymptotically  stable,  taking  a viewpoint  based  on  the  6p-equilibrium 
subspace  representation.  Finally,  we  use  simulation  examples  to  substantiate 
that  the  sufficient  conditions  proposed  by  other  authors  are  conservative. 

In  considering  MRAS  identification  using  stability-based  design 
methods,  most  authors  mention  the  problem  of  guaranteeing  parameter-error 
convergence  to  the  origin.  Lion  [40]  appears  to  be  the  first  author  to 
confront  this  aspect  of  MRAS  design  in  detail  in  the  literature.  He 
provides  sufficient  conditions  to  guarantee  asymptotically  stability  of 
parameter -error . For  the  continuous -time  MRAS  identification  problem 
posed  by  Lion,  convergence  of  state-error  implies  an  orthogonality  condition 
of  the  form: 


H(t)  - x£(t)  6p*  - 0 (2.4-1) 

— T r 

where  6p*  denotes  the  final  parameter-error,  and  x^  * [yM  ,y  , . . . ,yM  , 

M o **1  n-1 

uq,u^, . . . .u^]  corresponds  to  a phase-variable  state-representation  obtained 

from  "state-variable"  filtering  of  the  plant  output  and  input,  yw(t) 

M 

and  u(t),  respectively.  Since  x^t)  is  time-varying,  6p*  is  constrained 
to  the  subspace  E*  obtained  from  the  intersection  of  the  continuum 
of  hyperspaces  lH(t)  | te[t*,®)},  there  t*  denotes  any  time  after 


MMM 


■ 


I 


* 


I 


c 

[ 


which  e(c)  has  reached  the  origin;  i.e., 

E*  - 0 H(t) 
t«[t  ,(*>) 
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(2.4-2) 


Lion's  strategy,  then,  is  to  choose  the  MRAS  input  u(t)  in  such  a way  that 
the  only  possible  solution  for  (2.4-2)  is  E = f o} ; that  is,  6p  = o. 
Sufficient  conditions  given  by  Lion  to  guarantee  6p  = o are: 

1)  u(t)  is  a sinusoid  with  at  least  q*  = distinct 

frequencies  fw,,***,  w 
I q 

2)  The  phase  shift  cp(w)  of  y(t)  in  response  to  each  w^, 

i = 1 ,2 , • • • , q*,  must  satisfy  cp(wi)  + ^tt,  k=o,  ± 1,  ± 2,-*- 

Kudva  and  Narendra  [29,30]  start  with  a frequency  domain  viewpoint 
for  a discrete-time  multi-input  MRAS.  They  obtain  rank  conditions  for  a 
matrix  constructed  from  the  input  sequence  which,  if  satisfied,  are  sufficient 
to  guarantee  asymptotic  stability  of  the  parameter-error.  This  result  differs 
from  other  approaches  in  that  it  neither  requires  an  explicit  frequency 
content  condition  on  the  input  u(k),  nor  assumes  a particular  input  class 
such  as  sinusoids.  However,  how  the  matrix  to  which  a rank  test  is  applied 
is  obtained  is  not  clear  in  their  paper. 

The  most  recent  result,  and  an  advancement  over  Lion,  is  reported 
by  Kim  and  Lindorff  [24],  By  considering  an  arbitrary  state-variable  repre- 
sentation of  a single-input  linear  plant  model,  they  show  that  the  matrix 
error  [fiA(t)  6B(t)]  converges  to  o if: 
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1)  u(t)  Is  periodic  with  at  least  q*  distinct  frequencies 
{w^  • • • , w^*} , where 
min 

q*  = integers  {q|  qi  R/2J 

R = max  { R } 

i - 1 * ‘ " ,n 

R^  = number  of  non-zero  elements  in  the  i—  row 
^ of  [6A(t)  fiB(t)] 


2)  For  at  least  one  state  variable,  the  phase  shift  cdj(w)  satisfies: 

9j(«i)  ^ f11  for  a11  q*; 

k=o , i 1,  ± 2,  • • • 


3)  Given  the  set  of  non-zero  columns  of  6A(0),  the  corresponding 
set  of  plant  state  variable  {x^  } must  be  linearly  independent 

j 

in  the  steady-state. 

When  the  state  realization  has  minimal  dimension  = n,  q*  5 n— ^ -,  with  strict 
inequality  often  the  case.  This  compares  favorably  with  Lion's  condition 
q*  - ,rl..+  Jn  xhe  authors  also  extend  this  result  to  multi-input  systems . 

We  will  now  use  the  ELMA  of  Section  2.2  to  analyze  some  properties 
of  the  parameter-error  6j>(k) . We  first  consider  the  implications  of  state- 
error  convergence  and  v°(k)  convergence  on  the  adaptation  process,  concluding 
that  6p(k)  must  reach  a fixed  point  &p*  in  either  case.  Then  we  develop 
several  alternate  geometric  Interpretations  which  characterize  the  subspace 
in  which  <5p*  must  necessarily  be  contained. 

Consider  the  ELMA  under  the  condition  that  state-error  convergence 
has  occurred;  that  is,  e(k)  ■ 0,  e(k)  ■ 0.  From  (2.1-12,  2.1-14)  we  conclude 
that: 


! 


1 

I 
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v(k)  = v (k)  = 0 


6p(k+l)  = 6p(k)  = &p* 


(2.4-3) 


Thus  the  parameter-error  6p(k)  necessarily  becomes  a stationary  point 
<5p*  when  state-error  convergence  has  occurred. 

Making  use  of  (2.1-11,  2.1-13)  and  taking  c^  = - a^(k) , v°(k)  may 
also  be  expressed  in  the  form: 


v°(k)  =«  xj(k)  5j>(k) 


According  to  (2.4-3),  v°(k)  will  have  converged  to  the  origin  if  the  state- 
error  has  already  done  so.  However,  v°(k)  may  also  converge  to  the  origin 
prior  to  e(k).  We  may  show  this  by  algebraic  manipulation  of  already  avail- 
able equations.  Using  (2.4-5)  and  (2.1-14)  we  obtain  a recursive  equation 
for  v°(k) : 

xJ(k+l)F(k)X  (k)  _ 

v (k+1) 2—— 2 v°(k)  + xA(k+l)6p(k)  (2.4-6) 

1 + X*(k)F(k)xg(k) 

Using  (2.2-18,  2.1-10,  2.1-11),  we  obtain  a new  recursive  equation  for  e(k): 

o . 

e(k+l)  = F (k)e(k)  + G — — ^ (2.4-7) 

1 + Xg(k)F(k)*s(k) 

where 


0 1 0 ...  0, 

• • 

: o 

o ...  o 1 

a^(k)  ...  a^ (k) i 


is  a submatrix  in  the  S-subsystem  matrix  F(k).  Now  suppose  v°(k)  * 0, 


implying  X^(k)  6p(k)  ■ 0,  and  that  x^k+Jl)  remains  orthogonal  to  6p(k)  for 


e 
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all  8,  - 0.  Then  vo(k+i)  ■ 0,  dp(k+&)  ■ dp(k)  = dp*,  and  (2.4-7)  is  unforced 

e(k+l)  - Fg*  e(k) . (2.4-9) 

e(k)  will  be  asymptotically  stable  if  the  final  minimal  set  of  poles  (cor- 
responding to  the  eigenvalues  of  Fg*)  for  the  S-subsystem  are  stable,  and 
will  converge  to  the  origin  at  a rate  dominated  by  the  slowest  Eigenvalue 

of  F *.  It  is  interesting  to  consider  whether  it  is  possible  for  F * to 
o S 

have  unstable  eigenvalues.  From  our  analysis  thus  far  in  this  section, 
nothing  prohibits  such  an  occurrence.  We  will  study  this  question  in  more 
detail  in  Section  2.5. 

So  far  we  have  shown  that  dp(k)  is  a fixed  point  dp*  under  either 
of  the  two  conditions: 

1)  v°(k)  has  converged  to  0 

2)  e(k)  and  e(k)  have  converged  to  0; 

and  that  (2)  implies  (1),  but  not  the  converse.  We  now  develop  two  geometric 
interpretations  for  the  constraint  that  must  be  satisfied  by  6-p*.  The  first 
approach  is  based  on  (2.4-5),  the  second  on  a new  dynamic  equation  for  dp(k) . 

When  the  MRAS  adaptation  process  has  ceased  (v°(k)  = 0 for  all  k) , 
the  final  parameter-error  dp*  is  constrained  to  lie  successively  in  elements 
of  an  infinite  sequence  of  hyperspaces  [43],  (H(k+&)}  defined  by: 

H(k+£,)  - {dp*  e R2N+1  | xj(k+£)dp*  = 0,  l t [0,“)}  (2.4-10) 

where  k is  any  stage  after  which  v°  has  converged.  Let  (H*(k+Z)}  be  a 

2N4-1- 

sequence  of  subspaces  in  R defined  by 
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H*(k)  = R 


H*(k+£+l)  - H*(k+&)  n H(k+A) 


(2.4-11) 


where  D denotes  the  set  intersection  operator.  Clearly,  H*(k+Jt+1)  H*(k+£.). 

since  each  H(k)  is  a closed  set, 


H*  * lim  H*(k+Jt) 
Z-**> 


(2.4-12) 


is  also  a closed  set  [4].  Since  6p*  lies  in  every  H(k+Jl)  , it  also  lies  in 

every  H*(k+f.)»  and  finally  in  H*.  In  fact,  H*  completely  characterizes 
2JJ4-2  

the  subspace  in  R within  which  &p*  is  constrained  by  the  state  vector 
sequence  {^(k+Jl)  } . 

This  geometric  viewpoint  has  an  algebraic  counterpart.  Correspond- 
ing to  each  hyperspace  H(k+£)  is  a linear  functional  equation  (a  linear 
transformation  from  R^N+^  to  R) : 


«"(k+Jl)  6p*  = 0 


which  follows  from  (2.4-5).  Consider  the  linear  matrix  equation: 


(2.4-13) 


xJ(A)6p*  = 0 


(2.4-14) 


— T 

^(k+i) 
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and  J,  is  the  least  integer  l e [O,00)  such  that  p{XM(A)}  » for  all 

£.  i H.  In  other  words,  adding  more  stages  to  X^(Jl)  will  not  increase  the 
rank  of  X^(fc) • The  null-space  of  X^(&)  is  identical  to  the  subspace  H*  in 
(2.4-12): 


N{xJ(£)}  = H* 


(2.4-16) 


We  note  that,  if  X^His  full  rank,  p = 2N  + 1,  then  A/{xJ(£.) } = H*  = {0}, 
and  (2.4-14)  has  the  unique  (trivial)  solution  6-p*  = 0.  This  means  param- 
eter identification  has  been  successful. 

The  second  approach  to  a geometric  description  of  the  constraint 
subspace  for  6p*  is  based  on  the  following  autonomous  time-varying  linear 
equation  for  <$p(k) , derived  from  (2.4-5)  and  (2.1-14): 


(5p(k)-l)  - P(k)  frp(k) 


(2.4-17) 


where 


II 


I 

I 

l 

0 

G 

G 


P(k) 


2N+1 


F(k)Xs(k)  xj?k) 

1 + Xg(k)F(k)Xg(k) 


(2.4-18) 


When  v (k)  * 0 for  all  k,  6p(ld-l)  _ 6^>(k)  = 6p*.  The  following  proposition 
establishes  the  basis  for  the  subspace  in  which  6p*  is  constrained. 

Proposition  2.4.1 

Matrix  P(k)  has  an  eigenvalue  X ■ 1 of  multiplicity  m - 2N,  and 
full  degeneracy  q ■ 2N  if  and  only  if  xg(k)  j*  0. 
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Proof 


p{[Xl2N+r  p<k>  W - p{- 


F(k)xg(k)x£(k) 

1 + kg(k)F(k)3^(k) 


} 


since  F(k)  > 0 and  p{kg(k)xj[j(k)}  = 1.  Thus  A = 1 is  an  eigenvalue  of  P(k) 
with  degeneracy  q = (2N+l)-p  ■ 2N.  It  can  be  shown  that  P(k)  has  one 
eigenvalue 


~ 1 - xg(k)F(k)e(k) 

1 + kg(k)F(k)xs(k) 

where  eT(k)  = [eT(k)Io  ...  o].  Only  when  xc(k)  = 0 will  A = 1,  giving  A = 1 
a multiplicity  m = 2N  + 1.  Otherwise,  multiplicity  fpr  A = 1 is  m = 2N, 
and  degeneracy  is  full,  q = m = 2N. 

■ 

From  Proposition  2.4.1  and  (2.4-17),  6p*  is  an  eigenvector  for 
P(k),  and  since  A = 1 has  multiplicity  m = 2N  and  full  degeneracy,  2N  linearly 
independent  eigenvectors  for  A = 1 exist.  Thus  <5p*  is  constrained  at  stage  k 
to  lie  in  a (2N-dimensional)  hyperspace,  which  we  again  denote  H(k): 


H(k)  = span  (z^k)  ...  z2N(k)}  (2.4-19) 

where  (z^(k)  ...  z2N(k)}  is  the  set  of  eigenvectors  for  P(k)  associated  with 
A = 1.  Thus  6p*  is  cor  3 trained  successively  to  the  sequence  of  hyperspaces 
{H(k+£) | l > 0}. 

The  two  hyperspaces  (2.4-10)  and  (2.4-20)  are  identical;  we  show 
this  informally  with  the  following  argument: 

dim  N.{jtgX^}  * 2N,  thus  N{XgX^}  corresponds  to  a hyperspace  in 

2N+1  T'  * —i T 

R . Now,  (2.4-17)  implies  (x^x^dp*  = 0,  which  in  turn  implies  Xg(x^6p*)  = 

0 . Therefore,  x^6p*  ■ 0,  which  is  simply  (2.4-5),  the  basis  for  (2.4-10). 
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We  now  use  the  algebraic  characterization  for  the  6p*  constraint 
subspace  to  examine  conditions  on  the  input  sequence  for  &p*  to  be  constrained 
to  the  origin.  We  stated  earlier  that  (2.4-14)  has  the  unique  solution 
5-p*  = 0 if  and  only  if  X^Jl)  is  full  rank,  P = 2N  + 1.  In  order  to  examine 
when  is  full  rank  we  substitute  for  each  X^k+f.)  the  equivalent  ex- 

pression using  X^(k)  and  {u(k+l) , . . .u(k+&) } . Then  we  group  terms  into  four 
uubi.intriccs,  obtaining: 


XmW« 


‘11 


12 


X21  i X22 


(2.4-20) 


the  submatrices  X^,  X^j  X^^,  X ^ having  dimensions  Nx(N+l) , Nx(&-N)  , 
(N+l)  x (N+l),  and  (N+l)  x (&-N) , respectively. 


where 


‘11 


[[IN:0]xM(k)l[FM:FM(l)]xM(k):  '•••  : [fmVm(n)  ]xM(k) 


+ [h^(N-l)  ...  hM(0)  ] 


0 0 


0 u(k) 


'*  0 ° 

- ° u(k) 

0 u(k)  u(k+l)...  u(k+N-l) 


(2.4-21) 
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X12 


1FM1FS:'M(B)1»M(k)r  •••  :FMt‘NlF{l?M(N)lxM(k)1  + 


j t-w+i)  (N).  __  iFh  (M) :i,  (n) : ...  :h  (o)] 

M M . . M . m . .11 


0 u(k) 


X / 

u(k)  u(k+l)  ... 

u(k+I-(N+l)) 

u(k+l) 

u(k+I-N) 

u(k+N) 

u(k+l-l) 

X21  = 


u(k-N)  u(k-(N-l) ). . . u(k) 

! u(k+l) 

! u(k)  : 

u(k)  u(k+l)  ...  u(k+N) 


X22  = 


u(k+l) 


u(k+f,-N) 


u(k+N+l)  ...  u(k+£) 


where  F^,  FM(*)  and  h^*)  are  defined  in  Appendix  A and  (2.3-4). 

in  order  for  (2.4-21)  to  have  rank  p = 2N  + 1,  it  is  necessary  and 
sufficient  that  we  be  able  to  reorder  columns  so  the  first  2N  + 1 reordered 
columns  forms  a non- singular  submatrix.  However,  it  is  not  clear  what 
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necessary  and  sufficient  conditions  the  input  sequence  (u(k) u(k+£)} 

must  satisfy  in  order  for  this  general  reordering  to  be  possible. 

We  consider  instead  sufficient  conditions  for  (2.4-21)  to  have  full 
rank.  From  these  sufficient  conditions  we  are  able  to  make  some  partial  con- 
elusions  about  the  input  sequence.  The  sufficient  conditions  are  stated  in 
the  following  Lemma: 


Lemma  2.4.1 


The  matrix  X^(Jl)  is  full  rank,  p = 2N  + 1,  if: 
1)  X2^  is  nonsingular 

X, 

I 

has  rank  p £ N 


2) 


22 


3)  There  exist  N linearly  independent  columns  in 
they  span  the  N-dimensional  subspace  of  R^N+'*' 


fd 


such  that 


which  is  com- 
X 

■ 

plementary  to  the  (N+l)-dimensional  subspace  spanned  by 


:om 

X 

X 


11 
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Condition  (1)  implies  the  input  sequence  (u(k-N),  ....  u(k+N)}  must  be  such 
that  (2.4-24)  is  nonsingular.  Viewing  {u(k-N),  ...»  u(k+N)}  as  a vector 
element  e R^N+\  this  vector  "almost  always"  yields  a nonsingular  X^, 
except  when  UN  lies  on  a certain  nonlinear  hypersurface  defined  by  det 
(X^)  = 0.  For  example,  for  N = 1 and  N = 2,  respectively: 


det  (X21)  = u(k-l)  u(k+l)  - u2(k) 

det  (X^^)  ■ u(k-2)u(k)u(k+2)  + 2u(k-l)u(k)u(k+l) 

- u3(k)  - u(k-2)u2(k+l)  - u2(k-l)u(k+2) 


jpppip pmi 

P-  ] 

j 

Conditions  (2)  and  (3)  provide  little  insight  into  properties  which 

should  be  satisfied  by  the  input  sequence  (u(k) u(k+Jt)}.  However,  we 

make  the  following  remark.  In  an  effort  to  find  conditions  which  move  in  the 
direction  of  being  both  necessary  and  sufficient  to  guarantee  X^(Jl)  full  rank. 
Lemma  2.4.1  suggests  that  such  conditions  will  depend  both  on  the  input  se- 
quence and  on  the  M-subsystem  parameter  vector  p^.  The  latter  dependence 

A 

follows  from  definitions  for  F^,  FM(*)>  and  h^(*)*  Thus,  necessary  and  suf- 
ficient conditions  on  the  i*put  sequence  to  guarantee  X^Jl)  full  rank  must  be 
expressed  in  terms  of  p^,  and  would  be  dependent  upon  the  specific  plant  being 
identified.  This  is,  perhaps,  why  all  previous  investigations  into  this  prob- 
lem have  focused  on  sufficient  conditions  only,  at  a level  which  depends  only 
on  the  order  of  the  MRAS,  or  number  of  parameters,  but  independent  of  any 
plant  parameter  values. 

We  conclude  this  section  with  further  consideration  of  the  conditions 
proposed  by  Lion  and  by  Kim-Lindorf f , described  earlier,  which  deal  with  the 
frequency  content  of  periodic  inputs  sufficient  to  guarantee  parameter-error 
convergence  to  0.  When  the  Lion  and  Kim-Lindorf f conditions  are  each  applied 
to  the  ELMA,  the  same  sufficient  frequency  content  limit  q*  is  obtained  in 
each  case.  In  Lion's  case, 

* 2N  + 1 „ . . 

q*  1 — ^ = N + 1 

In  the  Kim-Lindorff  case,  we  consider  the  1-0  Delay  state  variable  representa- 
tion (2.2-8)-(2.2-12),  and  conclude  that  R = 2N  + 1 from  the  Nth  row  of 
[6f:6G] . Thus, 
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q*  = min  {q|q  £ R/2}  = N + 1 

From  simulation  experience,  however,  (see  Examples  2.2-1  and  2.2-2)  it 
appears  that  q*  = N may  always  be  sufficient  for  parameter-error  convergence. 

We  consider  now  two  more  examples  to  substantiate  this  conjecture,  first  for 
N = 3,  and  then  for  N = 4.  ' 

Example  2.4-1,  N = 3 
M-subsystem: 

a)  3 poles  at  z = (0.9,  0.4,  -0.5)  j 

b)  1 zero  at  z = (0.8) 

N 

c)  b =0 

o 

d)  b 2 chosen  to  provide  unity  gain  on  |z|  = 1. 

Thus  p^  = [-0.18,  0.29,  0.81  -0.36,  0.45,  0,  0] 

Initial  S-subsystem: 

a)  poles  at  z = (0,  0,  0) 

b)  S-transfer  function  numerator  = M-transfer  function 
numerator . 

Thus  pT(0)  = [0,  0,  0|  -0.36,  0.45,  0,  0] 

A-subsystem: 

a)  F(0)  = diag  [1000 1000] 

b)  c ( t)  = - a(t) 

Input  u(k) : 

sine  with  (amplitude,  period)  = (1,20)  + (3,7)  + (8,4) 
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Final  S- subsystem: 

PT(10)  = [-0.18,  0.291,  0.799:  -0.361,  0.451,  10-5,  10-5], 
corresponding  to  (poles,  zeroes)  = (0.9,  0.4, 
-0.50i:  0.8,  38500) 

Simulation  results  are  shown  in  Figures  2.4-1.  We  use  a sinusoidal  input 

u(k)  with  N = 3 frequencies,  and  find  that  parameter  identification  is  suc- 

cessful  to  within  10  after  about  3 seconds  (30  discrete  stages).  The 

figure  also  indicates  that  e(k)  is  converging  to  the  origin.  As  was  shown 

earlier  in  this  section,  after  v°(k)  has  converged,  the  motion  of  e(k)  is 

governed  by  the  final  poles  of  the  S-subsystem.  Since  one  pole  is  at  z = 

0.9,  a relatively  slow  transient  for  e(k)  should  be  anticipated.  Parameter 

trajectories  for  ( a, , b ) are  not  shown  in  the  figure. 

1 o u 

Example  2.4-2  N = 4 
M-subsystem: 

a)  4 poles  at  z = (0.9,  0.4,  -0.5,  0) 

b)  1 zero  at  z = (0.8) 

c)  b =0 

o 

d)  b^  chosen  to  provide  unity  gain  on  |z|  =1 

Thus  p1  - [0,  -0.18  , 0.29  , 0.8*.  -0.36  , 0.45  , 0 , 0 , 0] 

M • 

Initial  S-subsystem: 

a)  poles  at  z = (0,  0,  0,  0) 

b)  S-transfer  function  numerator  = M-transfer  function 
numerator . 

Thus  pT(0)  - [0,  0,  0,  O:  -0.36,  0.45,  0,  0,  0] 
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Figure  2.4-1  E1MA  Example  2.4-1 
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A-subsystem: 

a)  F(0)  = diag  [1000,  ....  1000] 

b)  c ( t)  * -a(t) 

Input  u (k) : 

sine  with  (amplitude,  period)  ■ (1,  20)  + (3,13)  + (7,  7)  + 
(12,  4). 

Final  S-subsystem: 

pT(10)  - [-0.8  x 10'3,  -0.179,  0.29,  0.8:  -0.36,  0.45, 

-6  -4  -4 

10  , -10  ,10  ],  corresponding  to  (poles,  zeroes)  = 

(0.9,  0.402,  -0.498,  -0.46  x lO-2!  0.8,  174.,  -174.) 

Simulation  results  are  shown  in  Figure  2.4-2.  Ue  use  a sinusoidal  input 

with  N = 4 frequencies,  and  find  that  parameter  identification  is  successful 
-3 

to  within  10  after  about  3 seconds.  Parameter  trajectories  for  (a2,  a^» 
b2,  b^,  bQ)  are  not  shown  in  the  figure.  ^ 

2.5  Analysis  of  Instability  in  the  Landau  MRAS 
Identification  Algorithm 

In  this  section  we  will  examine  the  instability  inherent  in  the 

Landau  MRAS  or  ELMA  identification  algorithm.  Although  the  algorithm  is 

structured  in  accordance  with  Popov's  hyperstability  conditions,  it  is  not 

generally  possible  to  satisfy  the  Popov  necessary  and  sufficient  condition. 

T 

on  the  linear  subsystem  of  Figure  2.1-3.  A vector  C = [c^  ...  c^]  must  be 


chosen  so  that  the  linear  subsystem  transfer  function 
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Coefficient 
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Figure  2.4-2  EIMA  Example  2.4-2 


N . " „ N-i 

z + I c. z 

1=1 

N 

N “ N-i 
z “ L a^.  z 

i=l  “i 


Z (SL,  +C.)ZN'2 
1=1  ^1  1 

N J N-i 
z - i a,  z 

i=l  ”i 


(2.5-1) 


T 

is  positive-real.  However,  since  the  M-subsystem  parameter-vector  a^  = 

[a^  ...  a^  ] is  unknown,  it  is  not  possible  to  evaluate  the  positive-real 
property  of  (2.5-1),  no  matter  what  choice  is  made  for  c.  A desirable  choice 
for  c would  be 

c “ - V (2.5-2) 

since,  in  this  case,  H(z)  = 1,  clearly  positive-real.  Landau  suggests  [35] 

A 

a strategy  where  first  one  obtains  an  estimate  for  a^  denoted  a^  by  a method 
with  relatively  low  computational  complexity,  such  as  generalized  least- 
squares  [2],  and  then  uses  this  estimate  both  for  the  initial  identification 
model  parameter  value,  a(0),  and  for  c;  that  is,  assign 


= - = - a(0) 


(2.5-3) 


A natural  extension  of  this  idea  (and  one  adopted  in  the  ELMA  identification 
algorithm)  is  to  specify  that  c track  the  value  of  a(k);  that  is. 


c(k)  = - a(k) , 


(2.5-4) 


the  strategy  being  for  c(k)  to  approach  -a^  as  the  MRAS  identification  pro- 
cedes  to  its  successful  conclusion.  The  success  of  this  strategy  will  depend 
on  the  starting  point  a(0),  and  may  present  a problem  particularly  in  those 


practical  situations  where  It  Is  not  feasible  to  obtain  an  adequate  initial 
estimate  for  a^. 

Choosing  c defined  by  (2.5-4)  does  not  guarantee  hyper stability 
for  the  MRAS.  Thus  we  would  like  to  find  the  region  in  Rn  in  which  c(k)  must  ' 

remain  in  order  to  guarantee  a hyperstable  MRAS.  This  MRAS  in  the  Popov 
configuration  no  longer  conforms  to  the  assumption  that  the  linear  subsystem 
is  time -invariant.  We  might  consider  viewing  (2.5-1)  as  a linear  time-varying 

■j 

system,  and  then  determine  under  what  condition  on  c(k)  (2.5-1)  satisfies 

j 

the  time-varying  version  of  the  Discrete  Positive  Real  Lemma  [35].  Based 
on  simulation  results,  we  expect  this  region  for  c(k)CRn  is  larger  than  the 
region  constraining  c*  constant.  For  c«  constant  it  is  required  that,  as 
part  of  the  positive-real  constraint  on  (2.5-1),  the  numerator  polynomial 
of  (2.5-1)  have  all  its  roots  inside  the  unit  circle.  However,  for  the  case 
(2.5-4)  it  has  been  observed  that,  during  a finite  number  of  stages  of 
stable  MRAS  operation,  the  roots  of  the  numerator  polynomial  may  lie  outside 
the  unit  circle. 

We  begin  this  section  by  discussion  the  difficulty  of  testing 
the  positive-real  condition  for  (2.5-1).  TVo  possible  approaches  for 
determining  whether  a transfer  function  is  positive-real  are  suggested. 

These  methods  are  of  use  only  for  analysis,  where  we  shall  assume 
the  plant  parameter  vector  aj^  is  "known."  Thus,  these 


68 


tests  would  not  be  useful  in  an  on-line  Identification  context.  We  then 
direct  our  attention  to  the  adaptation  gain  equation  (2.1.15)  to  observe 
what  effect  MRAS  instability  has  on  F(k) . Finally,  we  consider  whether  it  is 
possible  for  the  Landau  MRAS  identification  problem  to  be  restructured  or 
reformulated  in  a way  which  retains  hyperstability  but  avoids  the  potential 
instability  problems  related  to  the  positive-real  condition. 

Definition  2.5.1  [16]  A discrete-time  transfer  function  H(z)  = 

q(z)/p(z) , where  q is  order  m,  p is  order  n,  is  strictly  positive-real  if: 

1)  All  roots  of  p(z)  and  q(z)  lie  strictly  within  the  unit  circle 
of  the  complex  z-plane. 

2)  HCe^40^)  + H(e  ^wT)  > 0 for  all  u>  £ 0,  where  T is  the  fixed 
time  interval  between  discrete  stages. 


Since  RefHCe^10^)]  is  an  even  function  of  u>,  and  Im[H(e^W^)]  is  odd,  condition 
(2)  may  be  expressed  Re[H(e'^W'r)]  > 0.  For  example,  when  N ■ 1,  the  transfer 
function  (2.5-1)  becomes 


z + c. 


H(z) 


z - 


aM. 


1 + 


\ + C1 

^V 


(2.5-5) 


and  may  be  shown  to  be  positive-real  only  if  |a^  | < 1 and  |c^|  < 1.  In  this 
simple  case,  the  choice  of  c is  independent  of  a^.  When  N t 2,  the  analysis 
becomes  more  complicated,  and  the  positive- real  region  for  c,  given  a^,  will 
depend  on  a^.  A general  set  of  inequalities,  relating  c to  a^,  which  guar- 
antees (2.5-1)  is  positive-real  is  not  known  to  have  been  derived.  In  fact, 
condition  (2)  applied  directly  generates  a continuum  of  inequalities  involv- 
ing c and  a^,  parametrized  by  u.  However,  the  continuum  may  be  reduced  to  a 
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finite  set  when  one  recognizes  the  polynomial  form  of  the  inequalities.  For 
example,  when  N = 2,  condition  (2)  becomes 


y(c,aM,x)  = 2(c2-a2)x  + [CjU-a^-a  ]x 


+ [l+(a2-c2)-(a1c1+a2c2> ] > 0, 


(2.5-6) 


a quadratic  polynomial  in  -1  i x = cos  0)T  $ 1.  By  the  nature  of  the  quadratic 
form  on  the  interval  [-1,  1],  we  can  state  conditions  for  the  quadratic  to 
be  positive  on  this  interval: 

1)  yCc.a^+l)  > 0 

2)  y(c,a^,-l)  > 0 

3)  If  x min  5 1,  then  y min  * (c,  a^,  x min)  > 0,  where  x min  is 

the  value  of  x that  minimizes  yCc.a^,  ) for  fixed  c,  a^. 

Expressed  in  terms  of  c and  a^,  and  taking  into  account  the  stability  con- 
straint on  aM>  these  three  conditions  may  be  reduced  to: 
l')  1 + c±  + c2  > 0 

2)  l-c^+c2>0 

, a (l+c9)  - c. (1-a.) 

3 ) If  -1  < — - — * 1, 

4(c2-a2> 


-[c^l-a^  - a1(l+c2)]‘ 
8(c2-a2) 


+ [l+(a2-c2)  - (a^+a^)  ] > 0 


The  first  two  conditions  are  seen  to  reduce  to  a a subset  of  the  numerator, 
q(z),  stability  constraint  Inequalities.  Even  for  this  low  order  case,  N = 2, 
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i 

the  inequality  (3  ) cannot  be  solved  explicitly  for  c in  terms  of  a^.  For 
arbitrary  N,  although  the  problem  has  a theoretical  solution  it  is  not  a 
simple  matter  to  select  c,  even  when  is  known. 

The  only  feasible  approach  known  to  this  author  for  analyzing  the 
positive-real  condition  of  a transfer  function  is  to  apply  a test  to  specific 
numeric  examples  to  determine  whether  the  example  violates  conditions  speci- 
fied by  the  test.  This  approach  would  provide  a trial- and-error  method  for 
examining  the  c-sensitivity  question  stated  above,  and  would  require  that 
the  plant  parameters  be  assumed  known  during  the  analysis. 

Two  different  tests  can  be  proposed  in  this  context.  The  first 
involves  the  possible  application  of  the  discrete  positive-real  lemma  [16]  to 
the  given  transfer  function,  requiring  first  that  a state-variable  realiza- 
tion for  (2.5-1)  be  chosen.  This  lemma  gives  necessary  and  sufficient  con- 
ditions for  which  the  transfer  function  is  strictly  positive-real  in  terms  of 
a set  of  matrices  satisfying  three  matrix  equations . One  approach  in  applying 
this  lemma  would  be  to  develop  an  iterative  algorithm  based  on  the  three 
matrix  equations  which  would  converge  to  matrices  satisfying  the  conditions 
of  the  lemma  if  the  transfer  function  is  positive-real,  and  would  fail  to 
converge  otherwise.  This  author  is  not  familiar  with  such  an  algorithm;  the 
development  of  one  could  be  the  subject  of  future  research  in  this  area. 

The  second  approach  involves  application  of  algebraic  criteria 
obtained  by  Siljak  [55]  for  positive-realness  of  discrete-time  transfer  func- 
tions. The  Siljak  test  is  attractive  because  it  involves  finding  the  loca- 
tions of  roots  of  two  polynomials  with  respect  to  the  unit  circle.  The  test 
is  formulated  recursively,  and  although  the  second  stage  of  the  test  is  com- 
plicated, it  is  amenable  to  computation. 
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The  MRAS  exhibits  unstable  operation  by  the  fact  that  the  output 
error  e(k)  becomes  unbounded.  This  implies  that  y (k)  is  unbounded,  and 
thus  the  poles  of  the  S-subsystem  (determined  by  a(k))  must  lie  outside 
the  unit  circle.  Once  this  unstable  condition  occurs  in  the  MRAS,  the 
adaptation  gain  F(k)  responds  in  such  a way  that  adaptation  of  the  S- 
subsystem  poles  is  effectively  halted.  We  will  analyze  the  adaptation 
gain  dynamics  with  the  goal  of  obtaining  insight  into  a process  in  the  MRAS 
algorithm  which  is  linked  with  the  unstable  condition. 

The  adaptation  gain  equation  for  the  MRAS  identification  algorithm 
is : 


[F(k)x  (k)][F(k)x_(k)]T 

F(k+1)  = F (k) 1 5 

1 + xJ(k)F(k)x  (k) 


= FT(0)  > 0,  it  can  be  shown  that 


0 < F (k+1 ) < F(k) , 


(2.5-7) 


(2.5-8) 
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that  is,  F(k)  remains  positive  definite  for  all  k,  independent  of  x^(k),  and 
approaches  F = 0 in  the  limit  as  k + “,  independent  of  xc(k),  if  x = 0 for 

b O 

only  a finite  number  of  stages. 

We  shall  partition  F(k)  to  correspond  with  xg(k)  as  partitioned  in 
(2.1-2).  Thus 


F^k) 


Fi2(k) 


F(k)  = 


(2.5-9) 


F12(k)  F2(k) 


where  F^  is  N x N,  F2  is  (N+l)  x (N+l),  and  the  gain  equation  2.5-7  may  be 
decomposed  to: 


Fx(k+1)  * Fx(k)  - 


T X T 

Flxsxs  F1  + F1XSXS  F12  + F12XSXS  F1  + F12XSXS  *12 

1 + -sTfxs 


(2.5-10a) 


F12(kfl) 


F12(k) 

_ F1XSXS  F12  + F1XSXS  F2  + ?12XSXS  F12  + F12XSXS  F2 
1 + x^Fxs 


(2.5-10b) 


F2(k+1) 


F2(k) 


T T T T 

F12  XSXS  F12  + F12  XSXS  F2  + F2XSXS  F12  + F2XSXS  F2 

1 + xsFxs 


(2.5-10c) 
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where  the  time- index  k has  been  suppressed  for  convenience.  From  the  parameter 
adaptation  equation  (2.1-14),  and  using  the  above  partitioning  of  F(k),  that 
part  of  the  parameter  adaptation  affecting  the  S-subsystem  poles  (i.e., 
a(k))  may  be  written: 


a(kfl)  = a(k)  + 


F1(k)x°(k)  + F12(k)xJ(k) 


-T, 


1 + xg(k)F(k)xg(k) 


(2.5-11) 


When 


instability  occurs,  ||  Xg(k)||  , corresponding  to  the  S-subsystem  N- 


stage  delayed  output  sequence,  becomes  unbounded,  while  ||  x^(k)||  , correspond- 


ing to  the  (N+l)-stage  delayed  input  sequence,  remains  finite. 


Proposition  2.5.1 


In  the  limit  as  ||  x°(k)||  °°,  ||  x^(k)||  bounded. 


F(k+1) 


0 F2(k+1) 


We  will  not  prove  this  proposition  formally,  but  will  rather  develop  an  in- 
formal argument  using  (2.5-10)  for  the  case  N = 1.  Since  (2.5-10)  is  non- 
linear in  F,  it  is  not  a straightforward  matter  to  perform  a general  analysis 


of  this  system  of  dynamic  matrix  equations.  When  N = 1,  x°  and  F are 

w 1 


-1 


scalars,  Xg,  F^2,  and  F2  are  2 x 1,  1 x 2,  and  2x2,  respectively.  As 


|x°(k)|  terms  quadratic  in  x°  dominate,  and  (2.5-10)  becomes: 


f12(x°)2 

F,  (k+1)  = Fn  (k)  - - A = F1  - FL  = 0 (2.5-12a) 


V*s>2 
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F12(k+1) 


Fi2(k)  - 


F1^XS^  F12 
. o.2 

W 


F - F 
12  12 


(2.5-12b) 


I 

E 

C 

E 


F2(k+1) 


F2  (k) 


(x°)^F^  F 
( SJ  *12*12 

/ <\2 

W 


F2  - 


T 

F F 
12  12 


(2.5-l2c) 


Proposition  2.5.1,  in  conjunction  with  (2.5-11),  leads  to  the  observation 
that  when  the  MRAS  exhibits  instability,  S-subsystem  pole  adaptation  is 
effectively  halted.  As  ||x°(k)||  becomes  unbounded,  one  observes  in  simula- 
tions of  (2.5-7)  that  F^  and  F^2  converge  rapidly  to  0,  preventing  a(k) 
from  possibly  returning  to  a region  corresponding  to  stable  poles.  We 
should  emphasize,  however,  that  this  property  of  the  dynamic  adaptation 
gain  matrix  is  not  the  cause  for  MRAS  instability.  We  have  suggested  that 
this  instability  is  due  to  the  loss  of  positive-realness  in  (2.5-1). 

What  is  significant,  however,  is  that  the  MRAS  having  entered  into  an 
unstable  mode  of  operation,  coupled  with  the  property  of  F(k)  under  such 

conditions,  creates  a situation  which  appears  to  reinforce  the  instability. 

T 

One  should  examine  whether  using  a fixed  gain  F * F >0  might  improve  the 
stability  properties  of  the  MRAS  when  using  (2.5-4). 

We  have  concluded  that  potential  instability  of  the  Landau  MRAS 
is  due  to  the  difficulty  of  selecting  the  vector  c so  that  (2.5-1)  is 
guaranteed  positive-real  over  the  entire  time  of  MRAS  operation.  The 
vector  c was  introduced  into  the  MRAS  structure  in  the  definition 


I 


r? 
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v (k)  = c e(k)  + e(k) 


(2.5-13) 


In  concluding  this  section  we  will  examine  this  equation,  first  considering 
why  it  is  necessary  to  include  a term  involving  the  output-error.  Then  we 
will  consider  reformulating  the  MRAS  structure  in  such  a way  as  to  eliminate 
the  term  involving  the  unknown  vector  c.  By  doing  this,  we  eliminate  the 
cause  for  potential  instability,  and  obtain  a "robustly"  hyperstable  MRAS, 
one  whose  stability  does  not  depend  on  the  value  of  a design-variable. 


The  output  error  term  e(k)  = yM(k)  - yg(k)  in  (2.5-13)  may  be  ex- 


pressed: 


e(k)  = a^e(k)  + Xg(k)  6p(k+l)  = a^eCk)  + w(k) 


and  so  the  Popov  linear  subsystem  in  Figure  2.1-3  appears  as 


e(k+l)  = FMe(k)  + G2  w(k) 


(2.5-14) 


v(k)  = (a^c)  e(k)  + w(k) 


(2.5-15) 

(2.5-16) 


Thus  we  see  that  the  term  e(k)  in  (2.5-13)  is  responsible  for  introducing 
the  feedforward  element  in  the  Popov  linear  subsystem.  Excluding  e(k)  in 
(2.5-13)  would  yield  (2.5-16)  in  the  form: 


T 


v(k)  = (a^+c)  e(k) 


(2.5-17) 


We  make  use  of  the  discrete  positive-real  lemma  [16]  in  the  following  pro- 


position 


Proposition  2.5.2 


A SISO  positive-real  linear  system  of  the  form 


x(k+l)  = Ax(k)  + Bki(k) 


y(k)  = e x(k)  + Du(k) 


necessarily  contains  a non- zero  feedforward  element  D t 0. 


Proof : According  to  the  discrete  positive- real  lemma,  necessary 

and  sufficient  conditions  for  a linear  system  to  be  positive-real  are  that 
T 

there  exist  P = P > 0,  L,  and  W such  that 

T T 

1)  A PA  - P = -LL 

2)  ATPB  = C - LW 

3)  WTW  = D + DT  - BTPB 

If  D = 0,  condition  (3)  cannot  be  satisfied,  since  WTW  - 0,  and  -BTPB  < 0 
since  P > 0. 

■ 

This  proposition  implies  that  v(k)  necessarily  contains  a linear 
term  in  e(k),  given  the  nonlinear-subsystem  output,  z(k),  in  order  for  the 
linear  subsystem  to  contain  a feedforward  element.  This  is  in  contrast  to 
the  continuous- time  hyperstable  parallel  MRAS  structure,  where  it  is  possible 
for  a stable  MRAS  to  have  an  adaptation  equation  of  the  form 


p(t)  * F xg (t) [C  e(t) ] , 


i.e.,  simply  using  the  state  error  as  an  input  to  the  adaptation  subsystem. 
In  the  discrete-time  case,  the  term  e(k)  must  be  included  as  well.  One 
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might  interpret  this  difference  by  saying  that  the  parameter  change  Ap(k)  = 

F xsv(k)  rust  take  into  account  the  effect  of  parameter  change  on  the  system 
at  the  time  stage  k by  including  e(k)  as  part  of  v(k).  In  the  limit,  as  the 
sample  time  T -*■  0,  e(k)  -*•  e(k-l),  and  thus  e(k)  conveys  no  information  in 
addition  to  the  state-error  vector,  hence  is  not  necessary  in  the  continuous- 
time problem.  (Note  that  e(k-l)  is  an  element  of  the  state-error  vector, 
as  defined  in  2.1-5.) 

Discrete-time  MRAS  identification  algorithms  have  been  presented 
[19,  29]  which  are  asymptotically  stable  in  the  state-error  and  which  con- 
tain no  design  parameters  depending  on  the  unknown  plant  parameters  pM> 
Ionescu  and.Monopoli  achieve  a stable  parallel  MRAS  design  by  introducing  an 
augmented  error  signal  generated  by  an  augmentation  filter.  This  approach 
seems  to  be  related  to  the  augmented  state-variable  filter  approach  for  the 
continuous-time  MRAS  problem  [40].  Kudva  and  Narendra  propose  a stable 


series-parallel  MRAS  which  is  designed  from  the  Liapunov  viewpoint,  without 


requiring  knowledge  of  plant  parameters,  and  requiring  only  the  state-error  e(k)  = 


x^k)  = xg (k)  as  an  input  to  the  adaptation  subsystem.  The  Landau  MRAS 
structure  on  which  we  have  focused  thus  far  can  also  be  reformulated  into  a 
series-parallel  context,  resulting  in  a hyperstable  system  with  an  adaptation 
equation 

p(k+l)  = p(k)  + F(k)  ^(k)  e(k)  (2.5-19) 

The  equivalent  Popov  structure  for  this  series-parallel  MRAS  is  shown  in 
Figure  2.5-1.  The  linear  subsystem  is  always  positive-real  without  having  to 


FP-S3S7 


Figure  2 


and  G2  are  defined  in  Section  2.1 


5-1  Equivalent  Popov  Structure  for  Landau  Series-Parallel  MRAS 
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introduce  additional  design  parameters,  as  was  necessary  in  the  parallel 
MRAS  case,  eliminating  the  problem  of  potential  instability. 

This  stability  design  advantage,  however  is  offset  by  a disadvantage 
relating  to  measurement  noise.  The  adaptation  process  (2.5-19)  differs  from 
the  equivalent  equation  for  the  parallel  MRAS  in  that  it  uses  plant  output 
measurements  rather  than  identification  model  output  measurements  (xM  vs . 

Xg).  In  a practical  situation,  the  former  is  assumed  to  be  relatively  more 
noise- contaminated.  The  effect  of  this  noise  is  to  produce  a biased  param- 
eter estimate  in  the  series-parallel  case,  in  contrast  to  an  unbiased  param- 
eter estimate  in  the  parallel  case.  In  the  identification  of  a real  heat 
exchanger  plant,  it  was  demonstrated  [15]  that  this  is  indeed  what  occurs 
when  comparing  the  two  MRAS  structures. 

Thus  we  conclude  that  the  problem  of  potential  instability  of  the 
Landau  parallel  MRAS  Identification  algorithm  may  be  confronted  from  one  of 
several  directions.  One  may  retain  the  given  parallel  structure,  and  endeavor 
to  select  the  design  vector  c in  a way  which  retains  long-term  MRAS  hyper- 
stability. It  may  be  necessary  to  use  a pre-identifier,  such  as  a least- 
squares  method  or  a series-parallel  MRAS,  to  initialize  the  parallel  MRAS. 
Another  direction  would  be  to  consider  using  only  the  series-parallel  MRAS 
structure.  The  cost  in  taking  this  approach,  however,  is  accepting  noise- 
biased  parameter  estimates . 
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2.6  Discontinuous  Adaptation  in  MRAS 

Up  to  this  point  we  have  considered  continuous  parameter  adaptation, 
that  is,  adaptation  of  the  type  (2.1-14)  which  is  a continuous  function  of 
the  variables  xg(k)  and  v(k)  . In  the  context  of  the  MRAS  control  problem, 
our  goal  is  not  generally  to  achieve  convergence  of  6p(k)  = PM  “ p(k)  to 
the  origin  in  parameter-error  space,  where  p^  corresponds  to  reference  model 
parameters,  p(k)  to  the  controlled  plant  parameters , but  rather  to  provide 
only  state-error  convergence  to  the  origin.  Thus,  continuous  adaptation 
(i.e.,  the  capability  for  p(k)  to  attain  the  value  p^)  may  not  be  necessary. 

It  may  be  desirable  to  consider  simpler  adaptation  schemes  which  require, 
for  example,  less  extensive  multiplication  hardware  (or  multiplication 
processing  in  a digital  computer)  than  is  necessary  for  the  continuous  param- 
eter adaptation  algorithm  (2.1-14). 

In  this  section  we  will  consider  a discontinuous  approach  to  stable 
single-input  single-output  MRAS  control  systems  design.  This  viewpoint  has 
been  studied  primarily  by  researchers  in  the  USSR  [53,  65,  66,  26,  27,  28], 
but  some  research  activity  in  this  direction  has  also  occurred  in  the  U.  S. 

[45,  38].  In  this  previous  work,  Liapunov- techniques  are  used  both  to  ob- 
tain and  to  analyze  stability  conditions  for  discontinuous  adaptation  algorithms 
for  MRAS  control.  A recent  comprehensive  treatment  and  analysis  is  [8].  One 


feature  of  the  discontinuous  adaptation  MRAS  to  be  studied  here  is  that  the 


plant  parameter  vector  p(k)  may  take  on  one  of  only  a finite  set  of  values, 

the  order  of  that  set  for  the  ELMA  (see  Section  2.2)  being  2^N+^,  where 

- „ 2N+1 

p e R 


1 
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First  we  present  the  discontinuous  -parameter  adaptation  MRAS  design 
for  continuous-time  systems,  proving -stability  if  a set  of  assumptions  are 
satisfied.  This  problem  is  then -transformed  to -the  signal  synthesis  form, 
where  a discontinuous  control  signal r. rather  than  parameter  adaptation,  is 
used  to  achieve  the  asymptotic  stability -of -the -state-error.  Two  classes 
of  signal  synthesis  are  considered,  which -we  shall  call  "direct"  and 
"augmented."  In  both  the  parameter  adaptation  and  signal  synthesis  forms  we 
adopt  notation  for  the  adaptation  process -presented  in  [60]. 

Next  we  introduce  a measurable  disturbance-rejection  problem 
developed  in  [60].  The  control  design  is  of  the  variable  structure  or  dis- 
continuous type.  We  find  that  the  solution  for  this  problem  is  identical  in 
form  to  the  discontinuous  signal  synthesis  MRAS  design,  and  show  that  this 
MRAS  may  be  designed  without  resorting  to  the  Liapunov  approach.  Sliding 
mode  properties  of  this  MRAS  are  discussed. 

Finally,  the  above  developments  are  applied  to  the  discrete- time 
Landau  MRAS  or  EIMA.  It  is  proved  through -hyperstability  considerations 
that  discontinuous  parameter  adaptation  provides  stability  if  a set  of 
assumptions  are  satisfied.  The  concept  of  .sliding  on  the  discontinuous 
surface  for  the  discrete-time  case  must  not  be  viewed  literally,  as  was  done 
for  the  continuous- time  MRAS. 

Discontinuous  Parameter  Adaptation 

Consider  the  continuous  time  single-input  single-output  reference 
model  and  controlled  plant,  respectively: 
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E 

C 


m 

E.  b, 


,m-i 


»M<S> 


M 


1=0  n-mfi 

Sn  - E a S11"1 
1=1  ”i 


(2.6-1) 


HS(S) 


m 

E bs 

1=0  n-m+i 


,m-i 


Sn  - E a„  Sn_i 
1=1  bl 


(2.6-2) 


where  we  assume  m £ n. 


We  may  represent  the  two  systems  in  a non-minimal  state-variable  form  as: 

,(m) 

(2.6-3) 


Vm  + bmv 

yM  = ^ ® * * * 0] 


xs  = As(t)xs  + Bs(t)  V 


yg  = [1  0 ...  0]  xg 
T 

n • 1 • 


(m) 


(2.6-4) 


F(l) 


m 
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bt  - [B^-bIi  - [0  ...  0 b Jo  ...  01] 

l 4 n-m* 

? m taT  1,11  - [*n  •"  al-,bn  •"  ‘Wl-'Vi.1 


with  the  appropriate  subscripts  (M  or  S)  attached. 


(2.6-7) 


(2.6-8) 


Note  that  this  state-variable  structure  is  the  continuous- time  analog  to  the 


1-0  Delay  structure  in  Section  2.2.  By  defining  an  error  vector: 


*M“  XS 


(2.6-9) 


the  dynamic  error  system  becomes: 


e = ^e  + B2(6pT(t)xs(t)) 
where  6?(t)  = ~ PS(t) 


x_T(t)  = [x*jv(m)] 


(2.6-10 


(2.6-11) 


(2.6-12) 


Note  that  Xj^  and  x°  are  not  the  states  of  the  reference  model  and  plant. 

Theorem  2.6.1  There  exists  a vector  c e Rn  and  matrix  A = diag  [a^] 
such  that  the  parameter  adaptation  algorithm: 


Ps(t)  = Pg(0)  + A sgn  (xg)  sgn  (c  e) 


T 

1 1 
0 

i 0 

: o 


provides  an  asymptotically  stable  MRAS  control  system;  that  is,  lim  e(t)  = 0. 

t-x» 

The  notation  sgn  (Xg)  means 

— — x 

[sgn  (xgl)  ...  sgn  (Xg  )]  . 

n+nH-1 

T T 

Proof  Consider  the  scalar  function  V = e Pe,  P = P >0.  Then 
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V - eT(Ajp  + PA^e  + 2(cTe)  0^6p(t)) 


■ - eT  Q e + 2(cTe)(xg6p) 

where  Q > 0,  and  c is  the  last  column  of  P.  Since  is  assumed  stable, 

•p  j 

such  a positive-definite  pair  (P,Q)  exists.  Thus,  if  (c  e)Xg  6p)  can  be 
shows  to  be  negative,  V is  a Liapunov  function,  and  e = 0 is  asymptotically 
stable.  Substituting  the  algorithm  for  pg(t),  we  obtain 

T r T ®fmfl_  _ T 

(c  e) (x  6p)  = | c e | E |x  | (6p  (0)  sgn(x  ) sgn  (c  e)  -a.) 

i=l  l 1 Si  1 

If  - |6p^(0)|,  then  V ^ 0,  with  equality  possible  only  at  e = 0 g 

Thus,  such  a discontinuous  adaptation  algorithm  guarantees  asymptotic 

model-following  if  we  can  bound  each  element  of  the  initial  parameter  error 

6p(0).  This  condition  implies  that  one  must  know  bounds  on  the  unknown 

plant  parameters,  p^  - p^O)  - p.^  . This  design  may  exhibit  a stable 

min  ^max 

sliding  motion  on  the  surface  o = c e = 0.  Sliding  is  a property  whereby 
the  system  trajectory  e(t)  reaches  the  discontinuity  surface  a = 0 and  then 
moves  along  that  surface  while  the  control  variable  sgn (0)  oscillates  at  a 
high  frequency.  Necessary  and  sufficient  conditions  for  sliding  to  exist  are 


sgn  (a)  = - sgn  (a) 


(2.6-13) 


or  equivalently. 


a a < 0 


(2.6-14) 


In  this  case 

• 7>  < t t n I m+l  

o s C e » c Y + cn  XgfipCO)  - cn  E ajxg  |sgn  fa) 

i=l  ai 

Thus  (2.6-14)  will  be  satisfied,  and  sliding  motion  will  occur  when: 


a)  a =*  0 

T 

n+nrt-1  _ c JLfl 

b)  ai^sJ  >— — +Vp(0) 

i=l  i n 

in  some  neighborhood  of  a = 0. 


If  sliding  motion  occurs,  the  trajectory  for  the  error  vector  is  determined 
by  the  differential  equation 


o 1 o ...  o 


o ...  o 1 

-c1/cr,...c  i /c 
1 h n-1  n 


where  e = [e,  . . . e , ] 
l n-1 


This  result  is  easily  obtained  when  one  recognizes  that  during  sliding, 

T 1 n-1 

a = c e = 0,  or  e = - — E c.e. . It  may  be  shown  [7,  8]  that  for  c chosen 

n i=l 

as  the  last  row  of  P,  (2.6-16)  is  asymptotically  stable,  i.e.,  stable  sliding 
will  occur. 


Using  the  notation  in  [60],  the  discontinuous  parameter  adaptation 


algorithm  may  be  expressed: 


P(t)  = p(0)  + ¥ 


(2.6-18) 


U 
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The  system  matrices  Ag  and  Bg  are  assumed  now  to  be  constant  during  the 
operation  of  the  MRAS . For  the  direct  signal  synthesis  case,  adaptation  will 
be  achieved  directly  through  u(t).  In  the  augmented  case  , we  define 

u(t)  = v ( t)  + Au(t)  (2.6-24) 


with  adaptation  being  achieved  through  Au(t).  The  system  error  equation  for 
the  direct  and  augmented  control  cases  are,  respectively: 


where 


® = V + B2{6^S  + bM*M  - bSXS  - bS  u(m)} 

n-m 

e = + B2't^xs  + bTx^  - “ bS  Au^} 

n-m 

. 1 1 
Au  = Kg  - Xm 


(2.6-25a) 

(2.6-25b) 

(2.6-26) 


Equating  (2.6-25)  with  (2.6-10),  we  obtain: 


where 


o 

«i 

if  x°  o 

Si 

O 

-°i 

if  x°  0 

\ 

Si 

r 

1 

J , -l 

lf  V 

^ 1 

i 

1 

-1 

V 

if  ^0 

1 

i 

> 0 
< 0 

> 0 
< 0 


i -1  1,  ...,  n 


i 1,  ...,  iq  1 


(2.6-27a) 


umi  m»mm 


[ 0tsi  if  \°>0 

*J  - 

1 1 

-as  if  x a < 0 

i si 


6a 

0 > °an-(i-l) 

1 " b 

n— m 


1 > n-(i-l) 

M,  b 

i n-m 


X1  > V(i-j 

s.  - b 

i n-m 


where 


Au(m)  = sgn(bn_ra){^°  x°  + ^4  + ^TAu} 


0 o 

a if  x a > 0 

1 bi 


-a°  if  x°  a < 0 
1 Si 


\ if  \ a > 0 
if\a<0 


a_  if  Au,  a > 0 
L ) bi  1 

!i'  i „ - 

-ac  if  Au.  a < 0 
i 1 


FI  U 


D 

0 

0 

I 0 


o 

> 6an-(i-l) 

ai 

b 

n-m 

1 

^ 6Vd-i) 

“m. 

b 

i 

n-m 

1 

^ V(i-l) 

as. 

b 

i 

n-m 

iT  i iT  - 

The  terms  \J>  x and  ij>  Au  in  (2.6-27)  are  substituted  for  the  naturally 
bS  1 bs  - 

occurring  terms  x0  and  r Au,  respectively,  since  b„  is  part  of  the 

D o D D 

n-m  n— m ^ 

unknown  plant  parameter  vector  (2.6-8).  As  long  as  oC  , i = 1,  . . . , m, 

Si 

satisfy  the  bounds  indicated  above,  it  may  be  shown  through  Liapunov  stability 

analysis  that  the  stability  property  of  the  MRAS  is  retained. 

The  stable  discontinuous  signal  synthesis  MRAS  design  introduces 

constraints  on  the  plant  zeroes  and  the  discontinuty  surface  coefficient  cr 

which  do  not  appear  either  in  the  discontinuous  parameter  adaptation  design 

or  in  any  of  the  continuous  adaptation  designs . These  constraints  are  intro- 

T 1 T — 

duced  by  the  terms  bg  xg  and  bg  Au  in  (2.6-25),  and  become  significant  to 
system  stability  only  when  the  MRAS  exhibits  a sliding  motion.  We  analyze  the 
augmented  signal  synthesis  case  here;  the  direct  case  exhibits  similar  prop- 
erties. This  aspect  of  MRAS  design  was  first  considered  in  [64]  for  the  case 

when  c =1. 
n 

During  sliding,  the  condition  a *=  a = 0 and  (2. 6- 25b)  provide  us 

with  an  equivalent  augmented  control  [60],  denoted  Au^  : 

T 

. (m)  1 V . . T o . ..T  1 .T 

Au  = r l + oa  x0  + 6b  x,.  - b„  Au  j 

eq  d c a Mb 


(2.6-28) 
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This  control  corresponds  to  the  effective  continuous  value  of  Au^  ’ necessary 
to  maintain  motion  of  e(t)  along  0=0.  The  motion  of  Au,  defined  by  (2.6-26), 
during  sliding  is  determined  by: 


0 1 ...  0 


Au  = 0 • ' • 1 Au  + 


-Cn(^ 

n-m 


Zv  + 4 ♦ 6b^} 

n 


(2.6-29) 


Considering  Au^  as  the  output  of  this  subsystem,  the  transfer  function  from 

Au^  to  Au  is: 

eq  1 , 


H(S)  - 


(i  )s"  + Ja-tL  S--1  + ...  + A. 

n n-m  n-m 


(2.6-30) 


i.e.,  the  subsystem  poles  of  (2.6-30)  are  shifted  from  the  plant  zeroes  of 


(2.6-2)  by  an  amount  dependent  upon  the  value  of  (c^-1) . Thus  we  are  re- 
quired to  place  the  following  restrictions  on  the  plant  zeroes  occurring  in 
(2.6-2)  and  on  cn  in  order  to  guarantee  asymptotic  stability  of  the  MRAS 
during  sliding: 


a)  The  plant  numerator  coefficient  vector  b<,  must  be  such  that 
a positive  scalar  cr  exists  so  the  poles  of  (2.6-30)  are 
asymptotically  stable.  Necessary  conditions  for  the  existence 


of  c are  [61] : 
n 


.c 


i)  No  elements  of  bg  may  be  zero, 

T 

ii)  All  elements  of  bc  = [b  ...  b _ ] must  have  same 

o n n-m+1 

sign. 

b)  given  that  b<,  satisfies  (a),  cr  must  be  chosen  to  place  the 
poles  of  (2.6-30)  in  the  strict  left-hand  s-plane. 

To  illustrate  these  constraints,  we  consider  three  examples  for  the  case  m = 
3,  examining  the  effect  cr  may  have  on  various  systems  of  the  form  (2.6-30). 


Example  2.6.7 


a)  H(S) 


(— )S3  + 3S2  + 4S  + 2 
c 

n 


This  case  corresponds  (for  cr  = 1)  to  plant  zeroes  in  the  left-half  s-plane 

at  s = -1,  s = - 1 ± j.  Figure  2.6-2(a)  illustrates  that  for  cn  > 0.167,  the 

poles  of  H9s)  are  asymptotically  stable,  but  when  cn  - 0.167,  the  complex 

poles  move  into  the  right-hand  s-plane.  In  the  limit  as  c a 2nd  order 

system  is  obtained  with  poles  at  s = - .67  ± j.47.  Thus  constraint  (a)  is 

satisfied,  and  it  is  possible  to  design  a stable  discontinuous  MRAS  by  choosing 

c > 0.167. 
n 


b)  H(S) 


(~)S3  + S2  + 1.25S  + 6.5 
n 


This  case  corresponds  to  plant  zeroes  at  S = -2,  S = + .5  ± j/3,  yielding  an 

unstable  subsystem  H(S)  when  c **  1.  However,  if  c >5,  the  unstable  poles 

n n 

shift  into  the  asymptotically  stable  region  (Figure  2.6-2(b)).  In  the  limit 
as  cn  "*■  °°»  a 2nd  order  system  is  obtained  with  poles  at  S ■ - .625  ± j 2.48. 


(c)  Storting  ots  = -l,s  = l±j 


M-MSl 


Figure  2.6-2  cn- locus  of  Roots  for  3^  Order  Systems 
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Again,  constraint  (a)  is  satisfied,  despite  the  fact  that  the  original  plant 
is  non-minimum-phase. 


c)  H(S)  * ± 5 

(— )S'J  - S + 2 
c 


This  case  corresponds  to  plant  zeroes  at  s ■ -1,  s ■ + 1 ± j.  Since  it 
violates  both  necessary  conditions  in  constraint  (a),  the  unstable  poles 
remain  unstable  for  all  values  of  cr,  as  seen  in  Figure  2.6-2(c).  In  the 
limit  as  cn  °°,  an  unstable  2nd  order  system  is  obtained  with  poles  at  S = 

± 7i. 

m 

Relation  between  Disturbance  Rejection  in  VSS  and  the 
Discontinuous  MRAS 

Variable  structure  system  (VSS)  design  theory  [60]  has  been  developed 
to  solve  the  problem  of  rejecting  both  measurable  and  unmeasurable  input 
disturbances.  Consider  the  stable  single-input  phase-canonic  plant: 


x = Ax  + 


(f-w) 


(2.6-31) 


where  w is  the  control  and  f is  a measurable  disturbance.  A variable  struc- 
ture control  design  which  rejects  the  disturbance  f,  thus  allowing  x to 
approach  the  origin  asymptotically,  is: 


& 

* 


PM 


E 


IV 


w = f 


where 


a if  f o > 0 


3 if  f a < 0 


T 

a = c x 


a > 1 


0 < 1 


c chosen  to  provide  stable  sliding  motion. 


c = 1 
n 


When  sliding  occurs,  the  equivalent  control,  weq»  is 

T 

w - il>  f = (1  + - ~)f 
eq  eq  f 


and  the  motion  of  x along  o = 0 is  described  by: 


0 1 ....  0 


0 ...  1 X 


“C1  ‘ ' Cn-1 


where  x - [x,  . . . x , ] 
l n-i 


since  x = - E c.x. . As  x slides  toward  the  origin  along  o = 0,  we  see 
n i-1  1 1 

from  (2.6-33)  that  w approaches  f;  thus  this  design  exactly  counteracts  the 


disturbance  at  the  equilibrium  point  x ■ 0. 


The  existence  condition  for  sliding  is  0 O < 0 
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a * ic  Ax  + 


(2.6-36) 


aa  = 


[•cTAxxT<]  + (l-a)af  if  of  > 0 


[cTAxxTc]  + (l-B)af  if  of  < 0 


(2.6-37) 


and  so  the  existence  condition  for  sliding  is  satisfied  if: 


ct  > 1 + 


B < 1 + 


T 

vAx 


T* 

c Ax 


(2. 6- 38a) 
(2.6-38b) 


For  a particular  choice  of  a > 1,  B < 1,  the  existence  conditions  for  sliding 
will  be  satisfied  when  the  plant  state  and  disturbance  satisfy: 

T 

-d-B)  < (a-l)  (2.6-39) 

Thus  sliding  will  occur  when  ||x||  is  "small  enough"  or  if  |f|  is  "large 
enough";  i.e.,  the  larger  the  disturbance  magnitude,  the  greater  the  region 
along  a = 0 for  which  sliding  will  occur.  Figure  6.2-3  illustrates  this  idea 
for  two  constant  disturbances  of  different  magnitudes,  for  given  a and  B* 

An  important  feature  of  this  design  is  that  sliding  always  exists 
when  x * 0;  that  is,  the  disturbance  f will  not  force  x away  from  the  origin. 
At  the  origin  the  control  exactly  counteracts  the  disturbance  and  (2.6-39)  is 
always  satisfied  for  every  a > 1,  B < 1. 

If  we  examine  the  error  equation  (2.6-25)  for  the  signal  synthesis 
forms  of  the  MRAS,  we  see  that  they  have  the  same  form  as  (2.6-31)  by  making 
the  substitutions: 
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n-m 


b Au 
n-m 


(m) 


= w 


(2.6-40) 


and 


-To.  rT-1  ,T  1 

6a  XS  + Vm  ' bsxs 


6aTXg  + - bgAu 


= + f2+  f3  - f 


(2.6-41) 


We  will  consider  only  the  augmented  signal  synthesis  case  from  this  point. 

Following  the  design  principle  of  (2.6-32),  the  discontinuous  MRAS  control 

b Au^  will  reject  the  disturbance  (2.6-41)  when 
n-m 

b Au(m)  = ijV  + i|>2f.  + i{>3f- 
n-m  T 1 2 3 


where 


♦M 


a if  f±a  > 0 


6 if  fa  < 0 


a1  > 1 
81  < 1 

T 

a * c e 


(2.6-42) 


Since  we  cannot  measure  directly  the  complete  disturbance  signals,  f^,  f^. 

and  f_,  in  the  MRAS,  the  design  for  b Au^  must  be  modified  slightly, 
j n-m 

while  still  retaining  the  desirable  disturbance-rejection  property  of  the  VSS. 


' ... 
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The  following  theorem  suggests  the  solution  to  this  problem  for  the  MRAS. 

Theorem  2.6.2  Given  the  stable  plant  (2.6-31),  but  with  f = yf, 

f measurable  and  y unknown  but  bounded  by  Y . - y < y , then  the  control: 

rain  max 

i 

w = ip  f 

if  of  > 0 
if  of  < 0 

retains  the  disturbance- rejection  property  provided  by  (2.6-32). 


a 


•1, 


a > y 

i 

3 < y 


max 


min 


Proof : The  system  may  be  expressed: 


x = Ax  + 


0 

Y 


(f  - -) 
Y 


Assuming  that  the  control 


w = ^ = 


a > 1 if  af  > 0 


*- 


3 <1  if  af  < 0 


provides  disturbance-rejection,  then  the  control: 


I 
I 
I 
I 
1 
I 
! 
1 
I 
I 
I 
I 
I 
I 
I 
E 

[c 

0 

0 
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w = yipf  = ip  f 


where 


ja  > y if  af  > 0 

i 

8 < y if  af  < 0 

also  provides  disturbance  rejection.  Since  y is  unknown,  if  we  choose: 


a > y 

, J max 

ip  = 


6 

the  conditions  for  disturbance  rejection  are  still  satisfied.  ^ 

Applying  this  theorem  to  the  design  (2.6-42)  for  the  augmented 
signal  synthesis  MRAS,  we  obtain  a result  which  is  identical  to  the  design 
already  obtained  in  (2.6-27b)  from  Liapunov  methods.  Thus  the  discontinuous 
signal  synthesis  control  design  for  MRAS  obtained  from  Liapunov  stability 
analysis  is  equivalent  to  the  solution  of  a VSS  measurable  disturbance  re- 
jection problem. 

The  VSS  design  theory  may  be  extended  in  the  discontinuous  MRAS 

design  problem  to  provide  for  global  existence  of  sliding  for  all  e e Rn. 

This  may  be  accomplished  by  adding  another  term  to  the  control  equation 

T 

(2.6-27b)  of  the  form  [60]  e.  The  elements  of  \p  satisfy: 

/ 


a ^ if  ae ^ > 0 


8^  if  ae^  < 0 
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“t  - “x(ci-i  + V(i-1)  • ci(Vi  + *1)) 


6i  5 "1"(ci-l  + V(i-1>  • ci(Vl  " *1» 


(2.6-43) 


The  Discrete-Time  MRAS  with  Discontinuous  Adaptation 

The  continuous- time  development  for  MRAS  with  discontinuous  adapta- 
tion has  a counterpart  for  discrete-time  systems.  The  concept  of  sliding, 

however,  takes  on  a different  interpretation.  Given  a discontinuity  surface 

T 

in  the  MRAS  state-error  space,  cr(k)  = c e(k),  "sliding"  motion  of  e(k), 
if  it  exists,  does  not  occur  precisely  on  a(k)  = 0,  but  rather  the  motion 
stays  inside  a A-neighborhood  of  a(k)  = 0.  The  size  of  the  A-neighborhood 
will  depend  on  the  magnitude  of  the  discontinuous  gains. 

A discontinuous  parameter  adaptation  algorithm  for  the  Landau 
MRAS  or  ELMA  which  is  a counterpart  to  the  continuous- time  algorithm  of 
Theorem  2.6.1  may  be  shown  to  be  hyperstable,  given  that  a set  of  inequalities 
is  satisfied. 


Theorem  2.6.3  Given  the  Landau  MRAS  of  Figure  2.1-2,  there  exist 
a vector  c and  a matrix  A = diagfo^]  such  that  the  parameter  adaptation 
algorithm: 


p(k+l)  = p(0)  + A sgn  (xg)  sgn  (v) 


where 


v = c e + e 


provides  a hyperstable  MRAS. 


Proof : Transforming  the  MRAS  into  the  Popov  linear — nonlinear 

feedback  structure,  we  obtain  the  system  of  Figure  2.6-4.  The  linear  sub- 
system is  identical  to  that  for  the  continuous  parameter  adaptation  design, 
and  the  vector  c = -a^^  clearly  makes  the  linear  subsystem  strictly  positive- 
real.  We  must  show  that  the  nonlinear  subsystem  satisfies  the  Popov  in- 
equality: 


Z w(k)  v(k)  > - y' 


for  all  K > 0,  y > 0 
o 


The  nonlinear  subsystem  output  w(k)  may  be  expressed: 


w(k)  = Z |x  |[ot.  sgn(v)  - 6p.  (0)  sgn  (x  )] 
i-1  bi  1 1 l 

2N+1  _ _ 

Thus,  w(k)  v(k)  = | v(k) | Z [a.  - 6p  (0)  sgn  (x  ) sgn  (v) ] 

i=l  1 i 

A sufficient  condition  for  w(k)  v(k)  to  be  positive  is: 


* |6P±(0)  | , i = 1 2N  + 1 


Therefore,  the  inequality  * is  satisfied  whenever  A satisfies  **.  g 

Thus  in  order  to  Implement  a discontinuous  adaptation  discrete- 
time MRAS,  we  must  know  bounds  on  elements  of  the  parameter  vector  p^ 
(identification)  or  pg(0)  (control) . This  knowledge  enables  us  to  select 


elements  for  the  gain  matrix  A which  satisfy  the  condition  **  in  the  theorem 


Using  the  notation  in  [60],  this  adaptation  algorithm  may  be  expressed 
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p(k+l)  = p(0)  + 


a.  If  x v > 0 
l s± 


-a  if  x v < 0 
i s± 


ot  * |«Pt(0)| 


(2.6-44) 


Both  direct  and  augmented  signal  synthesis  versions  of  this  MRAS 
algorithm  may  be  derived  in  a way  analogous  to  the  continuous-time  case.  In 
the  direct  case,  u(k)  is  the  S-subsystem  input;  in  the  augmented  case, 
u(k)  *=  z(k)  + Au(k)  is  the  S-subsystem  input,  where  z(k)  is  the  H-subsystem 
input.  The  state-error  equations  are,  respectively: 


e(k+l)  = A^e(k)  + G2[6aTx*  + b^x^  - b^Xg]  - G2  bg  u(k) 

o 

e(k+l)  = ^(k)  + G2[6pTxg  - bg(xg  -:  x^)  ] - G2  bg  Au(k) 

o 


and  the  discontinuous  control  signals  become: 

, T . , T 


♦i" 


% 


M. 


sgn(bs 

V r . OT  O , . T 
){*  xs  + 

U 

a° 

if  x° 

v > 0 

i 

Si 

<-a° 

if  x° 

v < 0 

i 

5i 

6a. 

a°  > 

V 1 

?i 

bs 

, i = : 

o 

( 1 

—1 

\ 

if  Xj^j  v > 0 

K 

if  xj  v < 0 

(2.6-45a) 

(2.6-45b) 
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CHAPTER  3 

THE  MULTI-MODEL  MRAS 

This  chapter  Introduces  the  concept  o£  MRAS  structures  with 
multiple  reference  models,  a class  of  MRAS  which  has  apparently  not  appeared 
In  the  literature  before.  The  Idea  for  such  a MRAS  structure  was  suggested 
by  the  theory  of  dynamic  differential  games  [20] . Differential  games  Is  often 
viewed  as  an  extension  or  generalization  of  optimal  control.  A relationship 
between  the  single-model  MRAS  and  an  optimal  regulator  problem  was  developed 
in  Chapter  1.  This  relationship,  then  serves  as  the  motivation  for  jointly 
developing  a multi-model  MRAS  and  a dynamic  game. 

First  we  introduce  the  two-model  MRAS  (2M-MRAS),  emphasizing  the 
quasi-symmetry  of  this  structure,  and  develop  the  motivation  for  the  2M-MRAS. 
The  concept  of  a (2M-MRAS,  dynamic  game)  pair  is  then  developed.  It  is 
viewed  by  this  author  as  a logical  generalization  of  the  related  (MRAS,  optimal 
regulator)  pair  developed  in  Chapter  1.  Next  we  consider  the  question  of 
stability  for  a discrete-time  2M-MRAS  which  uses  the  Landau  hyperstable 
adaptation  algorithm.  A particular  regulator-type  (2M-MRAS,  dynamic  game) 
pair  is  then  considered.  Comparing  the  equilibrium  solutions  between  the 
2M-MRAS  and  the  dynamic  game  for  an  assumed  player  strategy  establishes  a 
procedure  for  Investigating  possible  equivalence  between  the  two  problems. 

By  equivalence  here  we  mean  that  the  two  equilibrium  solutions  are  the  same. 
Finally  we  analyze  the  effects  of  the  quasi -symmetric  elements  on  the 
equilibrium  for  the  2M-MRAS. 

3.1  The  2M-MRAS 


In  this  section  we  Introduce  the  two-model  MIAS  (2M-MRAS) , examine 
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Its  structure,  and  discuss  the  motivation  for  our  Interest  in  this  class  of 
MRAS. 

Figure  3.1-1  depicts  the  basic  2M-MRAS  parallel  structure.  It  is 
similar  to  the  one-model  MRAS  in  Figure  2.1-1,  with  the  addition  of  a second 
reference  model  (M^)  and  adaptation  algorithm  (A^).  The  S-subsystem  para- 
meter vector  p(k)  is  influenced  by  the  outputs  of  both  A^  and  A£.  A quasi- 
symmetry is  apparent  in  Figure  3.1-1  in  the  sense  that  the  2M-MRAS  is  sym- 
metric at  the  level  of  detail  depicted  in  Figure  3.1-1,  but  is  not  neces- 

4 A. , and  u.  4 u9  in  general.  These 
quasi-symmetric  elements  may  be  used  to  classify  the  2M-MRAS  into  subclasses. 
For  example,  two  subclasses  may  be  obtained  depending  on  whether  u^  = u^  or 
u^  4 U£.  Similarly  classifications  may  be  obtained  using  and  A^,  i =*  1, 

2.  First  we  define  two  partial  ordering  relations  [5],  which  may  then  be 
used  to  separate  the  2M-MRAS  into  various  subclasses. 


sarily  totally  symmetric;  M 4 M2,  A^ 


Definition  3.1,1  If  = u2  “ t*ie  subsystem  states  x^(0) 
(0) , and  the  respective  state  trajectories  satisfy  the  Inequality: 

R K 

E Xw  (k)Xw  (k)  £ E x^  (k)x^  <k)  for  a11  K - 0 
k-0  2 ”2  k-o”l  ^1 


then  we  say  that  satisfy  the  partial  ordering  M^  * M2,  and  that  M^  is 

''globally  uniformly  faster"  than  M2. 


Definition  3.1.2  For  u(k)  e U,  and  for  every  1M-MRAS  compatible 
with  A^  and  Ag,  if  the  1M-MRAS  with  Aq  and  Ag,  respectively,  satisfies: 


{ E eT(k)e(k)  | A } < { E eT(k)e(k)  | Aft}  for  all  K > 0, 

k-0  * k-0  13 


■1 
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where  e(k)  = x^k)  - xg(k),  then  we  say  that  (Aa,Ag)  satisfy  the  partial 

ordering  A - A0 
a U ° 

The  partial  ordering  for  (M]L»M2)  identifies  two  distinct  subclasses 
of  2H-MEAS,  one  subclass  where  the  transient  state  response  of  one  reference 
model  to  the  origin  is  globally  uniformly  faster  than  for  the  other  refer- 
ence model,  the  other  subclass  where  such  a property  does  not  hold.  For 
example,  when  and  are  linear  discrete-time  systems,  implies 

that  the  maximum  norm  for  all  eigenvalues  of  the  system  matrix  A„  is  less 


\ 


than  or  equal  to  the  minimum  for  all  eigenvalues  of  A^  . A norm  for  a com- 
plex eigenvalue,  X = a + jw,  might  be: 


1*1  - 


(3.1-1) 


The  partial  ordering  for  (A^,Ag)  provides  a means  for  classifying  the  2M- 
MRAS  in  terms  of  the  relative  "strengths"  of  the  two  adaptation  algorithms 
A^  and  Ag.  For  example,  for  the  Landau  adaptation  algorithm  (2.1-14),  if  two 
different  adaptation  gain  matrices  satisfy  F^Ck)  > Fg(k),  it  may  be  demon- 
strated through  simulation  that,  when  U = {u  = constant},  A > A.  . In 

a u p 

Section  3.4  we  will  consider  more  specifically  how  perturbing  the  2M-MRAS 
away  from  perfect  symmetry  will  influence  the  equilibrium  trajectories  of  the 
2M-MRAS  parameter  vector.  Both  the  effects  of  (A^A^  and  (M^M^  will  be 
considered. 

We  make  one  further  distinction  here  between  the  1M-MRAS  and  the 


2M-MRAS.  For  the  1-model  case  we  assumed  incomplete  knowledge  about  the 
Plant  (S-aubsystem)  parameters.  For  the  2-mddbl  case,  it  is  assumed  that 
plant  parameters  are  known  by  both  players; .uncertainty  here  for  each 
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player  is  restricted  to  the  effect  of  the  other  player's  reference  model, 
adaptation  algorithm,  and  external  Input  on  the  S-subsystem  state-trajectory. 

For  the  remainder  of  this  chapter  we  will  assume  that  the  M^-  and 
S-subsystems  in  Figure  3.1-1  are  linear  discrete-time,  and  the  adaptation 
subsystems  are  of  the  hyper stable  type  (with  respect  to  the  1M-MRAS) 

introd&ced  by  Landau  and  treated  in  Chapter  2.  Figure  3.1-2  depicts  such 
a system.  We  define  the  following  variables  for  this  2M-MRAS: 


x^  00  * ^x^T(k)  I x^T(k)]  * [yM  (k-n)...yM  (k-1)  j u£  (k-n) . . .^  (k)] 


I - U ! bj<3  - c 


M 


MX 


-V  *"aM 


. . .b„  3 


i,  i Mi  ""Mi 
1 1 n 


(3.1-2) 

(3.1-3) 


xo (k)  - [x°T(k)J  x*T(k)',  x?T(k)l 


■ [ygfr-n). . .yg(k-l)  i u^ (k-n) . . .u^ (k)  J u2 (k-n) . . .u2 (k)]  (3.1-4) 


P*(k)  - [aT(k)|  blT(k)  J b2T(k)] 

- [a  (k) . . .a. (k)  ! b* (k) . . V (k) . . .b2 (k)3 

u l | ii  l n u 

Xg  - [xgT(k)  I x^T(k)] 

^(k)  - [aT(k)|  biT(k)] 

® i (k)  - y^OO  -ys(k) 
e^(k)  - [ei(k-n)...ei(k-l)3 
vi(k)  " *£00 


(3.1-5) 

(3.1-6) 

(3.1-7) 

(3.1-8) 

(3.1-9) 

(3.1-10) 


The  concept  of  the  multi-model  MRAS  has  been  motivated  by  the 
theory  of  differential  dynamic  games.  An  n-player  dynamic  game  may  be 
defined  by  a plant  model  whose  inputs  are  partitioned  into  n distinct  groups, 
associated  with  distinct  players,  and  a performance  index  associated  with 
each  player.  (The  players  are  also  referred  to  as  control  agents  or  decision 
makers  by  some  authors.)  Each  player  desires  to  minimize  his  own  performance 
index  by  implementing  a control  policy.  However,  since  each  player's 
performance  is  affected  by  the  control  decisions  of  the  other  players,  all 
players  cannot  simultaneously  (and  independently)  minimize  their  respective 
performance  indices.  Thus,  unlike  the  (1-player)  optimal  control  problem,  it 
becomes  expedient  to  view  optimality  in  terms  of  a strategy  adopted  by  each 
player,  yielding  strategy-dependent  controls  and  system  trajectories.  The 
most  characteristic  strategies  in  game  theory  include  Nash  [57],  Stackelberg 
[56],  Pareto  [57],  and  minimax  [17]. 

In  Chapter  1 we  developed  a relationship  between  the  1M-MRAS  and 
a linear-quadratic  optimal  regulator  problem.  The  M-subsystem  reflected  the 
desired  performance  of  the  controlled  plant  (S-subsystem) , and  the  optimal 
regulator  had  an  equivalent  interpretation  in  the  optimal  control  context 
for  nominal  plant  parameters.  Thus  we  chose  the  M-subsystem  to  be  the 
optimal  regulator  obtained  using  nominal  plant  parameter  values.  In  the 
2M-MRAS  case,  each  M-subsystem  will  be  chosen  equal  to  the  regulator  obtained 
by  either  optimizing  that  player's  performance  index,  given  some  assumption 
about  the  control  strategies  chosen  by  the  other  player,  or  by  solving  a 
separate  dynamic  game  based  on  assumed  strategy  and  estimated  performance 
index  for  the  other  player.  We  will  now  consider  in  detail  how  the  (2M-MRAS, 
dynamic  game)  pair  is  developed. 
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Consider  the  linear-quadratic  regulator  game  defined  by: 

x(k+l)  = Ax(k)  + BjU^k)  + B2u2(k)  (3.1-11) 

Ji  ■ 2 Jo(xV  +2“l"iix  +2ujMijx  + 2ulGljuJ  ;+“tEiiut  + l “XjuJ>-  <2'1'12) 


(i,j)  = {(1,2), (2,1)}. 


X T I •p 

Considering  u ■ [u^  i u2],  we  rewrite  (3.1-11)  and  (3.1-12)  as: 

x(k+l)  = Ax(k)  + Bu(k) 

Ji  “ 2 k-O^^  +uTMix  +uTRiu)»  1 m l>2 


(3.1-13) 

(3.1-14) 


where : 


B - [B1  ! B2] 


Mi  = 

"Mll 

M12 

M2  = ! 

| 

"21 

>_ 

y 

R1  * 

«n 

T 

G12 

R2  = 

R21 

rT- 

G21 

T* 

GU 

R12 

G21 

R22 

(3.1-15) 


i: 

i: 


Solutions  for  this  dynamic  game  may  be  obtained  by  assuming  that  the  two 
players  jointly  adopt  one  of  the  game  strategies  mentioned  above.  Now 
consider  the  2M-MRAS  with  S-subsystem  transfer  function  matrix  defined  by: 

H (z)  - CcT(zI-A)'1Bi:cT(zI-A)"1B0] 

S 1 • z 

m 

where  C * [l  0 . . . 0] . 

The  M^-  and  ^-subsystems  must  still  be  chosen.  We  consider  several  alternate 
viewpoints  a player  might  take  in  assuming  something  about  the  other  player's 
control  strategy  or  performance  index. 

The  most  optimistic  viewpoint  by  player  (i)  would  be  to  assume 
player  (j)  is  operating  completely  in  the  interest  of  player  (i).  would 


| 

I. 
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then  be  obtained  by  minimizing  with  respect  to  both  u^  and  or 
T T T 

u = [u^  u^].  Using  (3.1-14),  and  assuming: 


Q±  > 0 


R±  > 0 

Qi  " MiRilMi  - 0 


(3.1-16) 


the  optimal  control  solution,  using  the  variational  approach  [18]  is: 


u*(k)  = -R“1{Mi+BTPiH"1Wi}x(k) 


(3.1-17) 


where  P is  the  unique  positive-definite  symmetric  matrix  satisfying: 


pi  ■ «i  - "KS + “IvA 


(3.1-18) 


and  where: 


H1  = I + BR;.1BTPi 


(3.1-19) 


W±  = A - BR~1Mi 


(3.1-20) 


The  closed-loop  optimal  plant  matrix  corresponding  to  (3.1-13),  (3.1-17)  is: 


(3.1-21) 


Thus  the  optimistic  viewpoint  would  result  in  choosing  according  to  (3.1-21). 
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A less  optimistic  viewpoint  which  might  be  taken  is  for  player  (i)  to 
estimate  the  performance  index  for  player  (j)  and  then  to  solve  a dynamic  game 
problem  based  on  this  estimate  and  an  assumed  strategy  for  both  players.  For 
example  an  assumed  Nash  strategy  would  yield  a closed-loop  regulator  as 
defined  by  (3.3-9)  in  Section  3.3.  This  regulator  would  then  be  chosen  by 
player  (i)  as  his  reference  model  M^.  Alternatively,  player  (i)  might  estimate 
directly  the  cost  strategy  to  be  chosen  by  player  (j).  For  example,  suppose 
player  (i)  assumes: 


Uj  - FjX.  (3.1-22) 

Then  player  (i)  would  solve  the  following  optimization  problem: 

x (k+1 ) * AjX(k)  + Biuj[(k)  (3.1-23) 

where:  A j = A + BjFj,  (3.1-24) 

Ji  * ? J0(xV+2uftiX  + UiRiiUi)  (3.1-25) 

where:  ^i  " Qi  + 2FjMij  + FjRijFj  (3.1-26) 

" Mii  + GijFj‘  (3.1-27) 


The  solution  to  this  problem  takes  the  same  form  as  (3 . 1-17)- (3 . 1-21) , with 
appropriate  substitutions.  This  viewpoint  is  similar  to  the  philosophy 
inherent  in  the  MRAS  approach  to  control  system  design.  For  the  1M-MRAS,  the 
uncertainty  in  plant  parameters  is  dealt  with  by  assuming  nominal  values  in 
order  to  obtain  a reference  model  from  the  optimization  process.  Then  the 
adaptation  loop  is  introduced  to  account  for  the  fact  that  actual  plant  para- 
meters may  not  be  at  the  nominal  values  assumed.  In  the  2M-MRAS,  would  be 
obtained  by  assuming  a nominal  Uj  (in  the  regulator  problem  here,  a state- 
feedback  form).  Player  (i)  would  then  assume  that  his  adaptation  algorithm 
compensated  for  the  fact  that  the  actual  Uj  might  differ  from  the  nominal  u^ . 
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As  a special  case  and  most  pessimistic  example  of  this  last  view- 


point, player (i)  would  assume  player(i)  is  operating  to  maximize  player(i)'s 
cost  J^;  i.e.,  a minimax  strategy.  Thus,  player(i)  would  obtain  u^  by 
maximizing  with  respect  to  u^.  The  resulting  solution  for  player  (i) 
would  then  be  used  to  define  the  M^-subsystem. 

In  Section  3.3  we  consider  in  more  detail  the  viewpoint  where  each 
player  estimates  the  other  player's  control  in  order  to  obtain  definitions  for 
and  1^.  We  consider  there  the  case  where  player (i)  assumes  u^*0.  Then 
(3.1-23),  (3.1-25)  become 


x(k+l)  *>  Ax(k)  + Biui 
Ji  “ 2 k-O^V  +2uiMiix  +uiRiiui> 


(3.1-28) 


(3.1-29) 


and  is  obtained  by  minimizing  (3.1-29)  with  respect  to  u^.  This  resulting 
2M-MRAS  is  then  compared  with  the  dynamic  game  where  each  player  adopts  the 
Nash  strategy. 


E 


3.2  Stability  Analysis  for  the  2M-MRAS 

We  shall  consider  the  2-input  1-output  linear  discrete-time 
2M-MRAS  of  Figure  3.1-2  and  (3.1-2)-(3.1-10),  where  A^  and  each  correspond 
to  a hyperstable  Landau  adaptation  algorithm,  studied  in  Chapter  2.  The 
hyperstability  of  the  1-model  case  was  established  by  transforming  the 
MRAS  structure  to  the  Popov  negative  feedback  form  of  Figure  2.1-3  consisting 
of  a linear  and  nonlinear  subsystem.  The  linear  subsystem  was  required  to 
be  positive-real,  and  the  nonlinear  subsystem  to  satisfy  a positivity 
inequality.  Since  the  nonlinear  subsystem  could  be  viewed  as  a linear  time- 
varying  subsystem  with  respect  to  the  input  and  output,  Landau  has  used 
the  time-varying  version  of  the  discrete  positive-real  lemma  [35] 
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to  establish  positivity  of  the  nonlinear  subsystem.  Thus,. our  first 
step  in  analyzing  the  stability  of  this  2M-MRAS  is  to  obtain  an  equivalent 
Popov-type  structure  consisting  of  a linear  subsystem  and  nonlinear  subsys- 
tem interconnected  by  negative  feedback. 

By  analogy  with  the  1-model  MRAS,  we  obtain  a Popov-type  structure 
whose  linear  subsystem  describes  the  dynamics  for  e^k)  and  e2(k),  (3.1-9), 
with  outputs  v^(k)  and  v2(k),  (3.1-10),  and  whose  nonlinear  system  describes 
the  dynamic  equations  for  parameter  errors  6p^(k)  and  6p2(k)  under  the  as- 
sumed adaptation  algorithm.  From  (3.1-8)  we  obtain: 


T 


ei°°  " ^ fi(k)  + 00 

- bj(k+l)^(k) 

(3.2-1) 

where 

fi(k)  " *M  00  " XS(k) 

and 

6p1(k)  ■ PM  - P1(k) 

(3.2-2) 

Thus  the  dynamic  equation  for  the  state  errors  becomes 

ei(k+l)  - I*et(k)  + B2e1(k),  (3.2-3) 

1 2T 
where  I is  defined  in  Appendix  A,  and  B ■ [0  ...  0 1], 

Lemna  3.2.1  The  variable  f±(k)  in  (3.2-2)  is  equal  to  e^k)  for 
all  k,  if  e^O)  - f1(0). 

Proof : By  expressing  the  M^-subsystem  and  S-subsystem  in  the  1-0 


Delay  state-variable  form  (see  (2.2-7)  - (2.2-13)), 


I 


Oc+l) 


\ 


XMi(k)  + 


u±(k) 


Xg(k+l)  = 


As  Ck+1)  Ag (1 ,k+l)  A2 (1 , k+1 ) 


xs(k) 


bJ  (k+1)  (k+1)  rUl 

,1  S2  1 
B2  0 „2 


we  may  obtain  the  dynamic  equation  for  f^k)  as 


fi(k+1)  = \i  Vk)  + ^ -As(k+i)]x°  + [^  (l)~Ag(l,k+l)  ]Xg 

- Ag(l,k+l)Xg  + [B^-Bg^(k+l))Ui  - Bg  (k+l)Uj 

T 

- A^f^k)  + B2(6i^(k+l)xs  -b^(k+l)^} 


T 

+ fi+5pi(k+l)Xg  -b^(k+l)j?} 


iifi(k)  + B2£i(k) 


Thus,  assuming  f^(0)  = e^O),  (3.2-3)  generates  a sequence  for  £^(k)  which 
is  identical  to  the  sequence  generated  for  f^(k) . 


Thus,  if  we  assume  f±(0)  - e^O),  and  denoting 

T 

w±(k)  - ^ (k)6pi(k+l)  - ^ (k)b^(k+l) 


(3.2-4) 


we  may  rewrite  (3.2-1),  (3.2-3)  and  (3.1-10)  as: 


I: 

L 


E±(k)  - a^  e1(k)  + wi(k) 
ei(k+1)  - ^ e±(k)  + B2Wi(k) 


(3.2-5) 


(3.2-6) 


vjL(k)  - w^k) 


(3.2-7) 


The  Landau  adaptation  algorithm  (2.1-14)  is  repeated  here  for 


reference: 


Pl(k+1)  « ?i(k)  + F±xs  (k)v±(k) 


(3-2-8) 


where  p^  is  given  by  (3.1-7)  and  > 0 is  assumed  constant  here,  rather 

than  the  more  general  time-varying  case  (2.1-15).  Then 


6pi(k+l)  - 6p1(k)  - FjXg  (k)v^(k) 


(3.2-9) 


Denote  the  nonlinear  subsystem  outputs  by  z^(k)  - -w^k).  Then: 


zA(k)  - -Xg  6p1(k+l)  + x^  bg (k+1) 


(3.2-10) 


where 


+ -xSl6pi(k)  + xsiFixs1vi(k)  + *sjFjxsjvj(k)  + *s  ^s(k) 


(3.2-11) 


I r 

■:  L 


The  combined  nonlinear  subsystem  outputs  may  be  written: 


' 


zl(i) 


z2(k) 


I -2' 


[0  ;*f]'  ^ 

' SJ\  s2 


<Sp1(W 


<Sp2(k) 


(3.2-12) 


i — T — I —x  ~ — 

*c  ^1  Xc  | X ^2  X. 

S1  1 S1  | s2  2 S2 
+ 


XS  F1  Xs  1 Xs  F2  Xs 

. S1  1 S1  | S2  2 S2 


V0  V, 

v,(k)J 


Combining  (3.2-6) , (3.2-7) , (3.2-9)  and  (3.2-12),  the  Popov  structure  for  the 
2M-MRAS  is  shown  in  Figure  3.2-1;  (the  time  index  1c  has  been  suppressed  for 
convenience).  The  structure  here  is  very  similar  to  the  1-model  case  in 
Figure  2.1-3,  with  the  exception  here  of  an  additional  forcing  term  in  the  non- 
linear sybsystem.  By  choosing  vi  as  in  (3.1-10),  (i.e.,  letting  the  design 
vector  - -aM  ) , the  transfer  function  matrix  for  the  linear  subsystem  is 
simply: 


1 0 

t 

0 1 


(3.2-13) 


which  is  clearly  strictly  positive-real.  In  order  for  the  2M-MRAS  to  be 
strictly  asymptotically  hyperstable,  the  nonlinear  subsystem  is  assumed  to 
satisfy  the  inequality: 


Z zT(k)v(k)  > -y2 
k=0  ° 


for  all  K > 0 
*0  ” 0 


(3.2-14) 


With  respect  to  _ as  state  variable,  v(k)  as  input,  and  z(k)  as  output. 
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the  nonlinear  subaystem  la  actually  linear  tine -varying,  aa  In  the  1-model 
caae,  but  with  the  additional  forcing  term.  If  we  neglect  the  forcing  term 
appearing  here,  we  might  expect,  due  to  Its  similarity  to  the  1-model  case, 
that  this  linear  time-varying  system  is  also  positive-real,  l.e.,  it  satisfies 
the  time-varying  version  of  the  Discrete  Positive  Real  Lemma  (TVDPR).  This 
subsystem  is  realized  by: 


U,B(k),C(k),D(k)} 


' 

To 

\\ ! 0 ' 

r _ i ,t“] 

\ 4 i 

\ 4n+2* 

1 

9 

,T  1 

l 

° ! F2XS2 

i 

[0  x ] 1 x? 

L s 1 s2  J 

-I  _ - i -T  : - 
XS2F2XS2 


_T  ~ _ | _I  _ — 

\Tl\  , 


I 


(3.2-15) 


We  shall  not  prove  that  (3.2-15)  satisfies  the  TVDPR  Lemma;  rather 
we  shall  assume  it  to  be  true  in  order  to  investigate  how  the  additional 
forcing  term  in  the  nonlinear  subsystem  influences  the  positivity  problem. 
Assuming  (3.2-15)  defines  a positive  system,  a sufficient  condition  for  the 
2M-MRAS  to  be  asymptotically  hyperstable  is: 

^(k)] 


[r2(k)ri(k)] 


|_v2  (k) 


> 0 for  all  k > 0 


(3.2-16) 


where 


,T  _ 

ri-*i  V 


This  follows  directly  from  (3.2-14).  A set  of  sufficient  conditions  and  a 
set  of  necessary  and  sufficient  conditions  for  (3.2-16)  to  be  true  are  given 
in  the  following  proposition. 


123 


Proposition  3.2,1 

Assume  that  (3.2-15)  satisfies  the  TVDPR  Lemma;  then  (a)  sufficient 
conditions  for  which  (3.2-16)  is  true  are: 


(r^>r2)  lies  within  the  region  defined  by  the  linear  constraints: 


(al)  rx  > 0,  r2  > 


(a2)  x1  £ 0,  r2  £ 


1 + *i 


I<*2+y1  (ot2— °1)  J 


(a3)  r2  > 0,  rx  > x + y 


2 [a1+Y2  (c^-ap  ] 


(a4)  r2  £ 0,  rx  < 


Y2  -2 


vheT.  Yl  - i 0 

°i  * [0  "4  11 


Spi 

SpJJ 


and  (b)  a necessary  and  sufficient  condition  for  which  (3.2-16)  is  true  is: 


^rl*r2^  within  the  elliptic  region  defined  by: 


(1+Y1)r^  + (1+Y2)r2  - (Y1+Y2>rir2 


- [a2+Y1(a2-a1)]r1  - [^-^(aj-o^)  ]r2  < 0 

Proof 

a)  A sufficient  condition  for  (3.2-16)  to  be  true  is  that  the 
separate  terms  r^  > 0 and  r^v  £ 0 for  all  k.  From  (3.2-4)  and  (3.2-7),  we 
obtain 
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I 

I 

I 

I 

1 

0 

Id 

G 


* « 
V1 

(1+Y2)  - Y2  a x - r2 

1 

" A 

v2 

"Yi  (1+Yi)  a2  - 

where  A»1  + y1  + Y2>0 


Thus 


r2*62rl~(1+Y2^r2  + l<1+Y2)al~62a2^ 

r2vl  " l - 0 

^{-(l+Yj)^  + BjT2  + [(1+Y1)a2-B1a1]} 

riV1  * A — 0 


d'2 

Since  A > 0,  these  two  equations  yield  the  sufficient  conditions  (al)  - (a4) 
directly. 

b)  By  considering  (3.2-16)  directly,  r^v2  + r2Vi  - ® becomes, 
using  the  definition  above  for  and  v2: 

+ Yirjr2  + [(l*Y1)tt2-Yiailri  ~ (l-hr2)r2  + Y2rlr2  + 1 (1+Y2)a1-Y2a2i * 


^ 0 


yielding  the  ellipse  constraint  directly.  Since  (b)  is  necessary  and  suffi- 
cient while  (a)  is  only  sufficient,  the  constraint  region  (a)  is  completely 
contained  within  the  constraint  region  (b)  in  the  (r^,r2) -plane. 

To  illustrate  the  constraint  conditions  of  this  proposition,  con- 
sider the  following  example  which  is  valid  for  a single  time  instant  k . 

o 

Example  3.2.1 

Let  Yj^  = Y2  = 1»  = a2  = 1 at  k = kQ.  Then  the  linear  constraints 


al)  rx  > 0,  r2  > 2^-1 


D 


are: 
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a2)  < 0,  r2  < 2^-1 

a3)  r2  > 0 rx  > 2r2~l 

a4)  r2  5 0 ^ < 2r2~l 

and  the  elliptic  constraint  Is: 

b)  2r^  + 2r2  - 2^^  - ^ - r2  < 0 

The  regions  satisfying  the  constraints  are  shown  in  Figure  3.2-2.  The  in- 
terior double-hatched  region  corresponds  to  the  sufficient  conditions  (a) 
while  the  circumscribed  ellipse  corresponds  to  the  necessary  and  sufficient 
condition  (b).  ^ 

Although  this  proposition  suggests  the  possibility  that  the  2M- 

MRAS  may  be  made  hyperstable  by  selecting  the  external  input  sequences, 

(ujOO.u^k)},  appropriately,  and  thus  {x^  (k) , x|  (k)  } , a straightforward  method 

of  doing  so  is  not  possible  because  the  region  constraining  (r^,r2)  depends 

on  the  unknown  parameter  error  vectors,  6p’1(k)  and  6p2(k),  imbedded  in  cx^k) 

and  a2(k).  Thus,  the  most  conservative  strategy  is  to  choose  Xg  and  Xg  to 

maintain  the  values  of  r^  and  r2  as  close  to  the  origin,  as  possible,  since 

(rl»r2)  ■ (0,0)  always  lies  in  the  constraint  region,  independent  of  and  a.^. 

This  strategy  is  possible  since  the  model  coefficients,  bu  and  b are  assumed 

**1  **2 

known  by  the  respective  players.  However,  we  should  remember  that  we  are 

working  under  the  (as  yet  unproved)  assumption  that  the  linear  time-varying 

r2 

system  is  positive  when  is  ignored. 

rl 


WWWUI 
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3.3  A Particular  (2M-MRAS.  Dynamic  Game)  Pair 

In  this  section  we  will  formulate  a discrete-time  dynamic  game 
with  quadratic  performance  indices  and  a 2M-MRAS  according  to  the  methods 
described  in  Section  3.1.  Adopting  a Nash  strategy,  the  closed-loop  Nash 
solution  for  the  dynamic  game  will  be  stated,  and  then  some  Nash  solution 
examples  will  be  compared  with  the  equilibrium  parameter  values  of  the 
corresponding  2M-MRAS  structure.  For  this  comparison  we  shall  adopt  the 
viewpoint,  described  in  Section  3.1,  that  player (i)  assumes  a nominal 
function  for  player(j)'s  control.  Thus  the  reference  models  and  M2 
are  obtained  by  minimizing  the  performance  indices  and  J ^ with  respect 
to  u^  alone  and  alone,  respectively.  To  simplify  the  minimization  process 
we  take  player (j)'s  assumed  control  to  be  Uj*0.  We  may  argue  that  this 
choice  is  no  less  reasonably  than  another  if  player (i)  has  no  apriori 
knowledge  about  u^ . 

The  two-player  dynamic  game  to  be  considered  is  as  follows.  The 
plant  is  defined  by: 


x(k+l) 

y(k) 

where  A has  the  observable 


= Ax  (k)  + B^  + B2u2 
* CTx(k)  + b*ux  + b^u2 
canonic  form, 


(3.3-1) 

(3.3-2) 


■ 
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. „T  r,  i 
and  B.^  ■ [b^  . 


. b*],  CT  - [1  0 . . . 0], 


The  performance  Indices  for  the  two  players  are  given  by: 

00 

- j ^(Q1y2(k)+Riiu^(k)+Rijuj2(k)),  i = 1,2 

j = 1.2 


(3.3-3) 


Using  (3.2-2),  we  obtain 


Ji  = 2 ^o(^TQi^2u1MiixfRiiu2+2uiGijuj+2ujMijXfRiju2) 


where  = Q^CCT 

Mii  - boVT 

Rii  ■ Ra + 

Rij  - *«  + (h°)^ 

Gij  ' b>o\ 


> 


J 


(3.3-4) 


(3.3-5) 


The  plant  parameters  a = [an  . . . a^],  b*  = [>n  . . . b^]  in  this  problem 
correspond  to  the  equivalent  parameters  in  the  2M-MRAS  S-subsystem.  The  closed- 
loop  Nash  solution  for  this  dynamic  game,  using  the  variational  approach,  is 
[11]: 


-R 


-1 


%<k>  • -fIMrVjiVx<k) 


(3.3-6) 


M. 


R11  R22 


R11R22  ' G12G21 
Mii  + b^k1h"1w 


where 
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H = 1 + J |TlRllBlKl  + T2R22B2K2J 
" ‘ 4 - A ITAkl  + T2E2^22J 


ti-  Bi  - YjjSi 


(3.3-7) 


J 


and  where  is  a matrix  which  satisfies  the  matrix  equation: 

h ■ Q1  - 1 {MiiRu[Mi-Gij8jjlMJJ  + 1SjRjjl‘GjiRiiMi+V) 


- i {GnR^i+Hi)lRji1("ir  i ‘Vl^-ViX+V  (3-3-8> 


+ GijRii(Mi"GijRjjMj)  + + ATKiH"1W 

Sufficient  conditions  on  A,  Qi£,  R^,  R^,  Mli,  , and  G for  existence 

of  solutions  to  (3.3—8)  are  not  known.  From  (3.3-6)  the  closed— loop  system 
equation  becomes: 

xN(k+l)  = rtx^k)  (3.3-9) 

Unlike  the  optimal  control  problem,  the  system  matrix  H_1W  is  not  guaranteed 
to  be  asymptotically  stable,  so  the  solution  is  practical  only  in  the  case 
when  this  variational  approach  yields  a stable  system. 

In  order  to  specify  reference  models  and  1^,  we  shall  adopt 
the  viewpoint,  described  earlier  in  this  section,  that  player (i)  assumes 
a nominal  control  u^  = 0,  and  then  minimizes  with  respect  to  u^  alone. 


4 
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With  these  assumptions  the  optimization  problem  becomes 


x(k+l)  = Ax(k)  + B^u^ 


(3.3-10) 


OO 

* y £ (x^x  + 2u^M^x  + R^u^) 


(3.3-11) 


where  A,  B^,  Q^,  and  R are  as  defined  for  the  original  dynamic  game. 
The  solution  to  this  optimal  control  problem,  using  the  variational  approach 


[18]  is: 


ui  ‘ -Ra[Mii+BipA1V* 


(3.3-12) 


where  = I + BjR^B^P 


(3.3-13) 


W±  = A - BiR"1M1 


(3.3-14) 


and  where  P is  the  unique  positive-definite  symmetric  matrix  which  satisfies 


the  matrix  equation: 


p = Q - MTr"1M  + W^P  H~^W 
i vi  iKi  i wiriHi  wi 


(3.3-15) 


The  closed  loop  system  equation  becomes,  from  (3.3-10)  and  (3.3-12): 


x (k+1)  « H~^WjX*(k) 


(3.3-16) 


and  is  guaranteed  to  yield  an  asymptotically  stable  system.  These 

optimal  systems  (3.3-16)  define  the  reference  models  and  for  the  2M- 
MRAS.  We  define  the  following  2M-MRAS  restrictions  so  that  the  2M-MRAS  and 


the  dynamic  game  have  the  same  regulator-type  structure.  First  of  all,  the 


parameter  adaptation  of  the  S-subsystem  must  be  restricted  to  the  a(k) 


portion  of  the  parameter  vector  p(k)  defined  in  (3.1-5).  If  adaptation  of 


> v»  J 
;•  ^ 
1* 


the  b (k)  portions  of  p(k)  is  allowed,  change  to  the  performance  index 
matrices  M^,  M^,  R^,  R^,  and  is  implied  due  to  changes  in  b*. 

Further,  the  convergence  of  bi(k)  to  some  fixed  point  for  the  2M-MRAS  has  no 
counterpart  in  the  dynamic  game  plant  equation.  In  the  game  context,  the 
parameters  (B^.b*)  are  assumed  to  be  fixed  known  quantities.  These  same 
fixed  values  will  therefore  be  assumed  for  the  M^- , M^-,  and  S-subsystems. 

We  will  refer  to  this  as  the  "pole  adaptation  only"  restriction  for  the 
2M-MRAS.  The  second  restriction  requires  that  the  external  inputs  and  u^ 
of  the  2M-MRAS  equal  zero  throughout  the  MRAS  operation.  This  corresponds  to 
the  fact  that  the  dynamic  game  under  consideration  is  assumed  to  be  of  the 
regulator-type . 

We  now  present  several  examples  where  we  compare  the  Nash  solutions 
of  the  dynamic  game  with  the  2M-MRAS  equilibrium  solution.  The  comparisons 
are  done  in  terms  of  the  Nash  equilibrium  plant  poles  and  the  2M-MRAS 
S-subsystem  equilibrium  poles.  This  comparison  is  followed  by  an  exact 
analytic  solution  for  the  S-subsystem  equilibrium,  in  Proposition  3.3.1.  The 
optimal  solutions  to  define  and  and  Nash  solutions  are  found  numerically 
by  iteratively  solving  (3.3-15)  and  (3.3-8)  until  a fixed  point  is  reached. 

Example  3,3.1 

s t 

Consider  the  1 order  dynamic  system  whose  transfer  function  and 

state  variable  representation  are  given  by: 


x (k+1 ) - -2x(k) +U]L(k) +u2(k) 
y (k)  = x(k) 

The  performance  indices  for  each  player  are  given  by: 


MM!***1*"!1' 


1 

I 

I 

I 

I 

I 

T 

l 
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Oft 

J1  " 2 k2o(x2(k)  + ul<k» 

J2  " \ Jo(’001  x2(k)  + u2(k))- 

In  terms  of  plant  poles,  the  two  optimal  solution  (and  M^- subsystem)  poles 
are  found  to  be: 

Mx  = -.382 
= -.5 

The  closed-loop  Nash  solution  yields: 


c 


H(z) 


Cz-.99)2  j (z-.99)2 
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’ 1.98  l" 

0 

-f 

0 

x(k+l)  = 

x(k)  + 

^(k)  + 

-0.9801  0 

1 

1 

u2(k) 


y(k)  = [l,0]x(k) 


The  performance  Indices  for  each  player  are  given  by: 

1 o' 

0 0. 

1.001  0 

0 0 


- J Z (xT 
k=0 


1 00  x 

J2~  2 1 


, 2. 

X + 


k=0 


x + u2) 


In  terms  of  plant  poles,  the  two  optimal  solution  (and  M^subsystem)  poles 
are  found  to  be: 

= 0.3728  ± j 0.3004 
^ = 0.8755  ± j 0.1096 
which  corresponds  to: 


I 

[ 


T 


[-0.21,0.79] 


- [-0.81,  1.74] 

The  closed-loop  Nash  solution  yields  poles  at: 
Nash  =■  0.3728  ± j .3003 


«rt*- 
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The  2M-MRAS  was  set  up  for  pole  adaptation  only,  - diag(100,100) , 

and  the  1-0  Delay  state  initial  conditions: 

=*  x^(0)  = [0,1,00] 

Xg(0)  = [0,1,0,0,0,0] 

The  equilibrium  poles  for  the  S-subsystem  depended  on  the  initial  pole  values. 
However,  in  all  cases,  the  equilibrium  poles  tended  toward  a compromise  be- 
tween the  poles  for  the  two  reference  models.  For  example,  when  the  initial 
S-subsystem  poles  lie  at  (.5,. 5),  they  move  to  final  values  (0.674  ± j 0.321); 
initial  values  at  (.9, .9)  converge  to  (0.653  ± j 0.325).  This  dependence  on 
the  initial  values  is  believed  to  be  due  to  the  fact  that  the  regulator- type 
2M-MRAS  allows  only  a short  transient  period  during  which  adaptation  may  take 
place  before  the  subsystem  states  and  state-errors  become  too  small  to  affect 
further  parameter  change.  The  stationary  value  a in  each  trial  represents 

D 

the  extent  to  which  as(k)  approached  the  true  system  equilibrium  before  adap- 
tation was  halted. 

The  S-subsystem  equilibrium  is  in  contrast  again  to  the  Nash  solu- 
tion which  falls  very  close  to  the  faster  of  the  two  optimal  solutions.  ^ 

The  properties  of  the  closed-loop  Nash  strategy  equilibrium  and 

the  2M-MRAS  equilibrium  in  the  above  examples  suggest  a basic  difference 

between  the  two  problems.  First  we  discuss  the  Nash  solution  property. 

Denote  = Qj/Rii  38  a mea8ure  °f  the  relative  weight  attached  to  the  state 

for  each  player  in  the  dynamic  game  as  compared  with  a normalized  weight  of 

2 

1 applied  to  the  control  energy  u^k).  Then  player  1 (since  > e2)  pro- 
vides a greater  proportion  of  the  regulation  control  in  the  sense  that  his 
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Nash  feedback  gain  is  larger  than  that  of  player  2.  Further,  player  2 pro- 
vides less  control  in  the  Nash  equilibrium  case  than  he  would  in  the  simple 

optimal  control  case,  yielding  an  improved  cost  J < J . This  is  because 

2N  2 

player  1 is  providing  "free"  control  from  the  viewpoint  of  player  2,  assisting 
player  2 to  achieve  his  goal.  In  Example  3.3-1,  the  solution  of  (3.3-19)  is 
= 4.2361  while  the  solution  of  (3.3-7)  is  = 4.2277,  implying  almost 
identical  control  efforts  for  player  1.  In  contrast  = 2.99  while  = 
0.004,  a significant  reduction  in  control  effort  for  player  2 in  the  Nash 
game  as  compared  with  the  optimal  control  solution.  In  the  limit,  as  e^/e^ 
approaches  °°,  the  Nash  transient  response  approaches  the  optimal  transient 


response  of  player  1.  In  our  examples  e^/e2  = 1000,  which  is  effectively 
00  as  far  as  this  phenomenon  is  concerned. 


The  2M-MRAS  equilibrium,  on  the  other  hand,  appears  from  the  ex- 
amples to  represent  a true  compromise  between  the  two  optimal  trajectories  as 
specified  by  the  reference  models  M^  and  M^.  The  "cooperation"  which  appears 
in  the  Nash  regulator  solution  does  not  occur  here.  The  goal  in  this  case, 
rather  than  taking  the  state  to  the  origin  as  quickly  as  possible  compromised 
by  the  necessary  control  effort,  is  to  take  the  S-subsystem  state  to  the  ori- 
gin along  the  exact  trajectory  specified  by  either  or  M2.  Thus,  it  is  not 
surprising  to  find  that  the  S-subsystem  poles  in  our  two  examples  reach  an 
equilibrium  lying  in  a region  between  the  poles  corresponding  to  and  M2. 

A more  precise  and  quantitative  description;  of  the  2M-MRAS  equilibrium  is 
given  by  the  proposition  below.  First  we  state  several  lemmas  which  will  be 
used  to  prove  the  proposition. 
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E 


Lemma  3.3.1  Given  the  set  of  linearly  Independent  vectors  in  Rn, 

{xq,  . . set  diadic  linear  transformations 

T 

^Xi  * xixi  I = 0»  • • .»  n - 1}  has  the  following  properties: 

1)  (R  , . . . , R , } spans  Rn 

o n-1 

n-1 

2)  n N - {0} 

i-0 

where  ^ denotes  the  range  space  of  X±  and  N±  denotes  the  null-space  of  X^. 
Proof: 

1)  The  range  space  R of  a diadic  linear  transformation  X = xxT 

has  dimension  * 1,  and  R = {y  e Rn  | y =,  ax, a e R}.  Given  that  {xq,  . . 

xn-1}  is  linearly  independent,  and  e R±,  where  a±  ? 0 but  otherwise 

arbitrary,  then  {aQXo,  . . .,  an_^xn_l^  is  also  a linearly  independent  set. 

Therefore  {R  , . . . , R ,}  spans  Rn. 
o n-i 

T 

2)  The  null-space  N of  a diadic  transformation  X - xx  has  dimen- 

n T 

sion  n-1,  and  is  equal  to  the  hyperplane  N ■ (y  £ R I x y = 0>.  If  y e W , 

o 

. . . , W . , then 
n-1 


y - 0 


But  {xo,  . . .,  xn_^}  is  a linearly  Independent  set,  so  y - 0 is  th-  only 

vector  which  may  satisfy  this  equality,  and  thus  lie  simultaneously  in  each  M 

n-1  1 

that  is,  the  intersection  fl  N 

i-0  1 ■ 
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The  next  two  lemmas  deal  with  the  linear  Independence  properties 
of  the  state  sequence  {xq,  x^,  . . x^}  generated  by  the  phase  canonic 

matrix  A.  We  separate  Rn  Into  two  regions,  denoted  G+ (A),  G (A),  which  are 
defined  as  follows: 

G+(A)  = {x  e Rn  | (A,x)  is  a completely  controllable  pair} 


G~(A)  - Rn  - G+(A) 


Lemma  3.3,2  The  two  regions  G (A)  and  G (A)  satisfy  the  following 


conditions: 


a)  G (A)  is  A- invariant;  i.e.,  AG~ (A)  c G-(A) 

b)  If  A is  nonsingular,  G+(A)  is  A- invariant 

c)  If  A is  singular,  AG+(A) C G~(A) 


Proof : 

a)  x e G (A)  implies  R * [x  ,Ax  , . . . , An  ^"x  ] is  singular. 

o o o o o 

ARq  = [Axo,A2xo,  ...»  AnxQ]  = [XjAxp  . . . , An  ^x^]  * R^  is 
also  singular.  Therefore  x^  e G~(A). 

b)  xq  e G+(A)  implies  Rq  is  nonsingular; 

■ ARq  is  also  nonsingular;  therefore  x^  e G+(A) 

c)  xq  e G+ (A)  implies  Rq  is  nonsingular,  but  ARq  = R^  is  singu- 
lar. Therefore  x^  e G (A). 

For  the  cyclic  matrices  under  consideration  (cyclic  since  A is  in 
phase-canonic  form),  G+(A)  + <p  (the  empty  set).  We  may  characterize  G~(A) 
by  applying  results  due  to  Kalman  and  Heymann  [23,14]  to  the  Jordan  Canonic 
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Form  of  A (which  has  only  one  Jordan  block  associated  with  each  eigenvalue) . 


The  results  are  rephrased  so  as  to  be  consistent  with  the  development  here. 


Lemma  3.3.3  [23,14]  Let  Aj  be  a cyclic  matrix  in  Jordan  Canonic 


Form,  with  Aj  = diag  [J^,  . . .,  J^],  3^  the  Jordan  block  associated  with 
eigenvalue  X . Partition  x e Rn  to  correspond  with  the  Jordan  blocks  of  A ; 

1 J 

T T T n 

x ■ [x^  . . . Xg]-  Then  x is  a generator  of  R (or  equivalently,  (AJtx) 


is  a completely  controllable  pair)  if  and  only  if  each  x £ 0,  where  x 

H 


is  the  last  element  of  x^,  i = 1,  . . .,  S. 


From  this  lemma  we  may  characterize  G (AT).  In  the  coordinate 

J 


system  for  the  Jordan  Form,  x e G (A^)  must  lie  in  a hyperplane  defined  as 


follows: 


Ht  = (x  | y^x  = 0;  y^  - [0  . . . 0 a 0 . . . 0],  a 0} 


(3.3-18) 


where  a appears  in  the  i. — position  in  y. 

JO 


There  are  Shsuch  hyperplanes  making  up  G~(AJ).  Therefore 


G"(A  ) = U H, 
i-1  1 


(3.3-19) 


With  respect  to  the  coordinate  system  associated  with  the  phase-canonic  re- 


presentation, A,  the  hyperplanes  are  simply  rotated  by  a linear  transforma- 


tion to  a different  orientation  in  R . 


From  (3.3-18)  and  (3.3-19),  "almost  every"  x e Rn  lies  in  G+(A  ) , 

J 

in  the  sense  that,  if  x e G^Aj),  i.e.,  x £ H for  some  i - 1 S,  then 


an  arbitrarily  small  perturbation  of  x in  any  direction  other  than  in  the 


hyperplane  ^ to  x ■ x + Ad,  A > 0 d t H^,  places  x e G (Aj). 


: r*  * 

f ' 

: - 
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Preposition  3.3,1  If  for  the  2M-MRAS  the  following  conditions  hold: 

a)  U1  * u2  " °» 

b)  pole  adaptation  only, 

c)  ^ * F2  = F * FT  > 0, 

d)  ag(k)  converges  at  k = kQ  to  ag  which  corresponds  to  stable 
S-subsystem  poles, 

e)  x°(kQ)  e G+ (A*) 


then  a. 


\ + % 


Proof : The  adaptation  equation  is  given  by: 

ag(k+l)  = ag(k)  + Fx°[v1(k)  + v2(k)] 

From  Proposition  3.2.1,  we  obtain  from  the  regulator  cs 


(1+Y)  -y  6a 


1 

1 + 2Y 


T o 

1 XS 


-Y  (1+Y)  6a i 


T 

where  y = x°  Fx°,  6ai  = a^-a^k) 
Substituting  into  the  adaptation  equation. 


o o 
2Fx  x 

a (k+1)  = [1°  - —4. 


n w vj 

Fxx„ 
S S 


- n - + rrtr 

When  the  stable  equilibrium  point  is  reached  at  stage  k , 


T 

FXSXS  (2aS  ■ (aML+aM2))  = 0 for  all  k > ko. 
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: 


[ 


f 


I 


i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

r 

l. 

0 

[ 

L 


Since  F 


X 

> 0,  2a^  - (a^  +a^  ) lies  in  the  (n-1)  dimensional  null-space 


of  X^  = x“(k)Xg  (k)  for  all  k ikQ.  Assumption  (e)  implies  (A^x^Ck^)  is  a 

completely  controllable  pair,  and  therefore  R = [x°(k  )A*x°(k  ) . . . 

o S o S S o 

A*11  ^x°(k  )]  = [x°(k  )x°(k  +1)  . . . x°(k  +n-l)]  is  nonsingular.  Based  on 

ddO  & O o O DO 

the  characterization  of  G (A)  in  (3.3-25),  we  claim  that  assumption  (e)  is 
not  significantly  restrictive. 

T 

Denoting  X = [x  (k  +l)x  (k  +1)],  and  using  Lemma  3.3.1, 
jL  b o bo  i 

^ n-i 

[2a  -(a,.  +a  )]  must  lie  in  the  intersection  (1  W. . But  since  (x_(k  ), 

S M2  1=q  1 So 

. . . x°(k  + n-1)}  is  a linearly  independent  set,  the  intersection  contains 

b O 

only  the  origin. 

* , * \ \ 

Therefore  2ag  - (aM  +aM  ) = 0,  or  ag  = ^ . 

This  proposition  implies  that  the  equilibrium  S-subsystem  poles 
may  be  obtained  from  the  characteristic  equation: 


n 1 v / 

z - T Z ( 


‘iVW2 


n-i 


= 0 


(3.3-20) 


3.4  Quasi-Symmetric  Effect  on  2M-MRAS  Performance 

The  2M-MRAS  quasi-symmetric  elements,  namely  the  adaptation  sub- 
systems and  reference  model  subsystems,  may  exhibit  the  partial  ordering 
properties  introduced  in  Section  3.1.  For  the  2M-MRAS  using  Landau  adapta- 


tion, the  notion  of  A^  > A^  corresponds  to  the  condition  where  the  adapta- 
tion gain  matrices  satisfy  F^  > F^.  In  the  1M-MRAS  context,  we  find  that 
the  convergence  rate  of  the  state-error  to  the  origin  is  greater  for  F^ 
compared  with  F^. 

The  partial  ordering  > M2  indicates  that  every  subsystem 


"I 


eigenvalue  has  a smaller  magnitude  than  the  M2  eigenvalue  of  least  magnitude. 
Thus,  the  transient  motion  of  x^k)  to  its  steady-state  value  is  more  rapid 
than  the  corresponding  motion  of  xM  (k) , given  that  the  initial  conditions 


of  both  subsystems  are  equal. 


In  this  section  we  analyze  these  partial  ordering  properties  of 


the  2M-MRAS  and  their  influence  on  its  performance.  We  restrict  our  analysis 


to  the  regulator-type  2M-MRAS  (u  = u2  = 0),  with  pole  adaptation  only,  and 


we  make  assumptions  about  the  form  of  and  F2  as  necessary  to  obtain  re- 


sults which  most  clearly  illustrate  the  main  ideas. 


During  operation  of  the  regulator-type  2M-MRAS,  we  shall  also  as- 


sume that  initial  states  (0)  = (0)  = x°(0),  and  that  x^  , x^  , and  x° 


all  converge  to  the  origin.  It  is  not  necessarily  true  that  x converges  to 

u 


the  origin;  the  2M-MRAS  may  become  unstable  for  certain  combinations  of 


(F1,F2)  and  initial  states,  in  which  case||  x°  ||  becomes  unbounded.  Since 


only  pole  adaptation  is  assumed  here,  we  are  interested  in  the  motion  of  the 


S-subsystem  parameter  vector  ag(k).  In  subsequent  equations,  we  will  delete 


the  time  index  when  no  confusion  is  possible. 


Adaptation  of  the  parameter  vector  a(k)  is  given  by: 


as(k+l)  = ag(k)  + FjXgV^k)  + F2x°v2(k) 


(3.4-1) 


From  Proposition  3.2.1,  recalling  that  u^  = u2  = 0 here,  and  denoting 


Y1(k)  = x°  (k) FjX° (k) 


we  may  express  v^  and  v2  as: 


. 

•: 


(3.4-2) 


r 


(1+Y2)  -y2 


1 + yl  + y2 


-Yx  (l+Yj.)  <Sa2 


(3.4-3) 


Using  (3.4-1)  and  (3.4-3)  we  obtain  the  following  equivalent  nonlinear  des- 
cription of  the  2M-MRAS  dynamics,  where  the  state  here  is  composed  of  x°(k) 

D 

and  ag (k) : 


0 1 . . . 0 


Xg(k+1)  = . 


0 ...  0 1 
aT(k+l) 


xs(k) 


(3.4-4) 


o <WXSXS 

a (k+1)  = [1° 1 ~ a„(k) 

S " 1+Y1+Y2 


(3.4-5) 


T T 

^ (l+Y2)Fl“Y1F2^xsXS  ^ + ^1+Y1^F2~Y2F1^XsXS  aM2 
! + Yi  + y2 


This  parameter  adaptation  equation  will  provide  the  basis  for  analyzing  the 
effect  of  the  quasi-symmetric  pairs  (F^,F2)  and  (a^^  ,a^)  on  the  performance 


of  the  2M-MRAS.  For  the  regulator-type  2M-MRAS,  given  that  stability  is  main- 


tained, the  primary  criterion  of  performance  will  be  the  equilibrium  state 
for  the  2M-MRAS.  Since  we  assume  x°(k)  converges  to  the  origin,  we  are 

* 

interested  only  in  the  dependency  of  the  equilibrium  parameter  vector  a on 

u 

(Fj.V  and  (aM1'aK2>- 

Define  E^(k)  and  E2(k)  as: 


* fl 

l-i 


i r^',4 

!*,  43 
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Ex(k)  = (F1+F2)_1[(1+y2F1-y1F2] 

(3.4-6) 

E2(k)  = (F1+F2)"1[(l+y1)F2-Y2F1] 

(3.4-7) 

The  inverse  exists  since  both  and  F2  are  positive-definite  and  symmetric, 
implying  (F^+F^  is  a^-so*  When  parameter  equilibrium  is  reached,  ag(k+l)  = 
as(k)  = a^,  and  (3.4-5)  becomes  (xg(k)  has  not  necessarily  reached  the 
origin) : 


o °T  * 
[xSxS  aS 


„ o o 

EiVs  V 


T 

, _ o o 
+ E2XSXS  V 


(3.4-8) 


This  equation  alone  does  not  yield  a unique  solution  for  a . In  order  to 

obtain  a unique  solution  in  terms  of  (F^Fj)  and  (aM  ,aM  ),  we  make  the 

”l  2 

following  assumption  about  the  adaptation  gain  matrices: 


F±  = diag  [ f ± ] ; f±  > 0,  scalar  (3.4-9) 

Although  this  restriction  eliminates  both  direct  coupling  among  the  elements 
of  ag(k)  and  weighting  of  one  element  of  as(k)  versus  the  others  in  (3.4-1) 
and  (3.4-5),  this  special  case  of  adaptation  gain  has  been  used  in  earlier 
periods  of  MRAS  theory  development  [32],  and  is  a reasonable  choice  when  a 
sensitivity  analysis  has  not  been  performed  for  the  S-subsystem  to  determine 
how  the  parameters  should  be  weighted  to  minimize  the  state-error  convergence 
rate. 


Given  (3.4-9),  and  E2  become: 


(3.4-10) 
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and  thus  E and  x°x°  commute; 

X b J 


T T 

_ o o o o _ 

EiXSXS  = xsxs  Ei 


(3.4-11) 


Using  (3.4-11),  we  may  rewrite  (3.4-8)  as: 


'Vs  1 (VElaM1'E2aM2)  * 0 


(3.4-12) 


This  equation  must  hold  for  each  stage  k > k , where  k is  the  stage  when 

o o 

ag(k)  = ag.  This  leads  to  the  following  proposition,  which  specified  how  the 
quasi-symmetric  elements  affect  the  S-subsystem  equilibrium. 


Proposition  3.4.1  If  for  the  2M-MRAS  the  following  conditions  hold: 

a)  ui  = u2  = °» 

b)  pole  adaptation  only, 

c)  ag(k)  converges  at  k « kQ  to  a^  corresponding  to  stable  S-sub- 
system  poles, 

d)  x°(ko)  e G+(A*),  where  G+(A*)  is  defined  by  (3.3-23),  then 

* 

ag  satisfies 

. 

as  • q + f2  • 

Proof:  Following  the  development  in  Proposition  3.3.1,  we  claim 

that  [Xg(kQ)  . . . x°(ko+n-l)]  is  nonsingular,  and  thus  (a^-E^a^  -E^^  ) must 
lie  in  the  intersection  of  null  spaces  for  the  matrices  = 

^xS^ko+1)xS  (k0+1)l»  1 = 0 - 1. 

n— 1 

But  H N.  — {0},  so 
1-0  1 

*S  ip-q- 


For  this  special  adaptation  gain  case. 


(3.4-9),  Proposition  3.4.1 


clearly  indicates  the  effects  of  non-symmetry  in  (F^,F2): 


For  f1/f2  = 1» 


lim 

fl/f2  * " 


lim 


Vv  ^ 


* 

a = 


V 


fl/f2 


(3.4-13) 


Given  f^  - f2,  the  effect  of  > M2  on  the  S-subsystem  eigenvalues  corres— 

■k  h 

ponding  to  a is  less  clear.  When  a « au  = aw,  a0  = au.  However,  for 
^ M S H 

> m2»  the  condition 

a + a 

* \ \ 

S 2 


does  not  necessarily  place  the  eigenvalues  of  the  S-subsystem  in  a region 
lying  between  those  of  the  and  M2  subsystems.  The  S-subsystem  eigenvalues 
may  overlap  into  either  the  or  M2  region,  depending  on  the  relative  separa- 
tion of  eigenvalues  within  each  group,  and  respectively. 

The  following  numeric  example  illustrates  the  results  of  Proposi- 
tion 3.4.1  and  compares  theoretical  results  with  results  obtained  from 
simulation. 


Example  3.4.1 


% 


nd  T 

Consider  the  2 — border  regulator-type  MRAS  with  a^  - [-.81,  1.8], 

[-.25,  1],  corresponding  to  poles  at  (0.9,  0.9)  and  (0.5,  0.5),  re- 


spectively. The  S-subsystem  is  initialized  at  a*(0) 

u 


[-.01,  .2],  corres- 
ponding to  poles  at  (0.1,  0.1).  The  following  results  are  obtained: 
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II 


I 

E 

E 


E 


fl 

f2 

* 

as~theoretical 

* 

Sg-simulated 

poles-theoretical 

poles-simulated 

100 

100 

[-.53,  1.4] 

[-.516,  1.38] 

0.7  ± j .2 

.69  ± j .199 

1000c 

100 

[-.8044,1.792] 

[-.7276,1.706 

.896  ± j .04 

.853  ± j .006 

100 

10000 

[-.2555,1.008] 

[-.309,1.09 

.504  ± j .0388 

.546  ± j .107 

When  f =*  the  theoretical  and  simulated  results  agree  quite  well.  However, 
in  the  other  two  cases,  > F2  and  F^  < F^,  the  simulated  poles  do  not  move 
all  the  way  out  (F^  > F^)  or  all  the  way  in  (F^  < F^)  to  the  theoretical  pole 
values.  It  has  been  found  that  as  and  f2  are  both  increased,  it  is  pos- 
sible to  reduce  this  error.  Thus,  we  conclude  that  this  error  between  the 
simulated  and  theoretical  values  is  due  to  the  regulator  property  of  the  S- 
subsystem;  adaptation  is  diminished  in  proportion  to  how  close  x°  is  to  its 
equilibrium  point  at  the  origin.  Large  f^  and  f2  cause  more  of  the  adapta- 
tion to  occur  early  in  the  trajectory  of  x°  before  it  becomes  too  close  to 
zero,  thus  improving  the  comparison  between  theoretical  and  simulated  results. 

Future  analysis  of  the  2M-MRAS  should  consider  the  more  general  case 

T 

of  adaptation  gains  F^  >*  F^  > 0.  In  addition,  the  cases  where  u^  and  are 
non-zero  and  where  numerator  adaptation  is  permitted  should  be  studied. 
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CHAPTER  4 
C01ICLUSI01JS 

4.1  Review 

In  this  thesis  several  topics  involving  a parallel  MRAS  have  been 
investigated.  We  briefly  summarize  here  the  main  results. 

For  the  discrete-time  single-input  single-output  MRAS,  we  introduced 
the  Input-Output  Delay  state-variable  representation  for  the  M-  and  S-sub- 
systems.  Controllability  conditions  for  this  representation  expressed  in 
terms  of  the  equivalent  transfer  function  were  found  and  used  to  produce  an 
M-subsystem  whose  order  is  less  than  that  of  the  S-subsystem.  This  was  done 
in  such  a way  that  the  Landau  adaptation  algorithm  remains  compatible  with  the 
reduced  order  subsystem;  thus  the  MRAS  retains  its  hyperstable  property. 

For  the  MRAS  applied  to  parameter  identification,  we  considered  two 
separate  problems.  The  first  dealt  with  asymptotic  stability  of  the  parameter- 
error,  the  other  with  hyper stable  design  realization  when  the  M-subsystem 
(plant)  parameters  are  unknown. 

The  asymptotic  stability  of  parameter-error  is  not  inherent  in  the 
hyperstable  design.  We  obtained  several  algebraic  and  geometric  representa- 
tions for  the  parameter-error  equilibrium  subspace.  From  these  representations 
we  were  able  to  state  necessary  and  sufficient  conditions  for  asymptotic  sta- 
bility of  the  parameter-error,  which  means  that  the  parameter-error  equili- 
brium subspace  contains  only  the  origin.  The  practical  use  of  this  result  is 
limited  by  the  fact  that  the  conditions  depend  on  the  unknown  plant  parameters. 

The  Landau  adaptation  design  guaranteeing  hyperstability  is  not 


I 


1 


directly  realizable  because  the  M-subsystem  (plant)  parameters  are  unknown. 
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We  proposed  a modified  adaptation  design  which  satisfies  the  positive-real 
hyperstability  conditions  asymptotically,  assuming  the  MRAS  remains  stable. 
We  suggested  a strategy  to  investigate  the  stability  properties  of  this 
modified  algorithm.  This  investigation  might  find  useful  a test  developed 
by  Siljak  for  positive-realness  on  the  unit  circle. 

Using  results  from  the  theory  of  variable  structure  systems  (VSS), 


we  have  established  the  stability  of  a particular  continuous-time  MRAS 
structure  which  uses  discontinuous  adaptation.  This  MRAS  was  first  shown 
to  be  equivalent  to  a general  VSS  which  is  designed  to  reject  measurable 
disturbances,  and  then  stability  results  for  VSS  were  applied  to  the  MRAS. 

A new  MRAS  structure,  one  with  multiple  reference  models,  was 
introduced.  Motivation  for  such  a structure  is  found  from  dynamic  game 
theory.  We  examine  the  stability  properties  for  the  two-model  MRAS 
(2M-MRAS),  given  that  the  two  adaptation  algorithms  are  of  the  hyperstable 
form  (with  respect  to  the  one-model  MRAS).  A 2M-MRAS  and  a dynamic  game 
have  been  jointly  developed.  Alternate  viewpoints  in  obtaining  the  M^- 
subsystems  of  the  2M-MRAS  have  been  described.  A procedure  for  comparing 
the  2M-MRAS  parameter  equilibrium  with  the  strategy-dependent  dynamic  game 
solution  has  been  developed  and  illustrated  by  comparing  a particular 
2M-MRAS  equilibrium  with  the  closed-loop  Nash  solution  of  the  dynamic  game. 
We  have  also  obtained  the  2M-MRAS  regulator  parameter-equilibrium  as  a 
function  of  M^-subsystem  and  A^-subsystem  parameters,  for  a particular 


adaptation  gain  structure. 
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4.2  Suggestions  for  Further  Research 

The  result  in  Section  2.3,  pertaining  to  the  reduction  In  M-subsystem 
transfer  function  order,  actually  maintained  the  same  system  order  for  both  the 
M-  and  S-subsystems  from  a state-variable  viewpoint.  However,  one  could 
construct  a MRAS  where  the  order-reduction  for  either  the  M-  or  S-subsystem  is 
explicit  in  the  state-variable  formulation;  that  is,  the  difference  in  order  is 
not  achieved  through  cancellable  transfer  function  poles  and  zeroes.  To  apply 
Landau's  adaptation  algorithm  in  this  context,  it  would  be  necessary,  first  of 
all,  to  define  a new  state-error  vector.  The  hyperstability  of  this  MRAS 
formulation  should  be  studied.  Preliminary  investigations  suggest  that  the 
difficulties  encountered  in  the  stability  analysis  for  this  system  are  similar 
to  some  of  the  problems  encountered  in  the  analysis  of  the  2M-MRAS  in 
Section  3.2. 

The  questions  developed  in  Section  2.5  concerning  the  hyperstability 
of  a modified  form  of  the  Landau  MRAS,  where  the  design  vector  c is  chosen 
to  be  time  varying: 

c(k)  - -a(k)  (4.2-1) 


should  be  investigated  further.  Use  of  the  Siljak  test  during  simulations 
would  be  helpful  in  determining  whether  it  is  necessary  at  every  stage  k to 
satisfy  the  Popov  positive-real  condition. 

The  results  in  Chapter  3 dealing  with  the  multi-model  MRAS  provide  a 
basis  from  which  generalizations  could  be  made  in  several  directions.  The 
S-subsystem  equilibrium  parameter  vector  for  the  regulator-type  2M-MRAS  could 
be  investigated,  assuming  more  general  forms  for  the  adaptation  gain  matrices. 
Further  generalization  could  be  attempted  by  assuming  nonzero  external  inputs 
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(e.g.,  two  different  constant  set-points),  and  considering  full  (both  poles 
and  zeroes)  parameter  vector  adaptation. 
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APPENDIX  A 

In  this  appendix  we  present  some  matrix  notation  and  definitions, 
and  then  develop  some  algebraic  propositions  related  to  these  matrices.  The 
material  here  supports  the  developments  of  Chapter  2,  serving  primarily  as  a 
central  repository  for  matrix  notation,  allowing  the  main  text  to  flow  more 
smoothly.  The  propositions  are  also  of  some  interest  in  themselves  in  rela- 
tion to  the  theory  of  linear  algebra,  and  in  the  Insight  gained  into  the 
structure  of  certain  matrix  forms. 


A.l  Algebraic  Structure  of  Two  Finite  Matrix  Sets 

In  this  section  we  present  our  notation  for  two  finite  matrix  sets, 

and  then  describe  their  algebraic  properties  under  the  binary  operation, 

matrix  multiplication.  The  first  set,  denoted  1N,  consists  of  those  N x N 

matrices  with  l's  on  a single  diagonal,  and  0's  elsewhere.  The  identity 

matrix,  1°  , is  a familiar  element  of  this  set.  The  second  set,  denoted  P,  , 

N N 

consists  of  those  N x N matrices  with  l's  on  a single  anti-diagonal,  and  zeroes 
elsewhere. 

Definition  A. 1-1 
XN  = £ lnte8er8*» 

k 1 
l£(k,j)  - { 

0 

Definition  A. 1-2 

PN  * e Integers),  where  matrix  elements  of  P^  satisfy 


jj^ 

where  the  matrix  elements  of  Ijj  satisfy: 

if  (j-i)  = k 
otherwise 


k 1 if  (N+l)-(i+j)  = k 

PN(i.J)  = { 

0 otherwise 


If  k > 0,  l's  appear  on  a superdiagonal  or  anti-superdiagonal. 

If  k < 0,  l's  appear  on  a subdiagonal  or  anti-subdiagonal. 

If  k = 0,  l's  appear  on  the  main  diagonal  or  main  anti-diagonal. 

The  two  sets,  1^  and  P^,  satisfy  the  following  algebraic  properties: 

Property.  1:  ijj  - 0 and  p|j  = 0 for  all  |k|  > N.  This  follows  since 

both  1,  j e [1,N];  thus: 


rr 

i 

i 

i 

i i 

1 1 

i 

1 1 
1 1 
i 

i 

- 

i 

r. 

E 


D 

L 

G 


c)  P°  Pk  - l""k 
* N N Tj 

d)  PN  PS  “ 

e)  PS  *N  PS  = ^ 

f)  P°  Pk  P°  = P~k 
' N N N N 

kl  k2  kl  "k2 

*>  p„  pN  -V1, 

h)  Pkl  ,k2  . ^ p-k2 

J N N N 

k2  kl  k2  ~kl 

i)  L‘iPi  = PZI  i 

J N N N N 


Property  3: 

kl  k?  (ki+kJ 

a)  V V = \ if  sgn(ki*k2)  = +1 


kl  k9  k9  k1  <k  +k  ) 

b)  hhf  + hh'-h 


Property  4: 

kl  k2  (kl'k2) 

a)  P 1 p • = I 1 1 

N N N 


if  sgn^*^) 


-1 


if  sgn  (k^k^  “ -1 
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kl  k2  ~k2  _kl  (ki -k9^ 

b)  PN  PN  + PN  PN  " h if  Sgn  (kl*k2>  = +1 

This  follows  from  properties  2(g)  and  3. 

It 

Properties  1 and  3(a)  imply  that  1^  has  a nilpotency  index  [10] 
equal  to  [N/k],  where  the  notation 

Min 

* - [N/k]  = In  | n > N/k> 

Integers 


The  elements  of  1^  and  P^  occur  in  the  development  of  Chapter  2 as 
either  right-  or  left-hand  linear  operators  on  general  matrices  or  vectors  of 
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compatible  dimension.  We  include  a qualitative  description  of  their  effects 

on  general  matrices  here  to  clarify  that  development. 

. _NxN 
Given  A R , 

— K 

1)  AI^  shifts  columns  of  A right  (+)  or  left  (-)  k steps,  filling 
with  0-columns. 

2)  I"  A shifts  rows  of  A up  (+)  or  down  (-)  k steps,  filling  with 

0-rows. 

ik 

3)  APn  pivots  A about  its  central  vertical  axis,  then  performs 


4)  P~  A pivots  A about  its  central  horizontal  axis,  then  performs 

IfA* 


SN  = ^ PN  {0},  and  the  binary  operation,  matrix  multiplication, 
denoted  *,  then  {S^,*}  does  not  form  a binary  algebra  or  semigroup  [5].  Prop- 
erties 3 and  4 show  that  • is  not  a closed  operation;  i.e,,  it  generates  ele- 
ments which  are  not  members  of  S . 

N 


A. 2 Matrix  Relationships  Involving  the  Companion  Form 

The  companion  form  for  a matrix  appears  often  in  problems  dealing 
with  state-variable  realizations  of  a LTI  system  [21,42],  and  is  also  a funda- 
mental form  in  studying  the  general  theory  of  linear  transformations  [10],  In 
this  section  we  develop  some  notation  and  relations  which  involve  the  compan- 
ion form,  and  which  also  appear  in  Chapter  2. 
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A. 2.1  The  Output  Delay  System 

Consider  the  difference  equation: 


y(k)  = E a y(k-i) 
i“l  1 


(A. 2-1) 


A state-variable  realization  which  occurs  naturally  from  this  equation  is: 


w(k+l)  - F w(k) 


y(k)  = a w(k) 


where  w (k)  » [y(k-N).  . .y(k-l)] 
aT  = [aN  ‘ * al] 


0 1 ...  0 


(A. 2-2) 


(A. 2-3) 


(A. 2-4) 
(A. 2-5) 


0 ...  1 


+ 0 


“M  * 


(A. 2-6) 


I 


The  dynamic  system  matrix  is  known  as  the  companion  or  1—  natural  form  [10], 
and  is  a fundamental  canonic  matrix  in  the  study  of  linear  transformations. 

We  will  refer  to  the  state  realization  w(k)  as  the  Output  Delay  State 
Realization. 

The  next  proposition  provides  the  basis  for  stating  the  necessary 
and  sufficient  conditions  for  observability  of  the  pair  (aT,F).  It  also 
provides  insight  into  the  algebraic  structure  of  the  companion  matrix  F. 

Proposition  A. 2.1 

a)  F has  an  eigenvalue  \ ■ 0 of  multiplicity  (N-k)  and 
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b)  p{F  } = k for  all  a £ (N-k)  if  and  only  if  the  elements  of 
T 

a = [a^  . . . a^]  satisfy: 

^aN*  * * * \+l^  “ \ ^ • • • a^)  arbitrary 

where  k may  take  the  value  k £ [0,N],  and  p{ • } denotes  the 
matrix  rank. 

Proof : Since  F is  cyclic  [10],  the  characteristic  and  mini- 

mal polynomials,  A(A)  and  iJ>(A),  respectively,  are  equal: 

A(A)  = iKA)  = AN_k(Ak  - l a.Ak_i) 

i=l  1 

Thus  F has  an  eigenvalue  at  A = 0 of  multiplicity  (N-k). 

The  degeneracy  of  F = 1,  for  A = 0,  since  the  Jordan  form  of  F has 
only  one  Jordan  block  associated  with  each  eigenvalue.  This  fact  follows  from 
the  cyclicity  of  F.  This  degeneracy  property  may  be  expressed: 

dim{W[F-AI^]  | A = 0}  = dim{W[F]}  = 1 

where  W[»]  = null  space  of  a matrix. 

The  generalized  eigenspace  [13],  W[F-Al°] ^N_k) , associated  with  A = 

0 has  the  following  property: 

dim{W[F-Al°]N_k  | A = 0}  = dim{M[F]N_k}  -N-k, 

and  further,  for  all  a > N - k: 


dimWfF-Al®]01  | A - 0}  - dim{W[F]a)  -N-k 
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Thus,  p{[F-Al]“  | A - 0}  - plF0}  » k for  all  a > N - k.  This  is  true  since 
p{Fa}  + dim{N[F]a}  » dim{P[F]a}  = N, 
where  D[»  ] denotes  the  domain  of  a matrix. 

p{F^}  = k for  all  a - N - k implies 


dim{W[F-Al°]a  | X « 0}  = N - k for  all  a > N - k. 


This  in  turn  implies  X m 0 is  an  eigenvalue  of  F with  multiplicity  (N-k) , 
and 

A (A)  =*  AN"k(Ak  - E a. Ak_i] . 

i=l  1 

Thus,  the  characteristic  equation  requires  (a  , ....  a.  .)  = 0,  a j*  0. 

N Kt X R | 

lr 

By  enumeration  of  F , k ■ (1,2,.  . . ,N) , it  is  possible  to  demon- 
strate that: 


e • • 

T 

a 

aT  F 
aT  F^1 


+ * 


(A. 2-7) 


When  k * N we  have: 


--  **?■*!*■  ■■  >-9SfsS*t WfJJfyW' 
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T 

a F 


;t  f*-1 


(A. 2-8) 


since  = 0.  Recognizing  that  F^  is  the  observability  matrix  for  the  pair 
(a1,!7),  we  have  the  following  proposition: 


Proposition  A. 2-2 

The  system  (A. 2-2,  A. 2-3)  is  completely  observable  if  and  only  if 


Proof:  It  is  sufficient  to  show  that  is  nonsingular  if  and  only 

if  &U  * 0. 

— ~V  F1*  nonsingular  implies  F**  has  no  eigenvalues  X ■ 0.  In  the 
context  of  Proposition  A. 2-1,  k » N,  and  we  have:  F has  no  eigenvalues 

X = 0,  and  pfF01}  = N for  all  a £ 0 imply  a^  J4  0. 

■■  ■■  - ajj  t 0 implies  the  case  in  Proposition  A. 2-1  where  k = N.  Thus 

a^  l*  0 implies  F has  no  eigenvalues  at  X = 0 and  p{Fa}  * N for  all  a £ 0.  In 

particular,  p{FN}  = N and  FN  is  nonsingular. 

A. 2. 2 A Matrix  Sequence  Induced  by  the  Input-Output  Delay  System 

The  2N  x 2N  matrix,  F,  defining  the  dynamics  of  the  1-0  Delay  state 
realization  (2.2-8),  induces  the  matrix  sequence  F^,  k * (0,1,2,.  . .),  as  a 
result  of  the  solution  to  the  1-0  Delay  state  difference  equation.  Imbedded 
in  F^  appears  the  following  N x N matrix  sequence,  which  we  denote  F(k): 


where  b **[bN  . . . b^  (A.  2-10) 

* /v 

Note  that  F(o)  is  undefined,  but  we  will  adopt  the  notation  F(o)  = 0 for 


),  it  is  possible  to  show 


Alternatively,  we 


F(k+1) 


Of  particular  importance  to  the  developments  of  Chapter  2 is  the 

A 

sequence  element  F(N): 


bTI^  + aTF(0) 

bTI^  + aTF(l) 


b1!^-1  + aTF(N-l 


Using  the  F(k)  notation,  we  may  express  the  sequence  FK  as 


(A. 2-14) 


In  particular,  for  k = N, 


0 I 0 

I 


(A. 2-15) 


The  observability  matrix  for  the  1-0  Delay  state  realization  may 
also  be  expressed  in  a very  compact  form: 

I bTl£  + aTF(0) 

PT  aT  1 bTI^  + aTF(l) 

P1  F aTF  ] . 


(P1)? 

(pW 


.v-1 


bTl^-1  + aTF (N-l) 

bTI°  + aTF(071 

b1!^  + aTF(l)  * 


i ; 


v-i 


bTlJ!_1  + aTF(N-l) 


..  . ,t ....  ...i< 


- * 

■ 


Thus,  using  (A. 2-15) 


(A. 2-16) 


l 
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